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HAUPTAUFSÄTZE 


Fine Integralgleichung für den rotierenden Schaufelstern 
und ihre Lösung‘). 


Von W. Kucharski in Berlin. 

Aufstellung der Integralgleichung für die Wirbelverteilung auf den Schaufeln; 
Transformation auf eine in der Tragflügeltheorie bekannte Form. Lösung 
für die Schaufelzahlen n=# nach der bekannten Methode durch Vergleich 
der Fourierkoeffizienten. Entwicklung eines anderen Verfahrens fürn >4. 
wo eine integrable Unstetigkeit auftritt: Abspaltung und geschlossene Er- 
ledigung der unendlich werdenden Funktion, Behandlung der beschränkten 
und stetigen Restfunktion nach der üblichen Methode. Näherungen für schwach- 
gekrümmte Schaufeln. Schaufeln,. die nicht bis in den Nullpunkt gehen. 
Teilweise Umkehrung des Problems. 


1. Einleitung. Das ebene Strömungsproblem des in idealer Flüssigkeit mit der kon- 
stanten Winkelgeschwindigkeit & rotierenden Schaufelsterns, der aus n geradlinigen, vom 


> 
: j r Ki BR: > . 
Nullpunkt ausgehenden, jeweils den Zentriwinkel es miteinander bildenden Platten von der 


radialen Länge I besteht, ist als Randwertaufgabe der Potentialtheorie mehrfach behandelt 
worden, insbesondere von W.Spannhake und Sörensen, nachdem ich eine erste behelfs- 
mäßige, angenäherte Lösung angegeben hatte. Ist y die Stromfunktion und x die Entfernung 
eines Punktes auf den Schaufelplatten vom Nullpunkt, so handelt es sich um Ermittlung einer 


Lösung von Ay=0 mit der Randbedingung y' — - .. x? auf den Schaufeln, da diese die 
Normalgeschwindigkeit »„—=w-x haben. Die Lösung ist bestimmt bis auf eine beliebige 
Zirkulationsströmung, die dann in üblicher Weise so gewählt wird, daß man auf den äußeren 
Schaufelenden tangentiales Abströmen erhält. Der Wert dieser Lösungen liegt darin, daß 
der radiale Schaufelstern als einfachstes, wenn auch sehr rohes Schema für eine Radialpumpe 


1) Der Inhalt der folgenden Abhandlung mit Ausnahme der Abschnitte 12 bis 14 wurde Februar 198 im Seminar 
für Mechanik an der Technischen Hochschule Berlin von Herrn Dipl.-Ing. A. Frieke vorgetragen, der an dem Zu- 
standekommen der Arbeit beteiligt ist. Herr Frieke hat die Lösung für » <4 in den Einzelheiten ausgearbeitet, einen 
Teil der übrigen Rechnungen kontrolliert und selbst durchgeführt sowie für die in Abschnitt 6 eingeführte Funktion 4 
die geschlossene Darstellung aufgefunden, wodurch die Auswertung meiner Ergebnisse wesentlich erleichtert wurde. 
Auch für manche Politur in mathematischer Hinsieht bin ich Herrn Frieke dankbar. Der Verfasser. 
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oder -turbine anzusehen ist. Den genannten Verfassern ist es gelungen, auch solche Fälle 
zu erfassen, bei denen die Schaufeln nicht bis in den Nullpunkt gehen und außerdem nach 
logarithmischen Spiralen gekrümmt sind, was eine bedeutende Verbesserung des Schemas mit 
Rücksicht auf die wirklich vorkommenden Verhältnisse bedeutet. 

Im folgenden gehe ich zunächst von geradlinigen Schaufeln aus, die sich bis in den 
Nullpunkt erstrecken. Die Lösung wird hier aber nicht nach den Randwertmethoden der 
Potentialtheorie gegeben, sondern durch Aufstellung und Lösung der Integralgleichung für 
die Wirbelbelegung auf den Schaufeln. 

Dies hat für die Auswertung manche Vorteile, denn bekanntlich kann man aus der 
Wirbelverteilung sofort die Geschwindigkeits- und Druckunterschiede zu beiden Seiten einer 
Schaufel hinschreiben, ohne auf die sonstigen Einzelheiten der Strömung einzugehen. Auch 
ist die Methode der Wirbelbelegungen für die näherungsweise Behandlung schwach ge- 
krümmter, von einem Radius nicht zu stark abweichender Schaufelformen von innerhalb 
dieses Rahmens weitgehend beliebigem Krümmungsverlauf brauchbar, worauf in Abschnitt 12 
dieser Arbeit unter sinngemäßer Umgestaltung eines in der Tragflügeltheorie üblichen Ver- 
fahrens eingegangen wird. 

Dann wird kurz der Schaufelstern mit nicht bis zum Mittelpunkt gehenden Schaufeln 
behandelt, wofür die Lösung nach der einen der entwickelten Methoden ebenfalls möglich 
ist; auch die Behaudlung schwach gekrümmter Schaufeln, die nicht bis zum Mittelpunkt 
gehen, ist mit den angegebenen Mitteln in der erwähnten ersten Annäherung durchführbar. 
Hieran schließt sich noch eine Erörterung für die teilweise Umkehrung des Problems, indem 
diejenige von der Geraden abweichende Schaufelform gesucht wird, welche zur Erzielung 
einer gewünschten Abweichung von der für die geraden Schaufeln gefundenen Wirbelver- 
teilung erforderlich ist. 

Es wird sich zeigen, daß bei n=4 die Lösung der Integralgleichung ohne jede Kom- 
plikation mit bekannten Mitteln möglich ist. Bei » >4 würde aber die bekannte Methode 
die Fourierentwicklung einer unendlich werdenden absolut integrablen Funktion verlangen. 
Um dies zu vermeiden, habe ich einen Lösungsweg entwickelt, bei welchem die unendlich 
werdende Funktion für sich in geschlossener. Form erfaßt wird, wobei für eine Fourier- 
entwicklung nur noch eme beschränkte stetige Funktion übrigbleibt. 

2. Aufstellung der Integralgleichung. Zur Aufstellung der Integralgleichung seien zu- 
nächst die Normalgeschwindigkeiten v, bestimmt, die von einem System von n konzentrierten 
Wirbeln mit der Einzelintensität m,, von denen jeder im gleichen Abstand & vom Nullpunkt 
auf einer Schaufel liegt, auf einer der Schaufeln „erzeugt“ oder „induziert“ werden; diese 
liege in Richtung der reellen Achse in der w-Ebene; es sei w=.r.+iy. 

Ein konzentrierter Einzelwirbel, der in einer 2,-Ebene im Abstand &” vom Nullpunkt 
auf der reellen Achse liegt, hat das komplexe Potential: 





m i 
fe) =; log (2, — &"). 
Der Übergang zur w-Ebene mit 
n 
w=y2, 
liefert dann das komplexe Potential F(m) der genannten n symmetrisch auf dem Kreis vom 
Radius £ liegenden Einzelwirbel von der gleichen Einzelintensität m, zu: 


mM, 


F(iw) = 


log (m" en) 


Die komplexe Geschwindigkeit dieser Wirbelanordnung ist: 


; n m, nn ww"! 
rd. —tid,y,= (m) ae =; 
- I w & 


Für die Punkte der reellen Achse folgt hieraus: 
N» +3 


"u vy—=M,' re n* 


ir Ir 


Jetzt werde 5 als Variable im Bereich O=5=1! auf der x-Achse aufgefaßt; ferneı 
werde an Stelle des einen konzentrierten Wirbels eine Wirbelbelegung m (£) auf den Schau- 
feln aufgebracht, derart, daß m (£).d& die Wirbelintensität auf der Länge d£ an der durch & 
gegebenen Stelle ist. Die Normalgeschwindigkeit an der Stelle x ist jetzt die Summe aller 
von den Intensitäten m (£)-d& herrührenden Anteile, d. h. es ist jetzt: 
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Dabei ist für das unendlich werdende Integral, wie überall in diesem Aufsatz, der Cauchysche 
Hauptwert zu nehmen. 

Das Schaufelsystem rotiert mit der Winkelgeschwindigkeit &; jede Schaufel hat also 
im Abstande x die Normalgeschwindigkeit x. Die Aufgabe wird hier nun so gestellt und 
gelöst, daß die Randbedingung v„—= © x, welcher die Gesamtströmung genügen muß, durch 
die von den Wirbelbelegungen m (£) auf den Schaufeln herrührenden Normalgeschwindig- 
keiten vollständig erfüllt werden Soll; zum Erhalt einer „Pumpen- oder Turbinenströmung* 
ist es dann nur noch notwendig, eine Quell- oder Senkenströmung von oder nach dem Null- 
punkt hinzuzufügen, wodurch die Randbedingung nicht berührt wird. M.a. W.: Die Normal- 
geschwindigkeit nach (1) soll gleich & x sein. 

Mit 

I 


on 7) 1 


a & EEE 


0 

erhält man so die für das vorliegende Problem zu lösende Integralgleichung; sie ist von 
erster Art und singulär. 

Die hieraus zu bestimmende Funktion m (£) muß außerdem die Bedingung tangentialen 
Abströmens an den Schaufelenden erfüllen; bekanntlich ist hierzu erforderlich: 

BE et ae oe 

Selbstverständlich ist: m =0 für E >1. 

Aus den erwähnten bisherigen Behandlungen dieser Strömung weiß man, daß das 
Problem durch die besondere Bedingung (3) eindeutig gemacht wird. 

3. Transformation in eine bekannte Gestalt. Durch einige Umformungen werde zunächst 
Gl. (2) auf eine in der Tragflügeltheorie bekannte Form gebracht, die z. B. in: Hamel, „Integral- 
gleichungen, Einführung in Lehre und Gebrauch“, Berlin 1937, S. 145 usw. behandelt wird. 


Mit 
re" Y: 2 
m|yz) (4) 
M (z) 
n—1 
2 n 
geht (2) über in die Integralgleichung für M (2): 
In 
"M (z) 0) ® 
\ dz (5). 
h Yy zZ n 
0 y n 
Hierin werde noch gesetzt: 
jr Ir 
z 5; 9% #3 
(6). 
M(w+ . |= Fi) 
womit (5) übergeht in: 
In 
£ £: F (m) , 
Far Im — fu) (), 
in 
worin: 2 
fi: Tr w (5) 


n 
2 alu + 5) n 
\ - 


eine bekannte, vorgeschriebene Funktion von « ist. Gl. (7) ist die verlangte Form. (Es sei 
noch betont, daß die hier eingeführten Größen y,z,w mit den zur Ableitung von (1) be- 
nutzten mit gleichlautenden Bezeichnungen nichts zu tun haben; Mißverständnisse sind wohl 
ausgeschlossen.) 


4. Lösung für n<4. Für n<4 ist nun die Lösung von (7) nach der bekannten Methode 
ohne Komplikationen zu gewinnen. Mit den neuen Variabeln y und 9, so daß: 
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1" ni 1" 
w=—z:.008y; u=—7 0089; 
1" 4 
F(- y cosy)- sin y— Yıy) Zt Er © & 
ferner mit a‘ ! 
23,9 
geht (7) über in: 
Y(y) a 
en N ER IST 
0 sin * oj2 E 
Zur Gewinnung von Y(y) setzen wir zunächst formal die Fourierentwicklung an: E 
6 & 
. 1 ” N 
Yy)=z>b+ \ BO a Sea 5 ee 
je _ı : 
tragen dies in (10) ein, berücksichtigen, daß nach einer bekannten Beziehung, die auch für j 
v—=ß0 gilt: ' 
n \ 
ll cosry 1 sinv« ? 
lea ir --7 a) Ah 
zn )coSp — COS y z sıng “ 
0 5 
und erhalten: 2 
y 2a-sing i i 
er ame Me 
2 P} |) ® 
‚—1i sın “ q j= 
’ i Lo > 
Mit Ausnahme der zunächst noch unbestimmt bleibenden Konstanten — b, sind also die | 
Fourierkoeffizienten der Cosinusentwicklung von Y(y) diejenigen der Sinusentwicklung der 
in (13) rechts stehenden Funktion, und die zunächst formal für Y(y) angesetzte Entwicklung 7 
geht in Ordnung, wenn diejenige nach (13) möglich ist, also im einfachsten Falle sicherlich £ 
dann, wenn g 
sing un 2n—4 s 
fig = sinkyR BB ie nn en 
im Bereich Null bis x einschließlich der Randpunkte stetig ist. Unstetigkeit kann nur für 
y=®0 eintreten. Von konstanten Faktoren abgesehen, verhält sich f(g) beim Nullpunkt wie i 
y'"K; dies ergibt: 2 
lim f(g)=0 für K<1 £ 
y>V 
=2 für K= (15). 
—»fürK >1 | 
Die Reihen (11) und (13) konvergieren also und stellen die Funktion im ganzen Intervall 
dar, wenn: 
A a re Ale 2 0 TS i 
wenn also 
a Er A a 1 a et Fe Ei 





Für diese Fälle erhalten wir also auf dem eingeschlagenen Wege ohne jede Komplikation 
die Lösung, d. h. für die Schaufelzahlen 


WE ES ea ie a Ar 
da nichtganze Zahlen, die eine vielfache, sich nicht schließende Überdecekung der Ebene be- 
deuten würden, physikalisch für das vorliegende Problem ohne Bedeutung sind. 

Unter Voraussetzung von (17) wird in bekannter Weise: 


a 
' 4a \ sing 
R-- . = — 
g sinKk @]2 
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Geht man nun wieder von Y(y) auf F(w) und schließlich auf m (£) zurück, so erhält 
man nach einiger Rechnung: 
’. 


=, 


n eN-ı / 1 
sb+ b„.cosvy), 
ur |— / 


m(d)= BD) "yıls)" ze | 

wobei mie ee a 

€ n 

cos y = 1—2 (5) 

Die noch unbestimmte Konstante b, dient zur Erfüllung der Randbedingung (3) für das 

tangentiale Abströmen am Schaufelende, wonach m(l)=0. Da der Nenner in (20) für £>1 

nach Null geht, muß zur Erfüllung von (3) auch der Zähler für £>1, also cosy—>--1, Null 
werden; hieraus folgt: 

1 
a b 


ee 1 
v—1 


Die endgültige Lösung für n= 4 ist also: 


2 


n en/2-ı \ ” 
m()=7z' u b,[ecosy —(-1N”] . . . . . . (22) 
) - A 


—1 
mit 
cosy—=1— 2 (£/1" 
und 5, nach (19). 
In den Fällen » = 1,2, 4 sind die b, von Herrn Frieke in geschlossener Form berechnet. 
Man erhält mit 





Pe de a Er ar 
am 1+6—20% 
—- +3 == für n=1; 
lo Vs C 
2am : 
— 91/1 _ PR rm >. 
g 2yl-i für n=2; 
x (24). 
2m 4 E Tr 4r R 
Fa er en für m 4; 
vi 
mit cosy=1--2{* 


In der Gleichung für n=4 ist nach „Knopp, Theorie und Anwendung der unendlichen 
Reihen“, Berlin 1931, Kapitel VIII, Gebrauch gemacht von der Summierungsformel: 


5:’_#r Eat se 
Lu 47° wu "En h. 
v—1 

5. Gedankengang einer anderen Lösung für n >4; Ansatz. Bei n >4, also für die 
praktisch vielfach interessanteren Fälle, wird die in Gl. (13) auf der rechten Seite stehende 
Funktion im Nullpunkt unendlich, hier also sicherlich nicht durch die sin-Entwickelung dar- 
gestellt. Ich sah daher 1937/38 bei Ausarbeitung dieser Abhandlung die Verwendung der 
bekannten Methode bei n >4 als zum mindesten zweifelhaft an; dies führte mich zur Ent- 
wickelung der jetzt darzulegenden Methode für e>n>4 oder -1<u<—1/2, wobei also 
die maßgebende Funktion im Nullpunkt unendlich wird, dabei aber absolut integrabel bleibt, 
wie man aus den an Gl. (13) anschließenden Überlegungen sofort ersieht. 

Nach einer kürzlich erschienenen Arbeit, auf die mich Herr Hamel während der Druck- 
vorbereitung des vorliegenden Aufsatzes freundlicherweise hinwies (Weissinger, Ein Satz 
über Fourierreihen und seine Anwendung. in der Tragflügeltheorie, Math. Z. Bd. 47 (1940), 
S. 16 bis 32, Heft 1), ist heute die Zulässigkeit des üblichen Verfahrens auch bei endlich 
vielen Unendlichkeitsstellen sichergestellt, wenn die Funktion absolut integrabel ist; meine 
früheren grundsätzlichen Bedenken waren also nicht berechtigt. Es dürfte aber nach wie vor 
die praktische Schwierigkeit verbleiben, daß die Fourierreihe für die Nachbarschaft der 
Unendlichkeitsstelle sehr viele Koeffizienten erforderlich macht, wenn diese auch nur einiger- 
maßen richtig erfaßt werden soll. Dieser Übelstand wird durch die: folgende Methode ver- 
mieden, die mir auch mathematisch einige bemerkenswerte Züge infolge der eigentümlichen 
dabei auftretenden Funktionen aufzuweisen scheint. 
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Zunächst werden zur Vereinfachung der Schreibweise neue Variable und Bezeichnungen 
eingeführt; 


y=s:l®"; =t-1*; | 
M(t-")=N(t); 
BR w AP 
o ER | 
A, 1" 2» re n 
womit (5) übergeht in: 
) 
'Ntt 
\ wie 
\s—t 
0 


Die gesuchte Funktion X (f) wird nun in zwei Teile zerlegt: 

NO=N,Md+N,() Me (27); 
dabei soll der Anteil N,(f) so bestimmt werden, daß er nach Einsetzen von (27) in die 
Integralgleichung (26) die auf der rechten Seite stehende, zur Unstetigkeit führende Potenz 
zum Verschwinden bringt und nur noch beschränkte und stetige Restbeträge stehen läßt. 
Auf diese Weise entsteht eine neue Integralgleichung für N,(f), die nun nach dem bekannten 
Verfahren der vorigen Abschnitte ohne jede Schwierigkeit lösbar ist. Ferner soll X, (f) die 
Bedingung 


Nd)=V0 TE 
erfüllen, und da selbstverständlich nach wie vor 

X\V=0 
sein muß, so folgt auch die Forderung 

ee... . Re ee Br En 0 26 


N,(f) kann dann also sicherlich nach der für » =4 angegebenen Methode bestimmt 
werden, die ja gerade für die Randbedingung (29) brauchbar ist. (Außerdem für m()=», 
was aber bei dem hier behandelten Problem aus physikalischen Gründen ausscheidet). 

Zu dem beabsichtigten Ziel führt der folgende Ansatz: 

rc tn. 54 3 
mit noch zu bestimmenden Konstanten e und o. Man sieht zunächst, daß er die Bedingung 
N, (1)=0 erfüllt; das Weitere ergibt sich folgendermaßen: 

6. Eine Integralbeziehung. Vorerst ist das Integral 

1 
"gu fo 
Ism0-\, a EN > 
0 


in geeigneter Weise darzustellen®). Ich gehe aus von 


1 reihe 1 
u [ fr" f" 
dt= dti— dt, 
s—t t 
s(i Ai 
s 
ir N 84: 


benutze für die Klammerwerte die geometrische Reihe, die bei < >0 konvergiert, führe die 
Integration gliedweise durch, was erlaubt ist, und gehe mit e nach Null. Es ergibt sich: 


Ms, u, =(l-+ u): Alu)-s"—(1+0)-A(0)-s’— Bis, u,0). . . (32). 
Hierin ist: 


1 \ 1 
3 .® by 
AT u a ar eh (33), 
’ l 


=. 1 


—n En, 
AO Tr " »”—(1+0) 
vi 


(34), 


2) Für rationale Exponenten, wie sie bei der hier interessierenden Aı wendung der Integralgleichung vorliegen, 
ist grundsätzlich eine geschlossene elementare Darstellung des Integrals möglich. Die Formeln werden aber so unüber 
sichtlich, daß der hier angegebene, auch für irrationale Exponenten brauchbare Weg für die Anwendung vorzuziehen ist. 
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LE 
_ R 
\ u 5 ° 


u RI — 6 
’ 1) 


B(s, u,0)=(6 — u) 


wobei zur Gültigkeit von (32) die Zahlen « und o den Einschränkungen unterliegen: 
u ee TATEN el a 
und 
BRUT TE Te ME 

Die erste Bedingung stellt die Konvergenz des Integrals im Nullpunkt sicher, die zweite 
die Konvergenz der auftretenden Reihen, die sich bei der obigen Zusammenfassung bis auf 
s—1 erstreckt. 

Für « sind beide Bedingungen bei dem hier behandelten Problem erfüllt, in andern 
Fällen zu beachten; bezüglich o ist auf die Einschränkungen bei der späteren Wahl geeigneter 
Werte Rücksicht zu nehmen. 

Die Einschränkung «> — 1 ist, jedenfalls für hydrodynamische Anwendungen, durchaus 
plausibel und dürfte wohl die schwächste Anforderung sein, die man im allgemeinen stellen 
kann. Faßt man (26) als die maßgebende Gleichung für einen ebenen Tragflügel oder eine 
Schaufel auf, so sagt diese Bedingung nichts anderes, als daß die Normalgeschwindigkeit auf der 
Schaufel mit Ausnahme endlich vieler Punkte bestimmt und beschränkt sein soll, und andere Fälle 
sind kaum vorstellbar. Die genannte Bedingung ist auch die gleiche, die sich aus der Weissin- 
gerschen Bedingung der absoluten Integrabilität der mit sin 9 multiplizierten rechten Seite 
der Form (10) ergibt, womit man auf die in Gl. (13) rechts stehende Funktion kommt. 
Multipliziert man die ganze Integralgleichung mit sing, so kann man links die Integration 
über p durchführen; geht man dann wieder auf die Form (26) und wählt für die Wirbel- 

1 
belegung X (f) einen Potenzausdruck, so sieht man, daß es im wesentlichen auf \N(f)dt an- 


0 
kommt; also die Gesamtzirkulation um die Schaufel soll beschränkt und bestimmt sein, auch 
dies eine fast selbstverständliche Forderung. Die Verwendbarkeit des Ansatzes (30) scheint 
also ziemlich allgemein zu sein. Auf die Konvergenzgüte der verbleibenden Reihe B nach 
(35) und ihre Bedeutung wird weiter unten noch eingegangen; auf jeden Fall ist B im 
ganzen Intervall beschränkt und stetig; ferner ist B=0 für s—=V0. 


7. Wahl einer noch freien Konstanten; die dann noch verbleibende Integralgleichung. 
Für den Ansatz (30) wird nun gewählt: 

a, on 
AHWAW) treten: 9. 
Damit geht unter Berücksichtigung von (27) und (32) die zu lösende Integralgleichung (26) 
über in: 


ce 


) 
N, (N) 
\ s—t 
0 

Das Glied rechts mit dem störenden Exponenten « ist fortgefallen; die normale Potenz- 
reihe für B mit ganzzahligen Exponenten kann zu Unstetigkeiten keinen Anlaß geben, solange (36) 
beachtet wird. Die Konstante o ist mit den durch (36) und (37) gegebenen Einschränkungen 
noch beliebig wählbar. Dabei sollen aber neue Unstetigkeiten nicht auftreten; ferner wird 
man auf möglichste Vereinfachung der weiteren Rechnung bedacht sein. 

8. Wahl der anderen Konstanten zwecks Ermöglichung der Lösung. Beide Forderungen 
werden radikal erfüllt, wenn in (39) auch noch das Glied mit s” zum Verschwinden gebracht 
wird. A(u) ist mit « gegeben und bei ganzzahligem n > 4 stets beschränkt sowie von Null 
verschieden; mit o+1=0, also o=—1 würde (s+1)- A (0) unendlich, wie Gl. (34) zeigt; 
außerdem wäre hiermit (36) verletzt. Das Glied mit s” kann also nur verschwinden, wenn 
A (0) =0 möglich ist. 

Durch geeignete Wahl von o ist dies tatsächlich zu erreichen. Eine Untersuchung der 
Reihe in (34) zeigt, daß A (0), als Funktion von o betrachtet, Nullstellen besitzt, sogar wie 
sich sofort ergeben wird — unendlich viele. Es sei also 


| +0 A(0) o ı Ad, 


dt=a,. I+n Alu)” a + u)» Ala) 


DIT EG MER : 5 





a ee ee Imre. ME 
ein passender Wert; dann wird aus (39): 
1 
we 4, i 
r r di=7 Fp)- Aka) DEAN. en gt (41) 


0) 
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mit der Bedingung: 
N,1)=0. 


9. Geschlossene Darstellung der Funktion A (0); ihre Nullstellen. 
Werte für o, wird dadurch erleichtert, daß die Funktion 
geschlossenen Ausdruck darstellbar ist. Für cotg « gilt nämlich die „Partialbruchverlegung“ 


(s. z. B. „Knopp, Unendliche Reihen“, S. 212): 
_. 
OMPAEE, ir, 
ax. otence=1-+2ır- 
| 
oder 
4 
2 eotgam. ı +9 . | 
x we” u — 
v—| 


Vergleicht man dies mit der Definition von A(6) nach Gl. (34), so ergibt sich: 


cotgr(l-+ 0) 
: we 
Aus (43) folgen sofort die Nullstellen von 4 (0) für: 
IR 8 
a 


A(o)=n 


+(o-+1) - 


unter Beachtung von (36) kommt aber nur in Frage: 
re, 
ve 


0 )» >) > 
also 
| 0 
9 5 tk; k=0,1,23,3----- 


10. Vollständige Lösung für n >4. Mit 


69 a [42] 


ei (dl +)Alu) 2.178 (0, 


u 


und 
1 
BE (u—v)(o,— vr) 
wird jetzt die zu lösende Gl. (41): 


1 vn 
"N,(d 3 . 
\ . dt=a;- P. Yy 8’ 
0 v 0 


(42), 
und dies ist nach der für n<4 benutzten bekannten Methode ohne Schwierigkeiten lösbar. 


Die Auffindung passender 
A(o) bekannt und durch einen 


(43). 


(4). 


(45) 


(46). 


Zur Lösung nach der bekannten Methode, die noch kurz skizziert werden sell, seien 


wieder eingeführt: 


cos yr 


siny- N, ( > j= Y(y) 


womit man erhält: 
rt 


Rn 
l Y(y a Bi 
5) \ a dy = = ; Y Yı 
=T) cos 7 cos y BHO Ka 
0 ' 0 


mit 
Y(a) =0 
Es wird wieder angesetzt: 


| - 
Yy)= 5 b,+ B h,-cosiw. 


=| 


1 1 
s——(1 cosgy); = ll cos y); | 


| 
| 


dann ergibt sich, daß zur Bestimmung der 5; die Entwicklung 


(+) an 
1 2 BER a 3 
3 \ b,-sinig=;! > ) 
I _ er 
4i—1 v—() 


',„. sin?’ ol2- sing 


(47), 


(48). 


(49); 


(50) 
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durchzuführen ist. Irgendwelche Unstetigkeitsschwierigkeiten gibt es nicht mehr; die Reihe 
rechts konvergiert gleichmäßig für alle 9; man darf gliedweise integrieren; die entstehenden 
Reihen für die b, sind ebenfalls konvergent. 

Man erhält: 


n 7 
2a 0 ua Ep r 
b=—- \ y+\sin y-sin’’ol2-sinäp-de. . .». . 2.2... 
e ” mn) h 
‚=Uv 1) 


Herr Frieke hat die Integrale ausgewertet; es wird: 


’ 114-1 EEE BE j 
bj: z | j Ib ty ya) + \ PA, ”)| re A 
worin: ce 
x (2r+1)-2r (2» —1)(2r — 2) 2Qi+1)-24 
PAN= 5, IWW -A-V & -W_-@i-j @i+1®-- @i 1) j 
(53). 
ER. f (2r+1)-2» (24+3)(2/4+2) 
| ıd +Y-AIH (2i+3”’— @i +1) 
Hiermit wird bei Rücksubstitution: 
\,D= 4 u.+ v 77 -Cos Ay (4): 
* oe‘ 2] t R > o j k x G / . . . . . 1 * 
4 1 
zur Erfüllung von N, (1)=0 führt wieder 
1 " nn 
2 b, — \ b; 1 )+ . j - 2 r 2 . IP»), 
4 1 
und hiermit ergibt sich: 
2b; -[eosiy —( 1%] 
e Ad, el .n 
== (fi 9 bee MH k 
ABER ZEN ayt-e BEL RIET.. 


< 


Führt man wieder alle Rücksubstitutionen aus, so erhält man schließlich mit = 1 : 


= —=2.tg : 1 En(l +00) - +2 yi zu’[0s = -+ 1) . \ by [eos y — (- 1)A]| (57). 


=] 
Damit ist das Problem auch für » >4 in der beabsichtigten Weise gelöst. 


11. Einige wichtige Einzelheiten; Zahlenrechnungen. Folgende Einzelheiten seien noch 
erwähnt: 

a) Trotz der Einschränkung (37) kann o, nach (44) noch unendlich viele Werte erhalten; 
das Eigentümliche ist, daß hierdurch der Verlauf von m (£) nach (57) nicht geändert wird. 
Durch Vergleichsrechnungen mit ,=—1/2 und o,=+1/2 wurde dies auch ziffernmäßig 
verifiziert. Das Problem ist eben, wie aus der Behandlung mit Hilfe der Potentialtheorie be- 
kannt, mit der Bedingung m (l)=0 eindeutig. 


b) Nach Abschnitt 6 und 8 hat sich die einfache Integralformel ergeben: 


1 


oo BE - N 13 
\.u=> ER SFETET" 
0 v‚—( 


Hierfür sind auch geschlossene Ausdrücke möglich. Z. B. ergibt sich für den kleinsten 


Wert 0,= — der sich auch für die Auswertung als’ zweckmäßig erwiesen hat: 


S%) 


— „ug Ss. 
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Man erhält dieses Ergebnis entweder durch einfache Umformung der Reihe, wodurch 
sie in die bekannte logarithmische übergeht, oder durch die Substitution z= y t in dem be- 
stimmten Integral und konsequente Bildung des Grenzwertes, wie in Abschnitt 8 angedeutet. 

Für die Anwendung im vorliegenden Problem sind diese Ausdrücke nicht zu empfehlen, 
da sie für s>1 unendlich werden, ebenso wie die in (32) nach Abzug der Reihe für o, noch 
verbleibende mit der zwischen — 1/2 und — 1 liegenden, sonst aber beliebigen Zahl u. 

c) Die Reihe für B (s, u, 0) konvergiert besonders für Werte von s nahe an I nicht gut. 
Für die Auswertung ist dies aber nicht wesentlich, da das entsprechende Glied gegenüber 
den anderen nur einen Bruchteil ausmacht. Man sieht dies bereits an dem Faktor o, — u; 
man wird natürlich o, möglichst in der Nähe von « und ebenfalls negativ wählen, womit sich 


in Übereinstimmung mit den bisherigen Auswertungserfahrungen der kleinste Wert o, = > 
als empfehlenswert ergibt. 
Wählt man für ein Beispiel den mittleren Wert «= j entsprechend » =, so wird 


an der für die Konvergenz ungünstigsten Stelle s=!: 


\ 23 dti=n 3 R: 
wobeı sicherlich E 
\ v s’ I a? 
R< (u 0,) „e ww ı 
Y 1 N 1 


der nach einer langsam konvergierenden Reihe zu berechnende Betrag macht also nur etwa 
ein Sechstel des geschlossen vorliegenden aus. Dies gilt für die ungünstigste Stelle s=1; 
für kleinere s liegen die Verhältnisse noch wesentlich günstiger; für s>0 geht auch die 
Reihe nach Null, hier bleibt allein der nach Unendlich gehende Anteil übrig, der geschlossen 
erfaßt ist. 

Bei häufigeren oder umfangreicheren Anwendungen wird man das Bestreben haben, die 
Abhängigkeiten von den drei Variabeln s, «, o, nach Möglichkeit zu trennen, so daß man im 
wesentlichen ein für allemal zu berechnende Funktionen von s erhält, die nun von Fall zu 
Fall mit leicht berechenbaren Koeffizienten in « und o, zu multiplizieren sind, wobei man o, 
natürlich fest wählen wird. Für solche Zwecke sei beispielsweise das folgende Ergebnis mit- 
geteilt; es wurde nach dem sogenannten Kummerschen Verfahren (s. hierzu „Knopp“, Ab- 
schnitt VIII; $ 33) durch sukzessives Abziehen und Addieren der hier naheliegenden Funktionen 


9 
\ - erhalten. 


n 


u 
l 

Es sei: 
0 1 E o ; 

y r E Ss N Ss 
CG=u+% . F,(s) \ ; („@o)= \ 3; 
(‚= wW"+uo, +60; Fr v un FEW 9) 

= ui Ei nn ’ ’— 
C,=e+0o+u®-+ 0} 
> = Fr 

2 
n=y ” ‚ \ . s 
7 R 2 Or, a) 
C, \ — A F,(s)= N 77 mm Per e—e) 
Pe en 4 u v » l 
n u ’ l 
n,v=0,1,2,3..... q=2,8,4,5..... .=@a188:. 
Dann wird: s 
no p=zn -23,q0=n 
\ ” er \ ie; 
= ARE; e -f uw" oc a Er a). 
(v u) (v 6,) D 4 +6, 1 fin 7 In (& 
‚—| p 0d,g=2 


Das gesteckte Ziel ist offenbar erreicht: die C, sind leicht berechenbare Koeffizienten 
in « und o,; die F,„ sind Funktionen von s allein, können also ein für allemal berechnet 
werden, so daß sich eine größere Rechenarbeit lohnt; f„_, ist allein von o, abhängig und 
kann ebenfalls ein für allemal berechnet werden. Das einzige noch von allen drei Variabeln 
abhängige Glied ist «”-,„; dies kann durch passende Wahl von n beliebig klein gehalten 
werden, wobei gleichzeitig die Konvergenz der verbleibenden Reihe beliebig verbessert wird. 

Ob es sich nach den früheren Bemerkungen über die Bedeutung des Reihengliedes 
lohnen wird, die Funktionen F, welche die Hauptarbeit machen würden, zu tabulieren, muß 
von Fall zu Fall beurteilt werden. 
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Die Form (a) habe ich bis » =4 nachgerechnet und dann das Bildungsgesetz erschlossen, 
ohne noch rückwärts den strengen Beweis durchzuführen. 

d) Auch die Schlußbeziehung (57) zeigt (in Übereinstimmung mit (37)), daß o,: 1 
unmöglich ist; denn hiermit würde das zweite Glied in der Klammer je nach der Wahl von o, 
in verschiedener Weise unendlich für {>0, was der Eindeutigkeit widerspricht. 

e) Die durchgeführten Zahlenrechnungen ergaben volle Übereinstimmung mit den ein- 
gangs erwähnten, auf anderen Wegen ermittelten Ergebnissen. Für die Konvergenz der auf- 
tretenden Reihen erwies sich o,= — 1/2 als besonders vorteilhaft. Für die Reihe in (57) ge- 
nügte im allgemeinen die Berechnung von nur 2 bis 3 Gliedern, wenn die entsprechenden b; 
bestimmt sind; hierfür waren mehr Reihenglieder erforderlich; die Rechenarbeit bleibt aber 
in durchaus erträglichen Grenzen, ohne Anwendung der Umformung (a). 

f) Für die Anwendung wichtig ist noch der Wert /' der Gesamtzirkulation um die 
Schaufeln, aus dem man nach Überlagerung einer Quell- oder Senkenströmung vom oder zum 


Nullpunkt sofort die Beziehung für "das Drehmoment und für die Leistung erhält. Es ist: 
1 
!=n-\m(£)-d®. (DS) 
u 
Die Integrationen lassen sich ausführen (für »=1,2,4 in geschlossener Form); man 
erhält: 
3 
= gn Fir ami, 
5. 
fezwu Pf une, 
BE - 5 (9). 
"»—-uf Meun=i, 
T 
n 2 | 2 a 2 u | 
"—_ tor . 2, (A, * 63 . 1)* - für >4 
R en 's “nn ‘ n(ito) n "ı \% a da \ i 2% 
4 L 





Die letzte Beziehung liefert für »=6,8,10 die auf 2 Stellen abgerundeten Werte 
P)o1?—=0,10,0,75,0,79. Für n>x gehen die beiden letzten Glieder in der eckigen Klammer 
nach Null, der Faktor von &l? vor der Klammer nach 1, so daß 


!>o! fürn>x, 


wie es sein muß. 


12. Näherungen für schwach gekrümmte Schaufeln; ein Beispiel. Es sei noch gezeigt, wie 
man mit Hilfe der dargelegten Methoden angenäherte Lösungen für allgemeinere Fälle er- 
halten kann. i 

Die » unter gleichen Zentriwinkeln liegenden Schaufeln seien jetzt schwach gekrümmt; 
y=y(x#) sei die Gleichung der Schaufelkurve in der w-Ebene. 

Ich setze nun voraus, daß y’ und y/x so klein sind, daß sie gegenüber 1 als von erster 
Ordnung kleine Größen angesehen und ihre Quadrate gegenüber 1 als Größen zweiter Ordnung 
gestrichen werden können. 

Auf diesen Schaufeln sei eine Wirbelbelegung aufgebracht. Dann ist sofort ersichtlich, 
daß die Normalgeschwindigkeiten, die dadurch auf den Schaufeln entstehen, nur um Größen 
zweiter Ordnung von denjenigen abweichen, die durch die gleiche Wirbelbelegung als Funktion 
von x bei y=0, also auf dem Radius angeordnet, auf dem Radius erzeugt werden. Denn in 
dem Ausdruck für die Normalgeschwindigkeit kommen die Abstände zweier Punkte vor, die 
um Größen zweiter Ordnung von ihren Projektionen auf die #-Achse abweichen, und der 
Cosinus kleiner Winkel von der Größenordnung y’ und y/x, der von 1 um eine Größe zweiter 
Ordnung abweicht. 

Das Analoge gilt für die Normalkomponenten der von der Rotation herrührenden Ge- 
schwindigkeiten der Schaufelelemente. Auch hierin kommen lediglich Quadrate von y/x 
und 4’ gegenüber 1 vor. 

Wenn also die Schaufeln in bezug auf y’ und y/x nur um Größen erster Ordnung vom 
Radius abweichen, so können sie bezüglich der Rotation in erster Annäherung wie radiale 
(rerade behandelt werden. 

Anders liegen dagegen die Verhältnisse bezüglich einer Quellströmung, die noch hinzu- 
kommen muß, damit das Schema einer Pumpe vollständig wird. Im Ursprung sei eine Quelle 
von der Intensität Q@ angeordnet, derart, daß im Abstande r eine nach außen gerichtete Radial- 
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geschwindigkeit Q@/r entsteht. Sind die Schaufeln radiale Gerade, so werden die Rand- 
bedingungen durch die Quellströmung nicht beeinflußt. Bei den nach der Voraussetzung 
dieses Abschnittes schwach gekrümmten Schaufeln bildet aber in einem Punkte x, y die Radial- 
geschwindigkeit Q/x (bei r=.x) mit der «-Achse den Winkel y/x, die Schaufeltangente da- 
gegen den Winkel y’. Die Radialgeschwindigkeit der Quellströmung bildet also mit der 


7. ’ 2 ® r [#] ’ 
Schaufeltangente den Winkel y’— y/.x, so daß eine Normalkomponente - x W—-ylx) ent 
steht. Diese ist von erster Ordnung und muß daher berücksichtigt werden. Die zum Aus- 
i ’ SE 
gleich erforderliche Wirbelbelegung mg mit den Bedingungen », — x (y — ylx) und my ()=0 


kann aber wieder auf dem Radius angeordnet werden, da hierdurch nur Fehler von zweiter 
Ordnung entstehen. 
Hiermit ergibt sich folgendes Problem: 
Die Schaufelform in der w-Ebene sei gegeben durch 
y=>a,: x. en er u IR ur er 
Die Reihe rechts habe je nach Bedarf endliche oder unendliche Gliederzahl. Es sei: 
yle= Na, 2"'<1 und y=}r-a,:x 3 ASIEN 4° 
so daß die Quadrate hiervon gegen 1 gestrichen werden können. 

Im Ursprung sei ferner eine Quelle von der Intensität @ angeordnet. Dann wird nach 
den vorstehenden Darlegungen unter Berücksichtigung von Gl. (1) und den dazugehörigen Er- 
läuterungen das hier interessierende Strömungsproblem gelöst durch eine Wirbelverteilung 
m (£), für welche folgende Integralgleichung gilt: 


l 
"m(HdE 1) an 
na" \ A nznict—: \ (r —- 1)-a,-x"”%, 
« A " JE —— 
0 
oder: 
1 
m(d) ee ) 4 l : 
\_ _;n di= n rt n i R, (r I) ar: yu-r+i re 
0 
mit m(l)=0. 


Da die Gleichung linear in m (£) ist, handelt es sich bezüglich des Summenausdrucks um die 
Lösung von: 


\ i (£) konst. (63), 


n d Fass ar 


wo p jetzt im Gegensatz zu früher ein von n unabhängiger Exponent ist. 
Zunächst sei geprüft, wann (63) nach der Methode von Abschnitt 4 ohne das Auftreten 
von Unstetigkeiten gelöst werden kann. Die Wiederholung der zu den entscheidenden 


6) 9 _9 
Gl. (13) bis (17) führenden Schritte zeigt, daß hier jetzt K, -— an Stelle von K--®- 
zu setzen ist. 

Gl. (63) ist also praktisch nach der einfacheren Methode von Abschnitt 4 zu lösen, wenn: 


2 


a, 98 u _n ‚ 
K=, =! d. h. wenn a ze oe ee (64). 


Wendet man dies auf die Glieder des von der Quellströmung herrührenden Summen- 
ausdrucks von (62) an, so ist p=n— r+1 zu setzen. Die Bedingung (64) gibt dann: 


a 2 NR 
Bezüglich der Anwendung der zweiten Methode (für an >4 bei radialen Schaufeln) er- 
—2 
gibt sich folgendes: In Gl. (26) ist «= - 5 „m ersetzen durch „=—£. Für u, gelten 


dann die Einschränkungen (36) und (37); d. h. die zweite Methode ist anwendbar, wenn 


£+ ganze positive Zahl einschließlich Null, 


und 
Be En a u 
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Ferner darf (1+ u,)- A(u,)=reotgr (1 + u,) nicht Null werden; für die Anwendbarkeit 


der zweiten Methode gilt also für die weitere Einschränkung: 


= -, n n = 


p 1 1 . 
"pe Ei Per, 5, 7 2 70.72 


n - _ = - o pr ö 


Für die Exponenten r in (60) und (62) ergibt dies die Einschränkungen: 


, 1 
1- ++ ganze a E95 ee rer 
was bei ganzzahligem r und » > 1 immer erfüllt ist; ferner 
ee ME re 
und 
r+#1+k- mit k=1,3,5-.:-- a ET en 


Dabei sind die Fälle nach (70) mit der einfacheren Methode zu erledigen. 
Zusammenfassend heißt dies: Für ganzzahlige r > 1 ist die Integralgleichung (62) nach 


. : . = S n . 
den entwickelten Methoden immer zu lösen; für r=1+7 kommt die Methode nach Ab- 


. n . R . . . 5 
schnitt 4, für1+5 >r>1 die zweite Lösungsart in Frage. Daß r=1 ausscheidet, bedeutet 


praktisch keine Einschränkung; dies heißt lediglich, daß die «-Achse im Ursprung tangential 
an die zu untersuchende Schaufelkurve zu legen ist. 
Das Vorstehende ergibt ein hübsches Beispiel für r —3. Gl. (62) geht hiermit über in: 


I 
m(d 1 1 
s —— D) . 
\zr am dE: n (w+2a,- 0) an-r 
v 


wobei die Schaufelform gegeben ist durch 
ya“, 


was bei a, >0 vorgekrümmte, bei a,<- 0 zurückgekrümmte Schaufeln im Sinne der Rotation 
bedeutet. Die Wirbelverteilung m (£) und die Gesamtzirkulation /' fällt dann so aus, als ob 
das Schaufelrad mit der Winkelgeschwindigkeit 


9,—=o+2a,:® 





rotiert. Bei r=3 (kubische Parabel als Schaufelform) ist also im Rahmen der getroffenen 
Annäherung die Randbedingung für die Radialströmung die gleiche wie diejenige für die 
Rotation. 

Die Schaufelzahl sei z. BB n=4. Dann ist nach (59): 


2 
Tan: ?= - (o+2a,:0)-". 
Die spezifische Masse der Flüssigkeit sei o=y/g; dann ist das Drehmoment pro Einheit der 
Höhe senkrecht zur Strömungsebene: 


M=o:27-0: 1'—40:0-(w+2a,0)- 1; 


die Leistung: 
L=M:-»=409:-o(»+2a, 0)- 1? 


und die sogenannte Bruttoförderhöhe des Pumpenbaues: 


L 2 © R g 
M= 9 gerad 


Bemerkenswert ist noch, daß für p=n keine Lösung von (63) existiert, wenigstens 
nicht nach einer der beiden hier behandelten Methoden. Physikalisch ist dies plausibel. 
Denn p=n entspricht, wie leicht zu ersehen, derjenigen Randbedingung für die Normal- 
geschwindigkeiten, die bei Vorhandensein eines konzentrierten Wirbels im Ursprung anzusetzen 
wäre. Ein solche ist aber schwer vorstellbar, wenn sich die Schaufeln, wie hier vorausge- 
setzt, bis in den Ursprung «= erstrecken. 
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13. Lösung bei nicht vom Nullpunkt ausgehenden Schaufeln. Zur Vervollständigung sei 
noch kurz auf die Verhältnisse eingegangen, die sich ergeben, wenn die » geraden Schaufeln 
sich nieht bis in den Nullpunkt erstrecken. 

Es sei also: 


wesen Er 5 Ta Re 
Die Integralgleichung lautet jetzt: 
Is 
m(d) TR. 3 
\- a di= n "gen 1 SE N ra 5 a ee (72). 
T, 


Die untere Grenze des Integrals ist jetzt !, anstatt bisher Null. 
Mit den Substitutionen 
"—- (mM—M).t+Hl, "le BETON 
1 (73) 
t=—>(1—-cosy); s—5(1-e0sg) | 


erhält man jedoch wieder die aus Abschnitt 4 geläufige Form der linken Seite mit der unteren 
Grenze Null; nämlich: 





\ rw) l nn? P1/ 
’ +(C N 2 I» ı 
\ 6089 coS y U a, -(a— sın? o/2) 
u 
mit 
[277 
a, 
I, \nj“ 
mo-lı (4) 
Pe a A GEEHRTE 
()' 
I» 0 
[24 (") . 
l; 
2 
u —] 
n 
Dies führt mit 
e N B, RR 
Yıy) ut \ b,.cosry. DE ee ee a 
„1 
auf die Entwicklung: 
” . a B Sr u; 
\ b,-sinvg = sing (a4 a 2 Er ET 
} l ; 


Die Funktion rechts ist jetzt bei a > 0 für alle « im ganzen Intervall stetig, so daß die 
Fourierentwicklung stets ohne Schwierigkeiten möglich ist. 

Die Methode von Abschnitt 4 ist also bei allen Schaufelzahlen ohne weiteres durch- 
führbar; und da die Konvergenz der auftretenden Reihen gut ist, spielt es für die Anwendung 
keine wesentliche Rolle, ob sich die jetzt auftretenden Integrale zur Bestimmung der b, 
geschlossen auswerten lassen oder nicht. Aus Zeitgründen mußte die nähere Bearbeitung 
dieser Integrale unterbleiben. 

Da jetzt zwei freie Schaufelenden vorliegen, steht es zunächst offen, ob man m(l,) oder 
m (l,) zu Null macht. Man wird jeweils die Austrittskante wählen, so daß bei einer Quelle 
im Ursprung m (l,)=0, bei einer Senke m (l,)—0 zur Bestimmung der Konstanten 1/2 b, in 
(75) zur Verfügung steht. Am jeweiligen anderen Schaufelende hat man dann m = x, also 
„Eintritt mit Stoß“. Bekanntlich kann man dies durch Überlagerung einer Zirkulation um 
den Nullpunkt von geeigneter Intensität auch noch vermeiden, so daß man dann den speziellen 
Betriebszustand für „stoßfreien*“ Eintritt und Austritt erhält. Auch die hierfür zu be- 
friedigenden Randbedingungen führen bezüglich der betreffenden Fourierentwicklung _zu 
keinen Unstetigkeiten. 

Die zweite Methode (n >4 bei radialen bis zum Nullpunkt gehenden Schaufeln) ist 
jetzt nicht erforderlich. 
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Auch jetzt sind die Annäherungen erster Ordnung für schwach gekrümmte Schaufeln 
möglich; die hierbei auftretenden Randbedingungen führen für die Fourierentwieklungen auf 
stetige Funktionen. 

14. Zur Umkehrung des Problems: Bestimmung der Schaufelform bei gegebener Abänderung 
der Wirbelverteilung. Zum Abschluß sei noch ein Ausblick auf die Möglichkeit einer Um- 
kehrung des Problems gegeben. Bisher galt die Schaufelform als vorgeschrieben; gesucht 
wurde die zugehörige Wirbelverteilung. Diese steht, wie schon erwähnt, mit den Druck- 
und Geschwindigkeitsunterschieden auf den beiden Seiten einer rotierenden Schaufel in 
unmittelbarem Zusammenhang, also mit Größen, durch deren Verlauf offensichtlich das Auf- 
treten von Ablösungen, Wirbeln und sonstigen Verlust- und Störungsursachen beeinflußt 
wird. Für die Anwendung auf die Bedürfnisse der Praxis liegt daher die Umkehrung der 
Fragestellung nahe: Gegeben sei eine Wirbelverteilung; welches ist die zugehörige Schaufel- 
form? Dies hat nur einen Sinn, wenn man gekrümmte Schaufeln voraussetzt: im Rahmen 
dieser Abhandlung, in welcher nur schwache Krümmungen annäherungsweise berücksichtigt 
wurden, kann es sich bei der angedeuteten Umkehrung der Fragestellung also nur um be- 
absichtigte Korrekturen der für die geraden Schaufeln erhaltenen Wirbelverteilung durch 
schwache Schaufelkrümmungen handeln. 

a a NEN ee: 
sei die zunächst noch nicht festgelegte Gleichung der Schaufelkurve; dabei sei y.e<1 und 
y'<1, so daß die Vernachlässigungen nach Abschnitt 12 zulässig sind. 

Die Integralgleichung (62) lautet hiermit: 


l 
\ m (£) 1: [02) l 1 ® & 
Ge — . > . (‘ ul) Ss). 
ar m n g"-: zn „ „ (9) 
u 
Es sei: 
FED... 3.2. 2 2 rer 
wobei 
I 
m,(S) 2 [7) 1 
\= md ae Se ar EEE 


d.h. m, (£) sei die bei geraden Schaufeln nach einer der beiden angegebenen Methoden zu 
bestimmende. also als bekannt anzusehende Wirbelverteilung. Hiermit verbleibt: 


m,(£) a f = 
Fa zn ds WU Yy a EHEN ch ER 59a. SE 


m,(£) werde nun als eine nach gewissen Gesichtspunkten zu wählende Korrektur von 
m,(£) angesehen; man will z. B. den Abfall von m (£) am Schaufelende weniger schroff ge- 
stalten oder die Geschwindigkeits- und Druckunterschiede auf eine gewisse Erstreckung der 
Schaufel etwas ausgleichen. Jetzt ist m,(£) eine bekannte Funktion von £, die man sicher- 
lich so wählen oder annähern kann, daß das Integral links in GI. ($1) als bekannte Funktion 
von x darstellbar wird. 

Gl. (81) ist jetzt eine Differentialgleichung erster Ordnung für y(#), d.h. für diejenige 
Schaufelkurve, mit welcher die gewünschte Korrektur der Wirbelverteilung erzielt wird. 

Mit Hilfe der Substitutionen (4) und (25) kommt man auf eine einfachere Form des 
Integrals; in Abschnitt 6 ist gezeigt, wie dieses für Potenzausdrücke im Zähler zu be- 
handeln ist, wobei sich Polynome und Potenzreihen in der an die Stelle von x tretenden 
Variablen ergeben. Damit eröffnet sich die Möglichkeit, y(#) für die gesuchte Schaufelform in 
entsprechenden Polynomen und Potenzreihen anzusetzen, die Koeffizienten dureh Vergleichung 
zu bestimmen und damit das angedeutete umgekehrte Problem in dem angegebenen Rahmen zu 
lösen, d. h. eine Schaufelform aus vorgegebenen Eigenschaften der Strömung zu bestimmen. 
Der unmittelbare Nutzen solcher Lösungen für die Praxis, der die Fragestellung an sich sehr 
nahe liegen muß, ist dadurch eingeschränkt, daß hier nur schwach gekrümmte, in der Nähe 
eines Radius verlaufende Schaufeln behandelt werden, die für die meisten normalen Zwecke 
reichlich extrem sind. Andererseits ist aber der Einfluß der Schaufelzahl im vorstehenden 
ohne Vernachlässigung erfaßt, und bei hohen Schaufelzahlen ergeben auch schwach gekrümmte 
Formen durchaus brauchbar erscheinende „Kanalformen*“ für die Relativströmung. Im 
übrigen sind die vorstehenden Überlegungen nicht so sehr darauf gerichtet, Ziffern und Formen 
für den unmittelbaren Gebrauch der Praxis zu erhalten; der Hauptzweck ist vielmehr die 
Erfassung grundsätzlicher Zusammenhänge an dem wohl einfachsten sich darbietenden 
Schema. 252 
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Zur geometrischen Deutung und Transformation der Grund- 
gleichung der ebenen kompressiblen Potentialströmung. 
Von M. Pinl (Messerschmitt- A.G.) in Augsburg. 


Das Geschwindigkeitspotential der ebenen adiabatischen kompressiblen Po- 
tentialströmung genügt einer quasilinearen partiellen Differentialgleichung 
zweiter Ordnung, die sich nach G. Braun als Eulersche Gleichung eines 
Variationsproblems auffassen läßt. Werden die beiden darin auftretenden physi- 
kalischen Konstanten A und u gleich I, so erhält man die Differentialgleichung 
der gewöhnlichen Minimalflächen. Entsprechend lassen sich die Integralflächen 
der Differentialgleichung für allgemeine, d.h. physikalisch vorgegebene A, u, als 
verallgemeinerte sogenannte Relativminimalflächen interpretieren, die Diffe- 
rentialgleichung selbst als das Verschwinden einer mittleren Relativkrüm- 
mung. Allen kanonischen Formen der Differentialgleichung entsprechen 
dann relativgeometrisch ausgezeichnete Parameterkurven auf den Integral- 
flächen. 


$ 1. Einleitung. 


Sind x, y unabhängige kartesische Koordinaten, so genügt das Geschwindigkeitspotential 
2 (x,y) der ebenen kompressiblen Potentialströmung der quasilinearen') partiellen Differential- 
gleichung zweiter Ordnung: 


Kt +mM)—-pilr—2pas+ß+tul Pe +) -Alt=0 .....M. 
Darin bedeuten: 
02 02 0° z 0’ z 0% z ’ 
vr q y' ’ es 7 L 8: dzdy’ du 7% ARE. 


Die physikalischen Konstanten 4, « sind durch 
„—1 I— 1 


i=,+—-m, p=—a a ee IE 2 
gegeben. x bedeutet den adiabatischen Exponenten, «a, und w, die Werte der Schallgeschwin- 
digkeit a=a(mw) bzw. des Geschwindigkeitsbetrags w in einem zur Normierung benutzten 
Punkt’). 

Für A=u=1 entsteht aus (1) die Differentialgleichung der bewegungsinvarianten 
Minimalflächen des dreidimensionalen euklidischen Raumes 

A+fr-Eras ti +PNt=0 . . . . .» 2 2... 

deren Lösungen man bekanntlich in integralloser Darstellung vollkommen beherrscht°). Dieser 

Umstand legt es nahe, die Lösungen von (1) als verallgemeinerte Minimalflächen aufzufassen, 

wie das gelegentlich von A. Busemann*) und H. Bateman’°) angedeutet worden ist. Um 


an die Differentialgleichung (1) invariante Methoden heranzubringen, behandeln wir sie nach 
dem Muster von G. Braun‘) als Eulersche Gleichung eines Variationsproblems. 


$ 2. Relativminimalflächen. 


+1 „ 


74 


. ' l ie ” . 
Es sei #1, Fran Ku! Am 1 77 und ww „= | „ga, wo a die kritische Schall- 
geschwindigkeit bezeichnet. Dann ist (1) die zu 
öl\atrdedy=d|\|[A +. +eWdzdy=0. .: ...... 8 


') Vgl. (zur Terminologie) Courant-Hilbert: Methoden der mathematischen Physik, Bd. II, Kap. 1, 8.2: 
Berlin 1937. 

2) Vgl. G. Braun: Die ebene kompressible Potentialströmung als Variations- und Eigenwertproblem, Annalen 
der Physik (5) Bd. 15/6 (1932), S. 645. 

») Vgl. z. B. W. Blaschke: Differentialgeometrie Bd. I, $ 107, S. 238, Berlin 1930. 


4) Vgl. A. Busemann: Aerodynamischer Auftrieb bei Überschallgeschwindigkeit, Atti dei Convegni, Fon- 
dazione Alessandro Volta Bd. 5 (1936), Quinta Seduta, 328 bis 360, insbesondere 348 bis 352, 355. 
5) Vgl. H. Bateman: Proe. Royal Soe. London (A) Bd. 125 (1928), S. 598 ff. 


6) Vgl. (9) 93 und (£). 
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; u a a 1 . 
gehörige Eulersche Gleichung’). Für = «=1, y=- geht (5) in 


s\\yl+p+fdaedy er a a 1 ET Fr 


über, d. h. in das bekannte Variationsproblem der Minimaltlächen. 

Die physikalische Bedeutung der Gl. (5) ist: das Integral über die 2y-te Potenz der 
Schallgeschwindigkeit ist stationär. Die geometrische Bedeutung der Gl. (5) ist: die „Relativ- 
oberfläche* \\a?’dxdy ist stationär. Im Sinne der relativen Flächentheorie*) wollen wir 
entsprechend die Integralflächen der Grundgleichung (1) als „Relativminimalflächen“ deuten. 
Führen wir allgemeine Flächenparameter ein und stellen die Lösungen von (1) durch den 
Vektor r= tx! (u, v), x? (u, v), #°’ (u, v)} dar, dessen partielle Ableitungen wir 

ÖL or 


eu: seh .. uw 
X, Lu du’ Ya LK dv zZ ER 


- En 


schreiben, so erhält der Normalvektor p =r,.% x, die Komponenten 


x 8 “ \ „ı 


z Kur Lu Lulu Kur Lu ) 
pP: nr >=, „1, a re 
In, In E Fi PA ' Can, In 
vr ı BIER a 


Für die Transformation des Doppelintegrals \\a’' dxdy erhalten wir 


> —=2(2 (u,v),x°(u,v),, p=— Pı TE a ° 
Ps Ps 
\ \ P(p)dude \ \» i+u Pi pP > a RE SEEE, (10). 
“. .. 3 








Der Integrand P (p)=|i+ u (p? + q?)]' genügt der Homogenitätsbedingung 
oo» \ 
Pir=p Pr tPp: Pra+tPpı Pas = P: (®», “35 m 7 7 m 3 CP | 1 
I 


Aus p gewinnen wir den Vektor n—p®"', dessen Komponenten 


Pi nn & 
ng FE en I 
in p; homogen von nullter Dimension sind. Nach (11) gilt 
nernte iu a. re 
Daraus folgt durch Differentiation 
RE > p Dy (Pa PD p D,) D, D ii D PD, ‘ 
v-n=—ın,N=— Bl , db: d: =6, Et; :: 


Zufolge (14) liegen die von den Vektoren üi„ und it, aufgespannten Tangentialebenen der 
Fläche 

FEREREIETN > no. 3 Ba 
in entsprechenden Punkten parallel zu den Tangentialebenen der Fläche x. Ferner gilt: 


\\Ppdudr=\\p®,dudr=|\\| Dix )dude=\\(e,1u,%)dude=0. . (16). 


Das bedeutet: der Vektor e=n=®, mit den Komponenten 
Du, = 2yupltu +", Pa=—2rnaitu +", Du. . (m 
ist der dem Variationsproblem (5) entsprechende Relativnormalenvektor, 
das Doppelintegral \\(e,x„,%,)dudr die Relativoberfläche und die Fläche 
e(u,v) die Eichfläche der Extremalflächen x). 
Jetzt betrachten wir das Variationsproblem 50=0. Seine Eulerschen Gleichungen 
lauten: 
RS HEN... AT NR R N PO 
Zufolge der Abbildung durch parallele Tangentialebenen gilt: 
eat, Bit. rn 
7) Ygl. (6) S. 652, 653. 
s) Vgl. A. Duschek: Sitzungsberiehte Akad. Wien IIa, Bd. 135 (1926), S. 1 bis 8; Bd. 136 (1927), 8. 265 bis 270, 
®%) Vgl. 9, 8.2, 4. 
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Das bedeutet für (18) 
KK -aAXD)=- a + UuXD)=0O .. . 0... (0). 


Wegen 8 „X e,=+0 ist somit das Eulersche System (18) mit dem Verschwinden der Invari- 
anten a-+d äquivalent. Andererseits gilt wegen (19) notwendig 


Re en Te 
Daher ist (18) mit dem Verschwinden der mittleren Relativkrümmung '°) 


a+d Br’ eu: X.) + (f, Lu» e,) 


RE 
u "Bea -u- 2 (f, tu, &v) 


äquivalent. Dabei bedeutet f einen willkürlichen Vektor, für den nur (f,x„,2.) #0 gilt. Die 
Darstellung (22) ist von f unabhängig, wie man mit Hilfe von (19) erkennt. 
« Für die Parameterdarstellung 

"== u, rmymu,mzln)) » .» 222 (BB) 


bedeutet das Verschwinden von H*, wenn wir die Komponenten von f mit k,, k,, k, be- 


zeichnen: 
k, 5) (D, x; 0 k, ) 1, (D,.)u 


k,, (Du), 11+16,,0, (8, )yi=0. 
k,, (Dos)x: q k,,P,(Py3)u 
Wählen wir f=$0,0,1}, so folgt: 
0,1,0 
(PD) tlP,),=0; bi) = 0,01 =1 
1,9, 4 
Setzen wir A +u(pP +) "= (p, Q), so gilt nach (17): 
(Dada=— 2ruler +Poor tg; (Brdn=  2rulwt+agns tage 
also nach (25) 
P+rlr+lpe +9 )s+lP +apalt=dV 
Für die Ableitungen g,, 94 erhalten wir: 
mM ri tu +, Mita Hd" 
und daher mit der Abkürzung [+ «(pP +) "”’=y(p,g): 
Y+rm=lituP td) PrVv(Bd; Pa t+49,= -2pqay(P,d); | 
P+4m zitat) lvo 
Setzen wir diese Ausdrücke in (27) ein, so erhalten wir für y=#0'') 
vd tae+P) pr 2pasti ta +N-h=0 . . . 0). 
Damit hat sich ergeben: die allgemeinen Integralflächen der Grundgleichung (Il) 
der ebenen kompressiblen Potentialströmung sind die Relativminimal- 
flächen des Variationsproblems ö0=0, dessen Relativoberfläche durch die 


2y-te Potenz der Schallgeschwindigkeit im ebenen kompressiblen Strö- 
mungsfeld gegeben ist. 


$ 3. Normalformen. 


Durch das im vorhergehenden abgeleitete Ergebnis ist der invariante Charakter der 
Grundgleichung (1) in Evidenz gesetzt worden. Es ist jetzt leicht, Gl. (1) mit geometrischen 
Hilfsmitteln auf bestimmte Normalformen zu transformieren. Dazu gehen wir wieder von der 
allgemeinen Parameterdarstellung aus und schreiben überdies u=«w', v=u?’. Sodann be- 
trachten wir die nichtsingulären quadratischen kovarianten Tensoren '?): 


10) Vgl. (8) S. 3. 


11) Für v=0 folgt durch Differentiation pr +4s=0, ps+qt=0, d.h. rt—s?=0, man hat also die Torsen 
- als Lösungen auszuschließen. 


12) Vgl. (8) S. 265, 266, 268. 
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Iıı u (a hi, F bhi,), gs u (a hi; Fr bh;,) dene (chi, +dh‘,) er Iaı ‚ Y32 (ec his + dh;,) (31) 
._ ak 1, __(aB,E,,6) 2 - un 
hi  — (e, L, 5 y,) R ka ‚= (e, e,. e,) a,p : R s R ; 5 5 . (32). 


Für ihre kontravarianten Komponenten gilt (a,ö=1,2; mod 2): 


* % h* 
Ja+ı,arı *12_ YJıa * Ia+rı1,as 1 


—— ( F N h aa — 
Iga 3] I a3 Ihe] 


u, 


Y 


rau —_ Kethatı kr— ki. 
kas| ’ kas 
Für H*=0 werden die Tensoren gs und küs proportional '*). Entsprechend lauten die in- 
varianten Formulierungen der Grundgleichung (1): 





ERIEE geht kei = isch =0] ut. te 





Wegen |gas\|kas|=50 können wir die relativisotropen Kurven '*) als Parameterkurven ein- 
führen und erhalten: 


ent =eir—=0, Gb ach BU: ==... ... (EB. 


Dafür können wir nach (32) schreiben 





2. 0, X2)=0, HE: - .» »» . nn 





und erhalten das Resultat: die im Überschallgebiet reellparametrige Normal- 
form (36)*) der (hyperbolischen) Grundgleichung (l) bestimmt ihre Inte- 
gralflächen in relativisotroper Parameterdarstellung. 
Von den relativisotropen Parametern gelangen wir zu den relativisothermen durch die 
Transformation: 
er en 0. dur du? ” TER Pr 
vV—=u tim, u - - : Img agdam SP el. 


In diesen gilt: 


= = Aulies = m ER = 
nel, nei Fr 7 Mu + hr) ee 


Dafür können wir nach (32) schreiben: 
(Eins Kira) + (Kuss &ı > Eo) =0 bzw.: AYKELX) | re 2) 


und erhalten das Resultat: die im Unterschallgebiet reellparametrige Normal- 
form) 





|1:6,x3): 0, u 920, 9,0 u 





der (elliptischen) Grundgleichung (1) bestimmt ihre Integralflächen in re- 

lativisothermer Parameterdarstellung. 
Wegen |h2s|=#0'”) können wir auf x Asymptotenparameter einführen. In diesen gilt 
= hihi =t .. . rl 


und wir erhalten die Normalform 





IM=h,=0, 9=0]... 





Schließlich erhalten wir für Relativkrümmungsparameter '*) die Normalform 





| ga + ga, =0, er rear ie 





13) Die quadratischen Differentialformen 2, = 9, ;du“d u?, = hr zdu“d u, 23 kr, zdutdu ? genügen der Re 
lation K* 2, — 2? H* 23 — 23 —=0, worin K* die relative Gaußsche Krümmung darstellt! 
1%) Vgl. (8) S. 1, 6; es handelt sich also um die Nullinien von 2, —=V. 
15) Vgl. (1) Kap. III, $ 2, 130 bis 133. 
16) Vgl. (15). 
17) Vgl. (11). 
18) Vgl. (8) S.3; die Gleichung der Relativkrümmungslinien lautet: bdw +f(d— aJdudr —cdv?=V, es i-t also 
(a—d)?+4beZ=0 vorauszusetzen. 


- 


‘* 
h 
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Umgekehrt bestimmen natürlich die Normalformen (42) und (43) der Grundgleichung (1) ihre 
Integralflächen in Asymptoten- bzw. Relativkrümmungsparameterdarstellung. 


S$S 4. Beziehungen zur Lieschen und Weierstraßschen Theorie 
der Minimalflächen. 

Die Normalform (36) der Grundgleichung (1) führt in der Theorie der gewöhnlichen 
Minimalflächen bekanntlich auf den Lieschen Satz über die Darstellung aller Minimalflächen 
als Schiebflächen isotroper Kurven, deren integrallose Darstellung die integrallose Darstellung 
der Minimalflächen ermöglicht. Für Relativminimalflächen gilt nach A. Duschek und L. Ber- 
wald'®): Relativminimalflächen sind dann und nur dann Schiebflächen relativisotroper Kurven, 
wenn die Eichfläche e (uw, ») des Variationsproblems von zweiter Ordnung ist, bzw. das Variations- 
problem vom Typus: 

r 3 
>| \ ae B;kPpipr d u d ® = BB; B;;=const, |B;r =F0, i, k: 1,2, 3 ö . R (44) 
# i.k=1 


. | i er h : 
Wegen y=F7 scheidet also das Variationsproblem (5) für eine Anwendung dieses Satzes aus. 


Die Normalform (40) führt in der Theorie der gewöhnlichen Minimalflächen bekanntlich 
auf den Weierstraßschen Satz über die Darstellung der auf isotherme Parameter be- 
zogenen Minimalflächen durch Potentialfunktionen. Eine Verallgemeinerung dieses Weier- 
straßschen Satzes verdankt man A. Haar’). Man kann die Eulerschen Gl. (18) als 
Integrabilitätsbedingungen des Systems 

Ber Be r E 


auffassen und erhält durch den Vektor {=X(w,v) die zu x adjungierte Extremalfläche des zu 
d\\P(p)duder—0 adjungierten Variationsproblems ö | \ Piö)dudr=0*"). Für adjungierte 
Extremalflächenpaare x und X gilt”): 


I), (de +4, (dar)? +4 (de?) = ausdudu’= 1, (dE”+i,(de? +4, (da), a ‚B=12 (4, 


wobei die Koeffizienten 4, und 4, durch?) 


; DPP,,®.®,. Pr A PS, BD, 0, Dr 
on DE Bee ee =] + Ba F 2188. ,. 
PıPzPs Pk PıP2 Ps db; k 
gegeben sind. Sind «', u? (in (46)) „quasiisotherme* Parameter, so gilt 
BE AT TI I 


oder ausführlich geschrieben: 
: Or! 2 * 0.r?\° , dx’ 2 F Ox'\? 2 Ri) x? 2 ; dx 2 
hy 5 =) + 1 + Mr Sy hy FE zn .(y = + il, a ’ | 
| it ker a PR: 


Od dd , darda” 
4 


7 s+ 4, + 4, = 
der du gu du? 


—=() 


In derartigen Parametern reduzieren sich die Differentialgleichungen (45) auf das System °*) 


Ri) x ; d ae! d r: er ; N) r? Ri) x’ a ; N) x’ 
du? 1du du dur’ du? du 
ö 8 RE 
9 - ; 02 09° ; 02° 08 j o.x° 
du Adu” du! du” du 79 


18) Vgl. L. Berwald: Über adjungierte Variationsprobleme und adjungierte Extremalflächen, Monatshefte für 
Math. u. Phys. Bd. 38 (1931), S. 89 bis 108, insbesondere 105, 106; (8) S. 7. 


2") Vgl. A. Haar: Adjungierte Variationsprobleme und adjungierte Extremalflächen, Math. Ann. Bd. 100 (1928), 
S. 481 bis 502, insbesondere S. 496 bis 500. 


2) Vgl. H.Gericke: Zur Differentialgeometrie von Flächen im »-dimensionalen euklidischen Raum, adjungierte 
Extremalflächen, Math. Z Bd. 46 (1940), S. 408 bis 459, insbesondere S. 428, 429; (20) S. 484, 485; (19) S. 92 bis 95. 


22) Vgl. (20) S. 494, 495. 
22) Vgl. (20) S. 495. 


24) Vgl. (20) S. 498, 499, 
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woraus z. B. für die Koordinaten .* der Fläche I die Gleichungen 





2 
Ok w 
kl» = —, k= 8 3 : ’ ° (51) 
*s 0 u“ 


N ck N rk Wi 
r FE Be E - Ri 04; 
)r Sc® + 4, + Ars „—V) in. = 
k kı d u' ka d u ; kı 


N 


du" 











folgen, an Stelle von 1x*=0 im Falle gewöhnlicher Minimalflächen. 


"Der Zweck der vorstehenden Ausführungen war, die geometrische Natur der in der (ras- 
dynamik grundlegenden Differentialgleichung (1) sowie ihre aus geometrischen Zusammen- 
hängen in begrifflich einfacher Weise hervorgehenden Transformationen zu beleuchten. Dabei 
sind Legendresche und sonstige Berührungstransformationen mit Absicht vermieden worden. 
In einer späteren Untersuchung soll versucht werden, die allgemeine Integration der Grund- 
gleichung (1) mit Hilfe der hier eingeführten geometrischen Hilfsmittel anzugreifen. 241 


Der schräg gestellte Ringspalt. 
Von Kurt Voelz in Braunschweig. 
Bericht aus dem Institut für Aerodynamik der Luftfahrtforschungsanstalt Hermann Göring Braunschweig *) 
Verf. zeigt, wie man den Einfluß der Schrägstellung der beiden Platten 
eines Kreisringspaltes berechnen kann und beschreibt eine für entsprechende 


Messungen geeignete Versuchsanordnung. In dem durchgemessenen Fall 
bestätigen die Messungsergebnisse im wesentlichen die Rechnungsresultate. 


A. Theoretischer Teil. 


1. Der parallele Ringspal. Für einen Ringspalt mit parallelen Begrenzungsflächen, 
d. h. einer Anordnung, bei der zwei kreisrunde Platten parallel gegenüberstehen und durch 
eine in der Mitte der einen Platte befindlichen Öffnung Öl einfließt und radial nach allen 
Seiten ausströmt (vergl. Bild 1), sind für laminare Strömung die Beziehungen zwischen 











2 dem Druck in der Mitte der Platten p;, 
ey — der Durchflußmenge in der Zeiteinheit @, 
DZ er .. 

der Zähigkeit des durchfließenden Ols «, 
l 1 
| | der Spaltbreite h 





IR, und der vom Ol auf die der Einflußöffnung gegen- 
überliegende Platte ausgeübte Kraft P 
-R; Bi 
Bild 1. 


bereits bekannt '!). 


Als Grundlage zu ihrer Berechnung dienen die Navier-Stokesschen Gleichungen, in 
denen wegen der großen Zähigkeit des Ols die Trägheitsglieder gegenüber den Zähigkeits- 
gliedern vernachlässigt werden können: 

gende — unit. . -. 2. 2 2. 

Wegen der Rotationssymmetrie fallen für den vorliegenden Fall sämtliche Glieder mit 

vp sowie sämtliche Differentialquotienten nach 9 fort. Ferner kann angenommen werden, 


daß die Komponente senkrecht zur Hauptströmungsrichtung so klein ist, daß sie samt ihren 
Differentialquotienten gleich O0 gesetzt werden kann. Aus dieser Annahme folgt: 


BE a a a Er 


so daß also ein Druckgradient senkrecht zur Hauptströmungsrichtung nieht vorhanden ist. 
Da unter den angegebenen Vernachlässigungen 
Wr, 
rot rotd=— —, , 
02° 
so erhält man 


*) Inanderer Form vorgetragen auf der wissenschaftlichen Tagung der Gesellschaft für angewandte Mathematik 
und Mechanik in Braunschweig am 4. Mai 1940. 


1) Vgl. Gümbel-Everling: Reibung und Schmierung im Maschinenbau. Berlin 1925, Seite #3 bis H. 
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=u NE ee a a we er RR 
or ' d2 


Man kann wegen (2) sofort über 2 integrieren und erhält wegen der Grenzbedingung, 
daß v, für z=0 und z=h verschwinden muß: 
1 dp 
h.z)- ae A ie 


v.„= (2? 12 
Aue 7 dr 


Die gesamte durch eine Ringfläche hindurchfließende Menge beträgt dann 
h 
; rl’ dp 


6u dr d). 


=2n r\ v„d 


0 


IN 


Daraus folgt durch Integration für den Druck unter der Bedingung, daß beim Austritt 
aus den Platten Atmosphärendruck herrscht, daß also, wenn der Halbmesser der Platten mit 
R, bezeichnet wird, der Überdruck p=0 für r=R, wird: 

buW R, 
p ch? n ee 

Wird der Innenradius der Platte, durch die das Öl eintritt, mit R, bezeichnet, so beträgt 
der Druck, unter dem das Ol eintritt, 


(6). 


buy, R, 


Ppi- r h°? n R, ’ (4). 
Man kann somit auch schreiben 
R, 
In — 
r 
p=Pi R (S). 
In „” 
R, 


Die Kraft P wird hervorgerufen durch den Druck p; auf die Kreisfläche mit dem Halb- 
messer R, und den veränderlichen Druck p auf die Ringfläche: 
R, 


P=n Rp+\2ap rdr. 
R, 
Durch Integration folgt daraus: 
>? >. 
P=3 2 Ep Re A 
R, 

2. Die Berechnung des Einflusses der Schrägstellung. In der vorliegenden Arbeit 
ist rechnerisch und experimentell untersucht worden, wie sich die obigen Beziehungen 
ändern, wenn die beiden Begrenzungsplatten nieht genau parallel, sondern unter einem kleinen 
Winkel ® gegeneinander geneigt sind (vgl. Bild 2), und 
inbesondere, wie der Druck p; und die mittlere Höhe Ah, bei 
konstantem P, « und @ bei Schrägstellung sich ändern. 

Eine direkte Lösung der Navier-Stokesschen 
Gleichungen konnte nicht gefunden werden, da jetzt außer 
der v,-Komponente auch noch die vy-Komponente in Um- 
fangsrichtung auftritt. Jedoch kann man eine Lösung ?’) 
durch Reihenentwicklung erhalten, wenn man 

1. sich auf kleine Neigungswinkel 9 beschränkt. Auch 
für diesen Fall läßt sich die v,-Komponente mit dem gleichen 

u! Recht vernachlässigen wie für die parallelen Begrenzungs- 
flächen. Die Geschwinäigkeitsverteilung über die z-Koordinate wird dann mır wenig von der 
für #—=0 berechneten abweichen, man kann also 

2. für die Geschwindigkeitsverteilung über die z-Koordinate ein parabolisches Profil an- 
nehmen, das in Analogie zu (4) gegeben sei durch 








1 
=; „(e- h,)gradp. 


?) Die Lösung ist nach einem Vorschlag von Prof. Busemann erhalten worden. 
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Man erhält daraus für die mittlere Geschwindigkeit 


h 
® 1 
D= h \pde=-— 


h? 
12 gradp NT a re 


Jetzt sind sowohl d wie auch p nur noch von den beiden Variablen r und 9 abhängig, 
man hat es mit einem zweidimensionalen Problem zu tun. Die Stromstärke ist dann gegeben 
durch 


q=hv. u SR Re ME De er Te |; ° 
Da im Inneren des Ringes weder Stromlinien entstehen noch enden können, muß 
j f a? 
divgq= diviz,, ER" EEE EEE | 
oder 
33h’ gradhgradp+h’graddivy=0. . . 2.2.2... (8) 


sein. Die Lösung soll durch Reihenentwicklung nach # erfolgen. Da man auch die Gesamt- 
kraft P und die durchfließende Menge 9 kennen will und beide von dem Vorzeichen von # 
unabhängig sind, kann man die Reihe frühestens nach dem quadratischen Gliede abbrechen. 
Bezeichnet man jetzt mit a den Gradienten der oberen Platte 


a Se En ae a ee A 
und mit h, den Abstand in der Mitte der beiden Platten, so ist 
h=h,+tua-r 
und 
u g#) 
Mit diesen Werten wird, wenn man noch p=p,+p, +p: +... Setzt, 


divgrad p, + div grad p, + divgrad p,—+.... 


Salar 
' 


3a 3a 3 ) 
u gradp, h gradp, + h gradp, +... 


Nun ist unter der Voraussetzung, daß ® klein ist, auch a klein, und man erhält, indem 
man Glieder gleicher Größenordnung gleichsetzt: 


a ee an et 
’ 3a 

divgradp, = h a le a EI N 
. 3a 3alar) m 

div grad p, = h gradp, + n: a 3 7 (17). 


Bezeichnet man nun mit i und j die Einheitsvektoren in Richtung der Zylinderkoordi- 
naten r und g, so ist, da für kleines ® gesetzt werden kann: sin®=z#, 


gradh=a=iPcosp —jPsing. 


Der Winkel p wird dabei von der Stelle aus gezählt, an der der Öffnungswinkel am 
größten ist. p, ist die Lösung für den parallelen Ringspalt. Es wird also nach (8) 


1 1 . > 
gradp,=pi' z F R,' (18). 
n R, 
Damit wird 
N 3 cosg 1 
div grad p; = u a ai R,' 
n R, 


Die Lösung dieser Differentialgleichung wird als ein Produkt zweier Funktionen f, (r) 
und g,(g) angesetzt. Es zeigt sich dann, daß g,(p)=cosp sein muß. Die übrigbleibende 
inhomogene lineare Differentialgleichung weiter Ordnung für f,(r) wird nach der Methode 
der Variation der Konstanten gelöst. Da p denselben Grenzbedingungen genügen muß wie p,, 
so muß p,=0 werden für r—=R, und r=R,. Das führt schließlich zu 











Z. angew. Math. Mech. 








38 Voelz, Der schräg gestellte Ringspalt Bd. 21 Nr.2 Apr. 1941 
' BR, 
39 ) „(MinR,-RinR, „n,. MR "ir 
7.73 h Pi aaa R: R? y rin? R: R: > 
*. Mi : ’ 
R, 
Daraus erhält man 
N R, 
139 1 R:ln R,—- R?in R, RR "RR, 
gradp, =il Pi °0sg RR Inr —1+ RR 
a o In >. 2 ı 2 1 
R, 
ML: 
2:43 I, [RnR-RıR ,_RR ge 3 
TuUTER  Fyaaaalıs KR: — R: "ET 
R, 
Nun folgt aus (14) 
a=iPcosp —jPsingy, 
also 
ar=rVcosq. 
Somit folgt aus (17) mit (18) 
' R, 
u Re... gi RinR,— RinR, 5 RR "TR 5 " 
y — m — . oo . =: 2 OS) « 
ae rt I Le 3 Yade R? — R? at ze er u Eee 


In R, 


Die Differentialgleichung ist von der Form 


Ip, =p(r)+g(r)eos2p. 
Setzt man als Lösung 
pP: = f.(r)cos2y + 9, (r), 
so erhält man für f,(r) und g,(r) wieder inhomogene lineare Differentialgleichungen zweiter 
Ordnung, die wieder nach der Methode der Variation der Konstanten gelöst werden können. 
Unter den entsprechenden Grenzbedingungen wird 


9.9 1 [., 1, RinR,—RlhR\,, 1R-R 
2) niPi R, rinr- 6 + RR ) +5 | R. In» 
"iR "IR, 
l R' In R, BR} In R, R BR} R3 | R, 
6 R, RR |"R, 
In „, 
h 


ı 
Die Beziehung zwischen der Gesamtlast P und dem Druck p; ergibt sich dann aus 
Er Rs 
P=rapR+\ \ptptp)rdrdg. 
oe Rı 


Die Summe aus den ersten beiden Summanden ist aus Gl. (9) bekannt. Da nun 
237 27 
\cospydg=lcos2ypd y=0, 
0" u 


an Ra 
so bleibt nur \ \g,rdrdg. Man erhält 
0 Rı 
an M—R o/1 ıR—-R RR R, 
p_ S ? EB Rd p) »2 rg u ı #3 2 
la: ice 00 a 7 2; nalen ir Se u ur 5 
In 7 In R. ö 


Die durchfließende Menge Q ergibt sich aus (10) und (11) zu 


in 


1 ri dp 
— 3 » ‘ 
Ze TAUFTAK 
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« 


4 ) BE; " 5 
Setzt man p=p,+tpı +p, und p, =fi- n ©0859 und 9, =f3 m: cos2g +5: g», so wird 


I), 


unter Vernachlässigung der Glieder höherer Ordnung, wenn man noch (19) und \cos’q 
u 


dgp=nr beachtet, 


dp dp af! dp 
2 5) «43 ® 2729? y s 3; “2 92 I1L37r 92 0 
[) 12 2arh: y +2n?9°’rh, p +3nh,r:: 2; ır’h,d d = 
oder ausgerechnet: 
ab 1 [,„3@e/i1R-R ,RR_®, 
GP RI TImlTE, RR ""R-E"R, 
nm nn 
Werden jetzt während der Versuche P und konstant gehalten, so ändert sich p; nach 
R, 
, br ri 9@li „ m, 1R-R BB,R, 
er u - saerreatrg r R-R"R 
j . In RB f 
1 


Da bei der Berührung der Platten #/h, stets den gleichen Wert, nämlich I/R,, annimmt, 
so folgt, daß in diesem Fall p;/P auch unabhängig von der Höhe wird. Der Druck sinkt 
also, wenn man von der Parallelstellung bis zur Berührung der Platten übergeht, stets um 
den-gleichen Betrag unabhängig von der ursprünglichen Breite h, des Ringspalts ab. 

Um die Änderung von h, zu erhalten, kann man den Wert von P/u @ berechnen. Löst 
man nach h, auf, so erhält man 

g I € 42 2 2 > 
PL Ri] PB... & (R+my 9 MR Ra)l. 
fi : gs m\a\" R:—_R! R,,| 

Die Gleichung kann am einfachsten durch sukzessive Approximation gelöst werden, 
indem man zunächst %k, für #=0 berechnet und mit diesem Wert in die Klammer eingeht. 

B. Experimenteller Teil. 

1. Die Versuchsapparatur. Zur Prüfung der Rechnung wurden Versuche mit Platten 
mit den Halbmessern A, =2cm und R,=5em unternommen). Die Apparatur hierfür ist 
in den Bildern 3 und 4 gezeigt. Die beiden waagerecht liegenden Platten von je 6 cm Stärke 








Dem - 














3ild 4. 





Bild 5 (links). 


befinden sich in dem Topf T. Die untere Platte hat in der Mitte eine kegelförmige Aus- 
bohrung von 4 em Durchmesser. Das Öl fließt von einer durch einen Gleichstrommotor an- 
getriebenen Wendelpumpe durch die untere Platte hindurch in diesen kegelförmigen Trichter 
und strömt von da radial nach allen Seiten aus. Nach dem Ausströmen sammelt es sich in 
dem Topf T und wird von da durch eine Saugleitung zur Pumpe zurückgeführt. In die 
Druckleitung war zur Bestimmung der durchfließenden Menge eine Meßstrecke mit zwei An- 
bohrungen eingeschaltet, die mit den Schenkeln des Quecksilbermanometers @ verbunden 
waren. Da der Druckabfall im Spalt proportional « @ ist und u wegen der nicht konstant 
zu haltenden Temperatur während des Versuches sich stark änderte, war es bei der Durch- 


3) Die Anregung zu den Versuchen gab Herr Dipl.-Ing. Eieke. Für Ratschläge bei ihrer Durehführung ist der 
Verf. Herrn Dr. Kerris dankbar. 
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führung der Versuche nötig, die Durchflußmenge durch Drehzahländerung des Motors nach- 
zuregulieren. Dies gelang sehr leicht, da auch der Druckabfall in der Meßstrecke der Druck- 
leitung proportional .@ ist und man so nur darauf zu achten brauchte, daß die Höhe der 
(uecksilbersäule stets konstant blieb. Durch ein Thermoelement wurde auch die Temperatur 
des Öls gemessen, und daraus konnte auf dessen Zähigkeit geschlossen werden. 

Die obere der beiden Platten ist durch drei Säulen an den Eisenträgern E, E, befestigt, 
deren beide Unterstützungspunkte durch die Stellschrauben Sch gehoben und gesenkt werden 
konnten. Damit war die Neigung der beiden Platten in zwei zueinander senkrechten Rich- 
tungen einstellbar. Das andere Ende des Trägers E, kragt etwas über. Hier wurden an den 
in der Abbildung sichtbaren Haken schwere Betongewichte angehängt. Kleinere Eisen- 
gewichte neben den Schrauben dienten als Gegengewicht, um ein Überkippen zu verhindern. 
Die auf der Platte aufliegende Last P konnte durch Rechnung aus den bekannten Beton- 
und Eisengewichten und der Hebelarmlänge ermittelt werden. Das Eigengewicht des Trägers 
und der Platte wurde dadurch gemessen, daß die Platte auf eine Waage gesetzt und die 
beiden Arme wie beim Versuch unterstützt wurden. 

Auf der Oberplatte waren vier Meßuhren aufgesetzt. Aus der Differenz der Anzeige 
dieser Uhren konnte die Neigung der Platte ersehen werden. Die Spaltbreite h, war durch 
ihren Mittelwert gegeben. 

Schließlich befand sich an der trichterförmigen Ausbohrung in der Platte eine Bohrung, 
die mit dem Federmanometer M mit einem Meßbereich von 7,5 at verbunden war, an dem 
also der Druck p; abgelesen werden konnte. 

2. Durchführung der Versuche. Die Versuche wurden so durchgeführt, daß bei 
einer bestimmten Belastung P die Neigung durch Drehen einer der Schrauben Sch geändert 


wurde, anschließend war mit Hilfe des Widerstandes W die Quecksilbersäule auf den vorher 


festgelegten Stand einzustellen. Dann wurden die vier Uhren sowie das Manometer M abge- 
lesen. Die Neigungsänderung geschah nach beiden Seiten. 









































nn N ’ ’ PER RR 
Zahlentafeli. | | | 
Mm | 
un | P | m | 1 im 
Be. kg at mm Abbildg. 
1. 53,8 1.496 0,396 o 
I. 1008 | 2788 036 D | 
II. 191 414 | 028 0 Po En 2 a | or 
IV. 1956 5.459 0,236 | | Versuch Nr. T = « 
V. 2416 6715 | 0,208 | | A 
VL | 3416 | 675 | 025 | + u ————— + 4 Ho 7 
VI. | 1491 | 4144 | 0292 | X | | | = en 
| | | n„ nyM=* 
Pe 2 3 4 5 5 7 3 
Besondere Sorgfalt erforderte - ; 7000. 
die Bestimmung der Nullage der Bild 5. 
Uhren jeweils vor Beginn der Ver- 7 —— 
Bor ENTE s rm, VersuchWr I =» 
suche. Sie war gegeben durch den gr 
niedrigsten Wert, den die Uhren beim - ii Ir DH. 
oe nr er. ee: 
Kippen anzeigten. = br _ 
Im ganzen wurden sieben Meß- 2 - - 
reihen aufgenommen, die in Zahlen- 
tafel 1 wiedergegeben sind. Die 
beiden Meßreihen V und VI und | 
ebenso III und VII sind mit gleichen Id 1 
(sewichten, aber verschiedenen Durch: | T an, ann 5 DE 
flußmengen aufgenommen worden. a ———— - 74 + | 
Da, wie oben gesagt, die Durchfluß- FORR 10 | 
menge*) während einer Versuchsreihe nn ERREE RE 
geändert wurde, hat es keinen Sinn, | | 
sie hier anzugeben. Statt dessen ist N - | 
die Höhe h, bei Parallelstellung an wm ’ ar ö CHE. u 
gegeben. _ Bild 6. 


*%) Die Zähigkeit des Öls betrug 1,5 bis 3,5 Poise. Die Durehflußmenge konnte zwisehen 15 und 45 eem/see ge 
ändert werden. 
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3. Versuchsergebnisse. Die Ergebnisse sind in den Bildern 5 und 6 wiedergegeben. 
In Bild 5 ist die mittlere Höhe über dem Neigungswinkel aufgetragen. Die durchgezogenen 
Kurven sind auf Grund der gemessenen P und «u @ sowie auf Grund der im ersten Abschnitt 
gegebenen Rechnung gezeichnet worden. Man sieht deutlich den Abfall von h, mit wach- 
sendem #9. Die stärke Streuung der Werte ist darauf zurückzuführen, daß mit den Uhren 
kaum genauer als auf '/ı» mm gemessen werden konnte. 

Zudem standen sie häufig nicht still, sondern zitterten infolge des unregelmäßigen 
Arbeitens der Pumpe etwas. Dabei entsprechen die besonders stark herausfallenden Werte 
dem Kippen nach einer Seite, so daß hier sehr kleine Unebenheiten in der Platte mitgewirkt 
haben dürften). Die gerade Linie in Bild 5 gibt diejenigen Werte wieder, die beim Berühren 
der Platten gemessen werden müssen; die Streuung der Werte um diese Gerade kann als 
ein Maß für die Genauigkeit der Messung von h, angesehen werden. 

In Bild 6 ist der Druck p; über der Neigung ® aufgetragen. Die Werte passen sich 
(len theoretischen Kurven sehr gut an, nur in der Nähe der Berührung sind Abweichungen 
vorhanden. Hier ist im allgemeinen der Abfall von A, auch nicht ganz so groß gemessen 
worden, wie er errechnet wurde. Beides ist darauf zurückzuführen, daß die Reihen, die für 
die Berechnung benutzt wurden, um so schlechter konvergieren, je mehr man an den Be- 
rührungspunkt beider Platten herankommt, so daß hier noch ein weiteres Glied berücksichtigt 
werden müßte. 


C. Zusammenfassung. 
In der vorliegenden Arbeit ist eine Möglichkeit gezeigt worden, den Einfluß der Schräg- 
stellung der beiden Platten eines Kreisringspaltes zu berechnen; die Rechnung ist für einen 
bestimmten Fall (R,=2cem, R,=5cm) durch Versuche im wesentlichen bestätigt gefunden 


worden. 226 
5) Die Versuche waren zuerst mit nur ? em starken Platten durehgeführt worden. Dabei ergaben sich voll 
kommen unzureicheude Ergebnisse, die auf die Durehbiegung der Platten, die höchstens einige !/,00 mm betragen haben 


dürfte, zurückzuführen waren. 


Schubmittelpunkt und Querschnittsverwölbung dünn- 
.wandiger Träger unterhalb der Beulgrenze *). 


Von H. Neuber in Braunschweig. 
Bericht aus dem Institut für Festigkeit der Luftfahrtforschungsanstalt Hermann Göring. Braunschweig. 


Die Schubmittelpunktachse wird als zugehörig zu einem neuen Haupt- 
achsensystem gedeutet. Aus ihrem Zusammenfallen mit der Drehachse 
(auch bei elastisch-nachgiebiger Einspannung) wird eine Schwerpunkteigen- 
schaft der Verwölbung hergeleitet. Bei „wölbfreien“ Querschnitten ist der 
resultierende Wandungsvektor für jede Schnittlinie geometrischer Ort des 
Schubmittelpunktes. Ein zugleich bestehender Zusammenhang mit Flächen- 
und Spannungsschwerpunkt gilt näherungsweise auch für nicht-wölbfreie 
Querschnitte und führt auf ein einfaches Verfahren zur Ermittlung des 
Schubmittelpunktes. 


I. Energetische Hauptachsendefinition. 


Die grundlegenden Beziehungen für den Schub- oder Querkraftmittelpunkt des Balkens 
wurden von C. Weber!) gegeben; danach ist dieser Punkt Angriffspunkt der Querkraft für 
drillfreie Biegung (Formänderungsdefinition). Aus energetischen Betrachtungen leitete ©. Weber 
den Satz ab, daß bei gewissen Einschränkungen der Schubmittelpunkt zugleich Drehpunkt 
des Stabquerschnittes relativ zum Einspannquerschnitt ist, wenn reine Torsionsbeanspruchung 
herrscht’). Später gab E. Trefftz*) eine energetische Beziehung durch Nullsetzen der vir- 
tuellen Arbeit der Querkraftschubspannungen an den Torsionsformänderungsgrößen (bzw. um- 
gekehrt) und zeigte, daß der so definierte Punkt mit dem Weberschen Querkraftmittel- 
punkt identisch ist. 


*, Vortrag auf der GAMM-Tagung in Braunschweig 1940. 


1) C. Weber: Z. angew. Math. Mech. Bd. 4 (1924), S. 334 bis 338 und Bd. 6 (1926), S. 85 bis 97; s.a. A. u. L. Föppl!: 
Drang und Zwang, 2. Bd., 2. Aufl., München und Berlin 1928, S. 78. 

2) Dabei ist die Einspannenergie vernachlässigt und der betrachtete Querschnitt als starrer Endquersehnitt 
des Stabes gedacht; die auf diese Weise vernachlässigten energetischen Beträge sind jedoch bei genügender Stablänge 
vernachlässigbar klein. 


3) E. Trefftz: Z. angew. Math. Mech. Bd. 15 (1935), S. 220 bis 225. 
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Der Gedanke der energetischen Zerlegung des Arbeitsausdruckes läßt sich übrigens 
ganz allgemein auf den gesamten Spannungsvorgang des Balkens anwenden, was ich zunächst 
der grundsätzlichen Bedeutung halber kurz erläutern möchte. Bei folgerichtiger Null- 
setzung der virtuellen Arbeiten zeigt sich nämlich, daß beim Balken zwei Systeme von je 
drei Querschnittshauptachsen existieren, und zwar bildet die Schubmittelpunktachse zusammen 
mit zwei weiteren in der Querschnittsebene liegenden Achsen (ich bezeichne diese als „Schub- 
hauptachsen“) ein ähnliches System wie Schwerpunktsachse und Biegungshauptachsen. 

Ist x, y, z ein kartesisches Koordinatensystem mit der x-Achse in Stabrichtung, sind 
ferner die Komponenten des im Stabquerschnitt übertragenen Kräftesystems P,= a, Py= b,, 
P.=b,, M,=b,, M,„=a,, M,=a,, allgemein a,, az, bß, by (i,k, ß,y=1,2,3), so ist die 
Spannungsenergie für 5 = (reiner Zug, reine Biegung): 


r 
” Na;a; A;x ar er ae NE. 
ik 
Hierbei ist, da für #0 nur Normalspannungen wirken, 
1 B ; k 
A, = >E \ 0% x dF FE er re 


F 
Es bedeutet hierin 0, die für a;—=1, a„=0(i=#h) eintretende Spannung; sie ist also gleich 
1 1 1 ER e . 
F bzw. 7 (J.2 - Jy2y) bzw. 7 Iu:2 J,y)(für i=1 bzw. 2 bzw. 3) zu setzen (F Querschnitt- 


Hläche, I,=|e’d F,I,=\ydF,I,„.=\yzdF,I=1,l T},.). A;, ist daher die virtuelle 


F H 
Arbeit der Spannungen zu a;—1 infolge der Formänderungsgrößen zu a; =1. 

Man erkennt nun leicht, daß für A,=4A,=0 die «-Achse durch den Querschnitts- 
schwerpunkt gehen muß (Achse I), für A,„=0 die y- und z-Achsen Biegungshauptachsen 
(II und III) sein müssen?). 

Für a;=0 (Schub, Biegung und Torsion) wird die Energie 
Bu .-.:. 


4. 
”. 


wobei 


I | ER 4 i 
00x + 21 + — eiy rk „tr: Tr:))AF (4 





I: 

F 
gesetzt werden muß, da jetzt Schubspannungen 7,. und r,. 
hinzukommen. Man erkennt hieraus’), daß für B,=B,=0 
die Trefftzsche Schubmittelpunktsbedingung erhalten 
für kurze Träger wird, (der Anteil der Normalspannungen verschwindet 







Schubhayptachsen wegen 0,—0 für $=1) d. h. die x-Achse als Haupt- 
achse im hier definierten Sinne entspricht der Schub- 


TRREERE mittelpunktsachse (IV). Ferner müssen die y- und z-Achsen 
für B,=0 eine ausgezeichnete Lage einnehmen; sie 
seien als „Schubhauptachsen“ (V und VI) gekennzeichnet. 

a IE Ihre Riehtung hängt von der Länge des Trägers ab. Für 


1% lange Träger liegen sie noch fast parallel den Biegungs- 
Bild 1. Lage der beiden Hauptachsensysteme hauptachsen, für kurze Träger weichen sie von diesen 
beieinem Winkelprofil(I Schwerpunkt, ITund a 2. ’ - PER. 2 En e . - 
he wesentlich ab. Bild 1 zeigt sämtliche Achsen für ein 
punkt, V und VI Schubhauptachsen. W inkelprofil. 

Beide Hauptachsensysteme lassen sich andererseits auch als Tensorhauptachsen 
deuten, da A;, und Bgß, als Komponenten zweier symmetrischer Raumtensoren aufgefaßt 
werden können. 


II. Der Schwerpunktssatz der Verwölbung. 


Die Schubmittelpunktsachse, die so als eine von sechs Querschnittshauptachsen erscheint, 
ist nach der Weberschen Formel durch die Torsionsverwölbung & festgelegt. Sind y und z 
jetzt Biegungshauptachsen, so gilt für die Koordinaten yr und z7 des Schubmittelpunktes 


4) Das Verschwinden des Zentrifugalmomentes I u. ist also hier nicht selbst Definition, sondern Folge der 
energetischen Definition. 
5) Durch Einsetzen der aus der St. Venantschen Theorie folgenden allgemeinen Ausdrücke für die Spannungen 


@ 3 (bzw. y) 


und Erfüllung der Nebenbedingungen a; =, ba=l: z 3 (bzw. 7) bei den Spannungen mit dem Index 3 (bzw. ;). 
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Bf; =; 

yr = ®zdF, =: ie 3 N © 

Y1 3 r,\ T #2, \ J or 
F F 


Hierbei ist # der Verdrehungswinkel; /, und /, sind Hauptträgheitsmomente. Die hierbei 
auftretende Verwölbung & ist auf den Schwerpunkt bezogen, d. h. für die Torsionsschub- 
1 . . 
spannungen WE und r,. gelten die Beziehungen: 
’ gr’ 
ER _ 9° +92 r\ H.: 
G *9 dy Mer A 7 


Da die bei Torsion eintretende Drehung jedoch nicht um den Schwerpunkt, sondern um den 
Drehpunkt erfolgt, welcher nach dem Satz von Weber für den einseitig eingespannten 
Träger mit dem Schubmittelpunkt zusammenfällt, schreibe ich die obigen Gleichungen in der 
Form 


BR Be. ' 


1, IE, 9, BR lı 05 TER . 

G ua, ro zT), G Txz Is v(y Y1 ru Et 
Jetzt bedeutet & die eigentliche Verwölbung, die sich auf die Drehung um den Schubmittel- 
punkt bezieht. Ich möchte sie als „natürliche“ Verwölbung kennzeichnen®). Die Identität 
der Gl. (7) mit (6) verlangt 


“ - 
ee “ 
. Kar - 


FE a er 
Der Unterschied besteht zwar nur in einer anderen Orientierung der (uerschnittsbezugsebene, 
hat aber zur Folge, daß jetzt eine bemerkenswerte Eigenschaft der Verwölbung hervortritt. 
Setzt man nämlich & für &° entsprechend (8) ein, so erhält man zwei Gleichungen, welche 
nunmehr von der Wahl der Achsen unabhängig sind, näınlich 

\EydF=0, Ei. ne ne I MR 


F v F 
d. h. es besteht der folgende Satz: Die natürliche Querschnittsverwölbung be- 
wirkt ein Heraustreten der Querschnittsfläche aus ihrer ursprünglichen 
Ebene derart, daß derSchwerpunkt eines schmalen Trägerstückes, welches 
von einer verwölbten und einer unverwölbten Querschnittsfläche begrenzt 
ist, mit dem Querschnittsschwerpunkt zusammenfällt. 
E III. Elastisch-nachgiebige Einspannung. 

Dieser Schwerpunktssatz läßt eine weitere anschauliche Deutung zu. Denkt man sich 
den Träger elastisch nachgiebig eingespannt, wie es in den meisten praktisch vorkommen- 
den Fällen angenommen werden muß, so erscheint es plausibel, die im Einspannquerschnitt 
bei Torsion entstehenden zusätzlichen Normalspannungen s, welche die freie Torsionsver- 
wölbung zu verhindern suchen, dieser proportional anzunehmen, d.h. 


o=cC 


Irre, 


(10). 


Hierbei ist & zunächst noch auf einen beliebigen Drehpunkt bezogen. 

Die Spannungen o bilden ein Gleichgewichtssystem und werden sich nur im Einspann- 
querschnitt und seiner unmittelbaren Nähe bemerkbar machen. Aus den Gleichgewichts- 
bedingungen folgen dann die Beziehungen 


EydF=0, \EzdF=0, 7 0 ME 
F F F 

Die beiden ersten Gleichungen sind mit (9) identisch, d.h. wir haben aus der Annahme (10) 
für zunächst beliebigen Drehpunkt zwei Gleichungen gewonnen, die mit jenen identisch sind, 
die für den Schubmittelpunkt als Drehpunkt gelten; mithin folgt aus (10) unmittelbar 
das Zusammenfallen beider Punkte, wie es auch dem Satz von Weber entspricht. 
Die dritte Gl. (11) verlangt noch, daß der Mittelwert der Verwölbung gleich Null gesetzt werden 
muß. Mithin ist £ gleich der natürlichen Verwölbung & (vgl. Anm. 6). Man erhält daher 
folgende Beziehung, die als Definitionsgleichung der natürlichen Veıwölbung 
dienen möge: . 


\eyar A\rzar „irar . u 5 
F I F 


’ 


Ir 


. 
— 
u 


y 
1. 


6) Hier und im folgenden wird unter der „natürlichen“ Quersehnittsverwölbung die auf den Sehubmittelpunkt 
bezogene Verwölbung verstanden: ihr über den Stabquersehnitt genommener Mittelwert sei, um eindeutig zu definieren, 
gleich Null. Wie nachstehend gezeigt wird, gilt sie insbesondere für Träger, die entsprechend GI. (10) elastisch nachgiebig 
eingespannt sind (praktische Fälle); vor allem verschwindet die natürliche Verwölbung exakt bei den unten behan 
delten wölbfreien Querschnitten. 
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Wir haben gesehen, daß die plausible Annahme einer der Verwölbung proportionalen 
Normalspannung im Einspannquerschnitt zu demselben Ergebnis führt, welches C. Weber 
auf anderem Wege abgeleitet hatte. Da ferner Gl. (10) für £=0 von selbst erfüllt ist (ebenso 
(9) und (11))), so folgt der Satz: Das Zusammenfallen von Schubmittelpunkt und 


Drehpunkt gilt ohne Einschränkungen, fallsdie Verwölbung verschwindet’), 


Nachstehend werde ich zeigen, daß es bei dünnwandigen Trägern eine Reihe von Quer- 
schnittsformen gibt, bei welchen die Verwölbung exakt verschwindet. Zugleich 
läßt sich ein einfaches Verfahren zur schnellen Ermittlung des Schubmittelpunktes angeben, 
das auch für nieht-wölbfreie Querschnitte anwendbar ist. 


IV. Die wölbfreien Querschnittsformen. 


In der Literatur ist nur der Kreis- bzw. Kreisringquerschnitt als wölbfrei bekannt; ein- 
gehende Untersuchungen des Verfassers ergaben jedoch, daß außer dem Kreisquerschnitt noch 
eine unendliche Mannigfaltigkeit dünnwandiger Querschnittsformen verwölbungsfrei ist, ins- 
besondere beliebig zusammenhängende Bereiche. Durch Betrachtung jener Torsionszustände, 
bei welchen eine Spannungslinie mit einer Linie konstanter Verwölbung zusammenfällt, läßt 
sich bei Beschränkung auf die nähere Umgebung dieser Spannungslinie ein stetiger Grenz- 
übergang zum dünnwandigen wölbfreien Querschnitt erzielen (die betrachtete Spannungslinie 
wird Wandungsmittellinie). Die ausführliche Ableitung wird an anderer Stelle veröffentlicht °) 
werden. Hier sei nur auf die hauptsächlichsten Ergebnisse hingewiesen. Für die Torsions- 
schubspannung des dünnwandigen wölbfreien Querschnitts gilt die Formänderungsbedingung 


ER ER IDEE Arad (13). 


wo F der senkrechte Abstand der Randtangente vom Schubmittelpunkt ist, ferner an jeder 
Stelle einer Wandung i (Spannung r;, Wandstärke d;) die Gleichgewichtsbedingung 


an ER LINE SER RE ER SEEEERE EN 
Für einen beliebigen Schnitt durch den gesamten Querschnitt folgt 


a a ei 3 A 
i 


Demnach gilt für die einzelne Wandstärke 
KOOLZEER: 23. an. Teer mar were, 
und für einen beliebigen. Schnitt 


FD a a AR EN - 




















4 


[Azi3.2] 


A219.3 \ 


Bild ?. Quersehnittsformen, die bei beliebiger Wahl der Wand- Bild 3. Der Schubmittelpunkt des wölbfreien 
stärken wölbdfrei sind (Dreieck, symmetrisches Trapez). Quersehnittes als Angriffspnnkt des resul 
tierenden Wandvektors. 


?) D. h. ohne die in Anın. 2 erwähnten Einschränkungen. Würde man zunächst eine beliebige Lage des Dreh- 
punktes annehmen (Koordinaten in bezug auf die Biegungshauptachsen Up und ?7) und die zugehörige Verwölbung 
! Null setzen, so würde =’ *@29 — 7.2) +50 (So ist Konstante) folgen. Aus (9) ergibt sich dann unmittelbar z,, ap’ 
„Ir: d.h. es kann nur die natürliche Verwölbung verschwinden und es besteht der obige Satz. 


Ss), Luftfahrtforsehung, demnächst. 





YıınA 


ea, & 


m — sin N er 


rn A 


:h. 


41 





Yı1laaA 


Z. angew. Math. Mech. a : . . ’ Fi e 
Bd. er Nr. ? Apr. 1941 Neuber, Schubmittelpunkt dünnwandiger Träger unterhalb der Beulgrenze 95 


ine 















geschlossener 
Dreieckquerschnitt 


ß Spannungsschwerpunkt 
7 























4 
oo Flachenschwerpunkt 
Vz 
un N Schubmittelpunkt 





Bild 4. Beispiele für wölbfreie Querschuittformen. 


\ 6-0 Winkelprofil 








Aus (16) folgt unmittelbar die Gesetzmäßigkeit für d, r 
. = . . z E Spannungsschwerpunkt 

zweifach zusammenhängende Bereiche. Bei konstanteı 
Wandstärke sind es Tangenten-n-Ecke. Dreieck und o/ X 
symmetrisches Trapez sind auch bei beliebigen Wand- / Flochenschwerpunkt 
stärken wölbfrei (Bild 2). Aus (17) folgt, daß bei Deutung pP 
der Wandstärken als Vektoren in Richtung der Wand- W, 
tangente der resultierende Vektor kein Moment in be- A, Schuömitenunkt 
zug auf den Schubmittelpunkt haben kann, also [0° N 
geometrischer Ort desselben ist. Die resultierenden °e. gr B. ERIEET 
(7 . . 8 . . 1? 1 5 Schubmitte tt, Flächense ‚p 
Vektoren zweier Schnittlinien, die z. B. mittels Seileck punkt und Spannungssehwerpunkt für Drei 
bestimmt werden können, schneiden sich daher im Schub- eck und Winkelprofil. 


mittelpunkt (Bild 3). Bei Anwendung auf jeden beliebigen Schnitt ist diese Be- 
dingungzugleichnotwendig und hinreichend für beliebig zusammenhängende 
wölbfreie Querschnittsformen. So ergeben sich unendlich viele Möglichkeiten der 
Formgebung (Beispiele zeigt Bild 4). 


V. Näherungsverfahren für den Schubmittelpunkt auf Grund des 3-Punkte-Satzes. 


Eine weitere Beziehung folgt aus (13) durch Integration für zweifach zusammenhängende 
(Querschnitte: 
Dryds Pryds $rzds Przds 


- : : = (18). 
Pris Prds Prds prds 

Hierbei ist „, z ein kartesisches Koordinatensystem mit dem Schubmittelpunkt als Ursprung. 
Die linken Seiten der beiden Gleichungen stellen die Koordinaten des Spannungsschwer- 
punktes dar (Wandungsmittellinie mit Belegung r), während die rechten Seiten die 1,5fachen 
Werte der Koordinaten des Flächenschwerpunktes liefern. Danach gilt der folgende 3-Punkte- 
Satz: 


„Der Schubmittelpunkt der zweifach zusammenhängenden wölbfreien 
Querschnitte liegt auf der Verbindungslinie des Spannungsschwerpunktes 
mit dem Flächenschwerpunkt, und zwar über den Flächenschwerpunkt 
hinaus in einem Abstand von diesem, der gleich der doppelten Entfernung 
beider Schwerpunkte ist.“ 


Bild 5 zeigt die Lage der Punkte für ein Dreieck. 


Die Anwendung dieses Satzes auf nicht-wölbfreie Querschnittsformen führt zu einem 
schnellen Näherungsverfahren für den Schubmittelpunkt. Sogar bei offenen Querschnittsformen 
ist dieses Verfahren brauchbar; dabei wird der offene Querschnitt durch eine Hilfslinie in 
einen geschlossenen Querschnitt verwandelt (Wandstärke der Hilfslinie ist Null, also ideeller 
Spannungswert unendlich; daher ist der Spannungsschwerpunkt zugleich Schwerpunkt der Hilfs- 
linie). Daß ein solcher Grenzübergang vom geschlossenen Querschnitt zum offenen durch Ver- 
schwinden einer Wandstärke in bezug auf den Schubmittelpunkt exakt möglich ist, zeigt 
insbesondere der Grenzübergang vom Dreieck zum Winkelprofil (Bild 5). Weitere Anwen- 
dungsbeispiele des Verfahrens werden an anderer Stelle folgen°). 219 
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Geometrische Darstellung eines Satzes von Hurwitz 
für Frequenzgleichungen fünften und sechsten Grades. 
Von Herbert Bilharz. 


Bericht aus dem Institut für Aerodynamik der Luftfahrtforschungsanstalt Hermann Göring, Braunschweig. 


Für reelle algebraische Gleichungen 5. und 6. Grades, die nur Wurzeln mit 
negativen Realteilen haben sollen, wird eine vektorielle Darstellung der 
dafür notwendigen und hinreichenden Hurwitz- Routhschen Bedingungen 
gegeben. Die Komponenten der Vektoren sind — soweit sie nicht selbst 
die Koeffizienten der “rleichung darstellen — einfach zu bildende Linear- 
formen aus diesen. 

Ein graphisches Verfahren führt mit Hilfe eines Kurvenblatts, welches eine 
bizirkulare Kurve #4. Ordnung und zwei zu dieser affin verjüngte Kurven 
enthält. rasch zur Verifizierung der Voraussetzungen des Hurwitzschen 
Satzes. 


1. Allgemeine Sätze über das Vorzeichen des Realteils der Wurzeln reeller al- 
gebraischer Gleichungen. 
Um die Stabilität der ständigen Bewegung eines dynamischen Systems von m Freiheits- 
graden (m >0, ganz) nachzuweisen, ist zu zeigen, daß eine reelle algebraische Gleichung n-ten 
Grades O<n=2m) 


fu («) Y Er ee a a a 


nur Wurzeln mit negativen Realteilen hat. Eine notwendige Bedingung für die Lösung 
dieser Aufgabe ist gleiches Vorzeichen aller Koeffizienten; für n =2 ist diese Bedingung auch 
hinreichend. 

Im folgenden sei a, als Diskriminante der (verkürzten) Bewegungsenergie positiv '). 
Dann gilt 

Satz 1 (Hurwitz)?’): Die sämtlichen Wurzeln der reellen Gleichung (1) mit a, >0 
haben dann und nur dann negative Realteile, wenn alle Determinanten 


a, 4, 0 0 0 0 0 
a, d, a, Ad, 0 0 0 

D,= 4, a, 4, dl, a, 4A, 0 (r 8... (2) 
Ayy-ı Aaron dy 


positiv sind; dabei ist au=0 für u>n und u<V. 

Genügt die Gleichung (1) den Voraussetzungen von Satz 1, so nennen wir sie im folgenden 
eine Hurwitzsche (Gleichung oder kurz eine H-Gleichung. 

Weil mit f„ (2) =0 auch f„() =". f„(lje)=0 eine H-Gleichung ist, sind auch alle 
Determinanten D,, die aus den D, durch Vertauschung von a; mit a„_; entstehen, positiv 
und umgekehrt. Ferner ist D„=a,„-D,„_, nach (2). Wir setzen 


ET A er ee A 


und nennen R,„ die Routhsche Diskriminante von f,(x). 
‚Jeder Gleichung (1) ordnen wir eine Substitutionsmatrix zu: 
y 0,0, 0... " 
>R,„= * BO EG RAR ARABIEN VE RE 
In In ww d u... (4) 
mit 


1) Für die numerische Rechnung normiert man (1) zweckmäßig auf aa —=1. Wir sehen hier von dieser Normierung 
ab, um die Homogenität in den später folgenden Formeln nicht zu zerstören. 

?, A. Hurwitz: Über die Bedingungen, unter welchen eine Gleichung nur Wurzeln mit negativen reellen Teilen 
besitzt. Math. Ann, Bd. 46 (1895). S. 273 bis 284. 
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b,=a,;, für v =O mod2 
Oo<vr<n 
a dy2, Alpe z. 
b,=- 228 für» =1mod2 (5) 
a, 2 
b„=0; an —=d für u >n. 


Aquivalent zu Satz 1 ist dann der folgende 
Satz 2 (Routh)°): (1) ist dann und nur dann eine H-Gleichung, wenn alle „Potenzen“ 


RR (a,) —=012,...,#0—-D 
positiv sind, also die Relationen bestehen: 
R® (a)=a,>0,R®(a)—b,>0, 
RY’(A)ERR (a ))=R(b,)= db, >0, 
b,b, bb 


NR (a,) N (b,) 1 r 3 >—V, 
Ferner gilt 
Hilfssatz 3: Fürv=1,2,...,n—1l ist 
1 * 
R(D,)= D ee TE BEN TE 


1 


Mit Hilfe von (6) erhält man unter Beachtung von a, = D,unda,-a, a,-a,—D, durch 
vollständige Induktion Satz 1 aus Satz 2. Setzen wir schließlich in (6) 


9,—R (D,) 
und betrachten diese 9, als Hurwitzsche Determinanten der Gleichung (n — I)-ten Grades 


n—1 a 


fu.) N) EEE IE N ae Bei. gr 


— n 1 


wobei die Koeffizienten die Zahlen b, der Matrix (4) sind, so erhalten wir den folgenden 
Satz 4: (1) ist dann und nur dann eine H-Gleichung, wenn die Gleichung (7) mit 
b, >0 diese Eigenschaft hat ?). 
Zusatz: Für die H-Gleichung (1) mit a, >0 sind alle b, in (5) mit O<r<n positiv; 
ebenfalls sind die b, positiv, welche analog den b, aus den Koeffizienten von f, (=) 
gebildet sind. 


2. Spezialisierung der allgemeinen Sätze auf Gleichungen 4., 5. und 6. Grades. 
Ziel der vorliegenden Arbeit ist, für n—=5 und n =6 eine geometrische Darstellung von 
Satz 1 zu geben und daraus eine graphische Methode abzuleiten, die geeignet ist, den rech- 
nerischen Aufwand abzukürzen. Wir erledigen zuvor den trivialen Fall n=4. 
I. Es sei n=4. In diesem Fall ist Satz 1 bekanntlich äquivalent zu 
Satz 5: Die biquadratische Gleichung 


f. (x) u a Re ae a 


„== 


ist dann und nur dann eine Hurwitzsche Gleichung, wenn sämtliche Koeffizienten und die 
Routhsche Diskriminante 


R,=a,a,a,—a,0a—ala, 


positiv sind. 

In einer (1, 3)-Ebene mit aufeinander senkrechten Achsen fassen wir a, und a, als 
Komponenten eines Vektors A auf, der mit der positiven (1)-Richtung einen Winkel a bilden 
1,4, — 4,4, 


2 und b,= a, 


ns . . . \s#f£ ( 
möge. Ebenso seien die a, und a, in (4) zugeordneten Größen b, = 
1 


3) E. J. Routh: A Treatise on the Stability of a given State of Motion, London 1877, S. 29 ff. Ein Auszug 
dieser Preissehrift ist in E. J. Routh: Die Dynamik der Systeme starrer Körper; Deutsche Ausgabe Bd. II, Leipzig 
1898, S. 207 bis 234, enthalten. 


4) Satz 4 findet sich in weit größerer Allgemeinheit bei I. Schur: Über algebraische Gleichungen, die nur 
Wurzeln mit negativen Realteilen besitzen. Z. angew. Math. Mech. Bd. 1 (1921), S. 507 bis 311: insbes. S. 311. Im 
folgenden wird im wesentlichen der Zusatz von Satz 4 benutzt. 
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als Komponenten eines Vektors ® mit der Neigung f gegen die positive (1)-Richtung in der 
gleichen Ebene eingezeichnet (Bild 1). Nach Satz 5 ist (8) dann und nur dann eine H-Gleichung, 
wenn die Ungleichungen 


ea De en 


L 


IV 


bestehen. _ 

Nach Satz 1 sind für n > 4 die Voraussetzungen notwendig, jedoch nicht mehr hinreichend, 
daß mit a, >0 die übrigen Koeffizienten und die Routhsche Diskriminante (3) positiv sind. 
Wir geben hierzu ein Beispiel: 

Die Gleichung fünften Grades 


e l r x l 
+ +2 +2 + cr t=0 


ist wegen 


re | | | 
"+5 :+20°"+28#° +50 +1l=(c+H) a" ->0+5 20° —ZicH+l) 
keine H-Gleichung, obwohl alle Koeffizienten und die Routhsche Diskriminante R,= 7/16 
positiv sind. 
Bevor wir zu den Fällen n=5 und n=6 übergehen, machen wir folgende Vorbemerkung: 
Aus der Analytischen Geometrie des Raumes ist bekannt, daß die Gleichung 


ee a ER Rare NA 
in einem rechtwinkligen, rechtsorientierten kartesischen (#, y. 2)-Koordinatensystem (Bild 2) 
einen Doppelkegel 8 mit der Spitze im Ursprung darstellt. Auf Hauptachsen £, n, £ trans- 
formiert lautet (9): 


x2-2--P=0 


2 2 


A 
L 
- Zz=0. 
») 


e. 


Irr 


Eine Schnittebene S: 

n = const. (> VÖ) 
schneidet 8 in einer Ellipse K vom Achsenverhältnis 1:y2. Die Projektionen von K (Bild 3) 
sind in Aufriß und Grundriß die Kreise K_, und K«, der Seitenriß ist eine Gerade unter 45° 
zur «-Achse und unter 135° zur z-Achse. Die Ungleichung 

Et ae, Er ee 
gilt für alle Punkte im Innern des elliptischen Kegels &. 

Il. Es sei n=5. Um die Voraussetzungen von Satz 1 für die Gleichung 


20,2 =0 a ee are ER 


‚= 


f,(&) 


nachzuweisen, müssen in den zu f,(x) und f,(#) gehörigen Routhschen Matrizen u. a. die 
folgenden Ausdrücke 
a,d,— d,d, a,a, - a,d, 
b,= E =. b,= - 
a, a, 
und 
a,a, — d,d, 


a, 


b,- 


nach Satz 4 positiv sein. 
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Für die Determinanten D, der Gl. (11) gelten dann folgende Äquivalenzen (z. T. unter 
Benutzung der vorangehenden): 


D, bu. D,»b, D,>O»ra,b, —a,b, >0: 
ferner muß 


D,=zaR(D)=a,.abb,-ab>0 .....2222.0..(12) 
sein. Wegen 


2 212 Ad, ,, 
a,a,b,b, -ab,=a,b, (a,b, —a,b,) — ; = >o0 
( 
i 
ist also 
D, a b,- 
Wir bilden in einem rechtwinkligen, rechtsorientierten kartesischen (1, 3, 5)- Koordinaten- 

system den Vektor 


j "=(y,,73,)5) 
mit den Komponenten: 
„=ab=4,.4,—4d,, 7: ab=a,0,— 44, %=ab=a,a,—a,a,. . (13). 


Nach Satz 1 ist also (11) genau dann eine H-Gleichung, wenn alle Komponenten von 7 
positiv sind und dieser Vektor nach (12) im Innern des Kegels 


8: W- )- =0 
liegt. 
Wir werden zeigen, daß sich die Berechnung der drei zweireihigen Determinanten in (13) 
umgehen läßt. 
In der zu (11) gehörigen Matrix (4) betrachten wir «a,, a,, a, als Komponenten eines 
Vektors A im (1,3, 5)-Koordinatensystem: ; 
A (a,,a,,4,), 
ebenso die den a,, «,, a, in (4) zugeordneten Koeffizienten b,—=a,, b,=a,, b,=a, als Kom- 
ponenten eines Vektors ® im selben System: 
B=(b,,b,,b,). 
Ist (11) eine Hurwitzsche Gleichung, so sind alle Komponenten von A und ® positiv. Wir 
bilden noch das Vektorprodukt 
vR p R ; : 
> EXB=E=Le,,6;, 6) 
mit den Komponenten 


G=4,4, —a,4,, .=4,4, — 40,4, v=04,— 4,4, 
Dann ist nach (13): 
= on v, =), > 0, c Yı <S 0 


und nach (12): 
o : C, cz > v0. 
III. Für n=6 habe die Frequenzgleichung (1) die Gestalt 
6 
=), =0 a ER 
) 1) 
In der zu (14) gehörigen Routhschen Matrix seien a,, a,, a, die Komponenten eines Vektors U: 
AU= (a,,a,,a,) 


und die ihnen zugeordneten Zahlen 


a,d, — d,d, 


die Komponenten eines Vektors ®: 
B—(b,.b,,b,) 
in unserem (1,3, 5)-Koordinatensystem. 


Für die H-Gleichung (14) müssen diese Komponenten alle positiv sein. Wir bilden noch 
das Vektorprodukt 
AXB=E =le,,r,,6) 
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mit den Komponenten 
.=ab,—a,b,, c,=a,b, — ab, c,=ab,—a,b,. 


Die Vorzeichen dieser Zahlen ergeben sich nach Satz 1 aus den Hurwitzschen Deter- 
minanten. Darnach ist 


D,=#e, >, D,=ub,>9; 


D,=a,- (a,b, -—ab)=—-ac,>0, 
also 
IE ee ER Dee 2 EA DER 
D,= -a,(b,c,+b,c,) >d, 
b, B _ © 
D: a,\e + E I = h, al=a,(,, —c)>0, 
also 
en nn -#-) ZUERETEEHEEE GEENELHEEE 5 
Ferner ist 
3 2 D, 2 
b, -(e, 6, )= (3 (b; + b, e,) 1, = cz D b, >, 
1 


und schließlich 
D,>0U re, >0. 


Hieraus ergibt sich aber aus (16) in Verbindung mit (15): 
e,<Vd. 
Es gelten also für den Vektor $% die Relationen 
und 
G—G5>0. 
Die Ergebnisse von II. und Ill. lassen sich daher zusammenfassen in 
Satz 6°): Die reelle algebraische Gleichung 
fünften Grades (11) sechsten Grades (14) 


mit a, >0 ist dann und nur dann eine Hurwitzsche Gleichung, wenn die Komponenten der 
Vektoren 


A = (a,,a,,a,) Y=(a,,a,,a,) 
B=(b,,b,,b,) B=(h,,b,, b,) 
E=-AXB 


in einem rechtwinkligen, rechtsorientierten kartesischen (1,3, 5)-Koordinatensystem folgende 
Vorzeichen haben 


) 9 © 
A| + - 4 








B| + + + n 
(0 + 
und überdies der Vektor &E nach (10) im Innern des Kegels 
a u A FE FE 
liegt, also die Beziehung besteht 
6, — 5 >0V re DE ent 


5) Nach Satz 6 brauchen die Vektoren A und ® nur bis auf Strecekungen (bzw. Verkürzungen) in sich bestimmt 
zu sein. Vermöge dieser Tatsache läßt sich die nachfolgend beschriebene graphische Methode selbst bei Koeffizienten 
verschiedener Größenordnung mit großer Genauigkeit durehführen. Z. B. gilt der Satz: (11) ist dann und nur dann 
eine H-Gleichung, wenn die Gleichung @-(wox5 + a3 x3 +as3x) + ?-(a X +a3x?+a,)=0 mit positiven « und 3 vom 
Hurwitzseben Typus ist. 
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3. Zeichnerisches Verfahren zur Entscheidung über das Vorzeichen des Realteils 
der Wurzeln bei Gleichungen 5. und 6. Grades. 


Nachdem in Satz 6 eine geometrische Formulierung des Hurwitzschen Satzes für n=» 
und n=6 gefunden ist, gehen wir jetzt dazu über, ein graphisches Verfahren für jenen Satz 
abzuleiten. 

Nach seiner Definition als Produktvektor steht & senkrecht im Nullpunkt auf der durch 
die Vektoren X und B aufgespannten Ebene &; durch (17) ist also die Orientierung von ® 
zu A gegeben; wegen (19) darf ferner 8 mit A nicht kollinear sein (im entgegengesetzten 
Fall wäre die Routhsche Diskriminante singulär). Die dritte Zeile in (17) legt denjenigen 
Teil des Kegels St fest, in dessen Innerem der Vektor & liegt. 

In der Projektionsebene tragen wir die Projektionen O A, ..., OBy der Vektoren A 
und ® ein (Bild 4) und sehen naclı, ob die Ungleichungen (17) erfüllt sind. Hierzu hat man 
das für die biquadratische Gleichung in Bild 1 angegebene Verfahren dreimal anzuwenden, 
nämlich auf die Vektoren O A, und OB, im Quadranten (1,3), auf die Vektoren OA, und OB, 
im Quadranten (3,5) und auf die Vektoren OA, und OB, im Quadranten (1,5). Die positiven 
Winkelrichtungen sind in Bild 4 durch Pfeilriehtung markiert. 
































ad 
(8) | [ (3) 
\ a \ 
je NN 
B. 
. par 
’ %; &, 
du e- = @) 88 
O-B,<a,<z Pr A Ss 
y-123 HP ’ 
Be u 
B; 
Bild 4. Bild >. 


Die Ungleichungen 


3 TI TE (= 1,2,3) 
sind dann notwendig und hinreichend dafür, daß die Vektoren A, Bund E=AXB die in (17) 
genannten Vorzeichen tragen. Das Erfülltsein von Formel (17) ist also mühelos aus der 
Zeichnung abzulesen. 

Es verbleibt die Ungleichung (19). Wir legen eine Ebene © senkrecht zur Achse des 
Kegels 8 in (18) durch einen vom Ursprung O verschiedenen Achsenpunkt M. Die Schnitt- 
ellipse K vom Achsenverhältnis 1:y2 ist im unteren Teil von Bild 5 in wahrer Größe ge- 
zeichnet. Ihre drei Projektionen sind die Gerade K, und die beiden Kreise K, und K... 
Die Schnittpunkte A’ und B’ der Vektoren A und ® mit © haben die Projektionen A’,,..., B%. 
Die durch A und B aufgespannte Ebene sei E (in Bild 5 schraffiert). Der Schnittpunkt des 
Trägers von € mit © ist — (". 

T sei eine Hilfsebene durch die Strecken O—- €’ und OM. Wir können und wollen 
dabei voraussetzen, daß 0 — C’ nicht mit OM kollinear ist, denn fällt € in die Kegelachse, 
so ist die Ebene € parallel zur Ebene & und der Fußpunkt F (s. u.) liegt im Unendlichen. 
Ist € senkrecht zur Achse, so fällt F mit M zusammen. Beide Fälle sind aber nach (17) 
auszuschließen. 

Wegen OM_LS& (nach Konstruktion von &) und O—C’ _LE (nach Definition von €) 
sind nach bekannten Sätzen der Stereometrie T_ LS und T_LE, und daher steht schließlich 
die Schnittlinie s (=Gerade durch A’ und B’) der Ebenen & und € senkrecht auf T. Wir 
zeichnen noch die Gerade durch — €” und M; sie schneide s im Punkte F. Da s__LT ist 
und — (”M in T liegt, ist s_L MF in F, also F der Fußpunkt des Lotes von M auf s. Ferner 


i 
ist wegen T_ LE auch O- " LOF. 











. 
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Diese Überlegungen gelten speziell für den Fall, dad O—C’ auf einer Mantellinie von 8 
liegt. Wir greifen denjenigen Fall heraus, für welchen diese Mantellinie in der Ebene T liegt 
und bezeichnen mit C* den Schnittpunkt der Mantellinie mit &, mit s* die Schnittebene der 
zu O(* im Ursprung 0 senkrecht stehenden Ebene und der Ebene &, schließlich mit F* den 
Schnittpunkt von s* mit der Verlängerung von MC* bzw. M— €. Alle diese Punkte F* defi- 
nieren eine bestimmte „Grenz*-Kurve @ in ©; also 


Definition: Die Grenzkurve @ ist der geometrische Ort aller Punkte F*, die Fußpunkte 
der Lote von M auf die Schnittgeraden der Ebene © mit den Ebenen senkrecht zu den Mantel- 
linien von 8 im Ursprung O sind. 

Denken wir uns für den Augenblick in der Ebene © ein rechtwinkliges kartesisches 
Koordinatensystem (X, Y) mit dem Ursprung M eingeführt (X sei parallel (3), die positive 
Y-Richtung sei MD) und den Radius der Kreise K, und K, mit o(> 0) bezeichnet. Das 
Oval @ ist eine bizirkulare Kurve vierter Ordnung vom Geschlecht Null. Ihre drei 
(reellen) Doppelpunkte sind der Ursprung M (mit imaginären Tangenten) und zwei Kreis-. 
punkte (mit reellen Tangenten). Die zwei ordentlichen Brennpunkte von @ sind die am 
Kreis um M mit dem Radius MD=M,O=oy2 invertierten Brennpunkte der Ellipse K. 
Wegen MC*.MF*=20? entsteht die Grenzkurve @ durch Inversion der Ellipse K an 
eben jenem Kreis, hat also die Gleichung: 


G: (X?’+-Y)% —-207”’—40X’=0 mit X’+Y?>0. 
Die Gleichungen der Kurven @ , und G, sind dann 
(+27 -40. 7? -40.X°’=0 mit £°+7°>0, 


mit den Koordinatenursprungspunkten M, bzw. M, und den Achsensystemen (X, Y) in der 
Aufriß- bzw. Grundrißebene, X parallel zur Richtung (3), Y in Richtung (5) bzw. (1). Die 
Einheiten auf sämtlichen Achsen sind dabei einander gleich und sonst beliebig wählbar‘). 
. . u ie ra — i y. me .. T . ._—— 

Analog wie oben ist schließlich 0 C* 1 0 F* |also der Winkel (C* O F* =5) und s*_ LM F* 


in F*. Liegt nun F außerhalb der Strecke MF* (also außerhalb @), so ist der Winkel 


(- CO FF)< —, also (- ("O F*)< (C*O F*); daher liegt (” innerhalb der Strecke MC*, 


also innerhalb X, d. h. der Vektor & liegt im Kegel ft. Liegt andererseits F innerhalb M F* 
(also innerhalb 6), so ist (—- CO F*) > (C*O F*); daher liegt —(” außerhalb M C*, d. h. der 
Vektor E liegt außerhalb des Kegels 8. Im unteren Teil von Bild 5 sind diese Winkel in 
wahrer Größe konstruiert. Zusammengefaßt ergibt sich der 


Satz 7: Wenn das vom Kegelachsenschnittpunkt M gefällte Lot die Schnittgerade s 
der Ebene € (durch die Vektoren A und ® aufgespannt) mit der Kegelschnittebene & im 
Punkt F trifft und dieser Punkt außerhalb bzw. innerhalb des Ovals @ liegt, so liegt der 
Produktvektor E=AXB innerhalb bzw. außerhalb des Kegels &. 

Damit ist also der geometrische Nachweis der Bedingungen von Satz 6 einfach zu 
führen. Wir sprechen ihn in folgender Regel aus. 


Regel: In einem Kurvenblatt, welches das Oval @ und seine drei Projektionen @,, @«, 
G@s—=Kg enthält, zeichne man die Projektionen der Vektoren X und B: OA,,..., OBg ein, 
und sehe nach, ob sie die Bedingungen von Bild 4 erfüllen (Nachweis der Vorzeichen- 
relationen (17)). Ist die Forderung (17) verwirklicht, so zeichne man zunächst den Aufriß 4’, 
und B’; der Schnittpunkte von A und ® mit der Ebene ©. Die Ungleichung (19) ist erfüllt, 
wenn die Gerade durch die Punkte A’, und B’, die Kurve @, nicht schneidet. (Diese 
Untersuchung kann man auch in der Grundrißebene durchführen; der Aufriß ist aber kon- 
struktiv vorteilhafter.) Wird jedoch diese Kurve geschnitten, so lote man die Schnittpunkte 
auf das Oval @ herunter. Man erhält die Gerade s. Vom Punkt M ist schließlich das Lot auf s 
zu fällen. Liegt der Fußpunkt F dieses Lotes außerhalb G, so ist die Ungleichung (19) 
dennoch erfüllt; liegt andererseits F innerhalb oder auf G, so ist (19) nicht erfüllt, also die 
Ausgangsgleichung (11) bzw. (14) nicht vom Hurwitzschen Typus. 

Aus dem Vorangehenden ergibt sich, daß für n=5 das Verfahren rein geometrisch ist, 
während für n=6 zuvor die Ausdrücke b, und b, — also zwei zweireihige Determinanten — 
zu berechnen sind. 229 


6%) Zur Geometrie der bizirkularen Kurven 4. Ordnung vgl. die ausführliche Darstellung in G. Herglotz: 
Über die analytische Fortsetzung des Potentials ins Innere der anziehenden Massen. Jablonowski Preisschrift XLIV:; 
Leipzig 1914, S. Tft. 
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Über eine Verzerrungsfunktion in der Nomographie. 


Von Hans Knobloch in Berlin. 


Die Darstellungen der Nomographie enthalten zum Eingang und zur Ab- 
lesung Teilungen, denen jeweils bestimmte Abbildungsgesetze zugrunde liegen. 
Es wird gezeigt, wie die Abweichung einer solchen Teilung von einer Ver- 
gleichsteilung. beispielsweise einer linearen. zahlenmäßig angegeben werden 
kann. Anwendung bei Doppelleitern. 


Einführung. 


Die zum Eingang und zur Ablesung dienenden Teilungen nomographischer Darstellungen 
werden aus Gründen der leichten Herstellung i. a. mehr oder weniger verzerrt aufgetragen. 
In vielen Fällen jedoch wird von der Person des Benutzers her die Forderung gestellt, daß 
die Eingänge und die Ablesungen bei graphischen Tabellen möglichst an linearen Teilungen 
erfolgen sollen. 


Nun lassen sich beliebige Funktionen natürlich nicht allein mit linearen Teilungen dar- 
stellen; die Anlage beispielsweise einer Doppelleiter läßt sich jedoch stets so vornehmen, daß 
die lineare Teilung als Optimum angestrebt wird. Da sich die folgenden Überlegungen 
wesentlich einfacher gestalten bei Verwendung spezieller zusamnıenfassender Begriffe, so 
sollen zunächst zwei Funktionen definiert werden: das Verzerrungsmaß und die Verzerrungs- 
funktion. 


Die: Grundbegriffe. 


Die Teilung a mit dem Schritt 9, = |a„+, —- a„|, welche durch die Abbildungsfunktion . (a) 
auf die (in « linear geteilte) Trägergerade x gebreitet wird, entfernt sich offenbar um so weiter 
von dem Charakter einer linearen Teilung, je stärker zwei aufeinander folgende Teilungs- 
intervalle 


Axrlan-,) = |la,) - 2 (au- )|; 4x (a) = |* (ay4 1) — x (a,) 


voneinander abweichen. Dabei kann die Differenz zweier aufeinander folgender Teilungs- 
intervalle, bezogen auf die Länge des einen Intervalles, als Maßstab für die Größe der Ver- 
zerrung gelten. Auf ein solches Maß der Verzerrung einer Teilung hat nun der Abbildungs- 
maßstab keinen Einfluß. wohl aber die Größe des Schrittes 9, = a„+,— a,|=const. Es 
läßt sich leicht einsehen, daß bei Verdoppelung der Schrittlänge auch der eben beschriebene 
Maßstab der Verzerrung auf das Doppelte anwächst, daß also das zu definierende Ver- 
zerrungsmaß der Schrittlänge ®, proportional sein muß. (Man denke sich z. B. auf dem 
Rechenschieber die Teilstriche zwischen 1 und 2, 2 und 3 usf. ausgelöscht.) Dies zu betonen 
ist deshalb wichtig, weil man sich zur Aufstellung praktisch brauchbarer Abschätzungen ge- 
zwungen sieht, von den Differenzenformeln zu Differentialbeziehungen überzugehen, d. h. ge- 
wisse Aussagen auf infinitesimale Bereiche zusammenzuziehen. Als „Verzerrungsmaß* V..n 
der Funktion x (a) sei daher definiert: 


Vo = lım a . ne; (a + da) > Ix (a) ® x” (a) 


=— :9 . ’ . . : ö i 
Ja>\ Ja Ix (a) x (a) q a (1) 


(x’ (a) und x’ (a) bedeuten hier stets die Ableitungen nach dem zugefügten Argument, deren 
Existenz vorausgesetzt wird.) Der Quotient aus den ersten beiden Ableitungen der abbildenden 
Funktion x (a) werde als „Verzerrungsfunktion“ vx.., bezeichnet: 


x” (a) B 


d Ei 
Ur” x' (a) rt da [log x (a)] 
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In Tabelle 1 ist eine Übersicht Tabelle 
über die Verzerrungsfunktionen ; 
einiger häufig auftretender Abbil- Abbildungs- | Verzerrungs- 
dungsfunktionen wiedergegeben. funktion funktion Bemerkungen 
Man entnimmt daraus u. a., daß (a) Uxım 
für einen Rechenschieber zur Lö- 
sung des pythagoräischen Lehr- . dan 
satzes (Teilungen 2 =.«?) ebenso { x R 2 
2a ( ge . ) 40 2, z=a® Bei @>0 ist die Abbil- 
wie für den üblichen Multipli- Ex j 27 
kations-Rechenschieber (Teilungen rn a. Wiek When, Se 
ge u a - die Verzerrung über alle 
x—=]log.a) das Verzerrungsmaß im ee en Me 
ganzen Bereich konstant ist, wenn wächst ferner mit dem 
man den Schritt #9, proportional « Betrage von n—1. 
wählt. 
Es ist nun nützlich, sich ei- 
nige allgemeine Beziehungen an- : 
v 8 ns ur ®, proportional a ist das 
zumerken. Für die Umkehrungs- ih en 
f ktio } &) ilt erzerrungsmäatk ım gan- 
unktion a=a(x) 8 zen Bereich konstant. 
(Multiplikations-Rechen- 


DR 
i 2.2 ee s 
Gm = nr (3). schieber.) 








Lineare Teilung. 


z= log a Bei «>60 ist die Abbil- 
dung nicht möglich. Für 


(Beispiele hierzu bilden die Zeilen 3 3 1 Verzerrungsfunktion im 
und 4 der Tabelle.) ganzen Bereich konstant. 
Eine weitere wichtige Be- Bei + a/2 ist die At 
ziehung findet man, wenn man ae FREE BR Fee 
rn . bildung nicht möglich, bei 
z. B. die Gesamtverzerrung einer Ei 
s : - @>0 ist sie unverzerrt. 
gegebenen Funktion $= (a) durch 
Zwischenabbildung auf einen Trä- 
ger x in zwei Teile aufspaltet. Die Teilung « werde durch «= .r (a) auf den Träger x ab- 
gebildet; die Ablesung erfolgt über = f(x) so, daß = (x (a)) = (a) resultiert. Dann ist 
3 Val?) 
B3 (a) =Vx (a) + Vila) x (a) = 
(@) N ( ‚+ (x) ( ) a’ (8) 
bzw. 


Ux(e)* d «A u Va (9) ° d pP .- 0"; (x) ° d.r 0 
. w (4). 
d (log x’ (a)) + d (log a’ (£)) + d (log pP’ (#)) =V 
Ganz allgemein besteht im Falle von » geschachtelten Funktionen x, x, ... für die 
also =” = = (#”+®P) mit r—=1,2,... n gilt — die Beziehung 
n—| 
N 2m («P+d)- da"+» + tin) (au) dx? —=0. 


vl 
Anwendung bei Doppelleitern. 

Es seien nun die an einer Doppelleiter auftretenden Verzerrungen betrachtet. Vor- 
gegeben sei eine Funktion $=P(a), welche durch die Abbildungsfunktionen # = x, (a) und 
x=x,(ß) auf dem (gemeinsamen) Träger x als Doppelleiter dargestellt werden soll. Mit 
P=Pß(a) ist auch die Verzerrungsfunktion vs(.), für die Gesamtverzerrung vorgegeben, welche 
nun nach (4) so in vx,(«) und vs) bzw. vr,» aufzuteilen ist, daß dabei die Verzerrungs- 
maße Vx,(« und Vx,(3) der beiden abbildenden Funktionen x, (a) und x,(ß) möglichst kleine 
Absolutbeträge annehmen. Man hat also nach (4) und (3) 

P’ (a) 
7 


N 
3 


u ı-. . , 
V3la) tx (a) Un) Pla = 9. Vx; (e) Vx,(3) U EEE 


a 
Da nur monotone Funktionen P(a) betrachtet werden sollen, kann man P’(a) im ganzen 
Bereich als positiv ansetzen. Für vx,«=vs(«) ist dann vx,(9=0, für vx,«) 0 dagegen ist 
vx,(3)50. Läßt man nun die Funktion vx,(«) von:.®;(«) zu immer kleineren Funktionswerten 
und schließlich zu vx,(« 0 übergehen, so wachsen gleichzeitig die Funktionswerte vx,(s) von 
Null bis auf endliche Beträge. Im Bereich 0< vx, (© < v3(a) (bzw. O< -- va,(e) < — 93(a)) EXi- 
stiert also ein Teilbereich, in dem sowohl vx,(.) als auch vx,(3) (für jede Stelle x bzw. «) in 
ihren Beträgen kleiner sind als ”3(«). Dies bedeutet, daß sich die Gesamtverzerrung einer 
darzustellenden Funktion stets aufteilen läßt in zwei Teilverzerrungen von geringeren Be- 
trägen. 
Die Aufteilung der Gesamtverzerrung möge nun so vorgenommen werden, daß an jeder 
Stelle x des Trägers die Beziehung gilt 
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r er Ve. a= Ve). 
Wegen (5) ist dann 
VB (a) 


I ’ 
1+ d5 ß (a) 


Dadurch ist nach (2) auch die Abbildungsfunktion «, (a) festgelegt mit 


®x, (a) = < VD 3 (a) 


2 


x, = eonst. 


a.+2. [fear te. : 
2, (a)=x,+4-\le da, (), 


wobei die Funktion vx, («) nach Gl. (6) bekannt ist. Der Notwendigkeit, daß .r, (a) eine monotone 
Funktion darstellen muß, wird dadurch genügt, daß durch richtige Wahl des Vorzeichens 
von (a) die Ableitung f’(a) positiv gemacht wird. Die Funktion «#,(P) findet man aus 
B=Pß(x (a)). An die Genauigkeit der nach (7) erforderlichen Quadraturen, welche nach einem 
beliebigen Verfahren der angewandten Mathematik durchgeführt werden können, werden 
übrigens keine großen Anforderungen gestellt, da die Genauigkeit der Doppelleiter, wenn die 
Funktion x, (a) bzw. a(x) einmal festgelegt ist, lediglich von der Genauigkeit der Berechnung 
P(a(x)) bzw. &,(ß) und der Herstellung der Zeichnung abhängt. 

Als Beispiel sei die Beziehung zwischen Fallhöhe und Fallzeit eines Körpers im luft- 
leeren Raum, also die Gleichung %—=a? behandelt. Die Lösung liefert nach (6) und (7) 


2,(a)=r,+4- [k-a—log(i+k-a)] 


), = const. 


? ) 
mit k=2. ze. Für 
; p=h = Fallhöhe, 
+) 
a | <= t = Fallzeit, 
4 
KO m h 5000 m, 
”) 
nV = =in, 


Hd, = 2W m 


ist diese Funktion (also =) in Bild le für gleichmäßig verteilte Verzerrung dargestellt. 


Es ist hier 1 <0,0I6 = 69, 


VYx w= 
DR; 
(1490 t) 


R Rh 
H- 5000 
_F 4800 
4600 
H- 4400 
_F #200 
H 4000 
+ 3800 
H 3600 
_F 3400 
H 3200 
"F- 3000 
H 2800 
H- 2600 
H 2460 
H- 2200 
- 2000 
7800 
H- 7600 
7400 Bi q@ [A202] 

Fr 7200 Bild I (links). Bild 2. 


i i .t2 ie zie + de 
- 7000 Die Beziehung h = 9 + Die Beziehung « han 
s a unverzerrte ?-Teilung. 
r 800 m Mai , . 
a unverzerrte h-Teilung. bh unverzerrte «-Teilung. 
- 500 b unverzerrte t-Teilung. c gleichmäßig verteilte Ver 
Meter ce gleichmäßig verteilte Ver- zerrung. 
zerrung. 




















a [420.1] 
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Als weiteres Beispiel mag die Funktion = 5 dienen. Man findet Bild 2c mit 
a= 30,35, 40, .... 100 und 5= 10, 12, 14,... 34, wobei 
r [74 
2, (a)= %, +4 -arcig 50 
, , 10 . 
Vx, (a) — Ve.) = 3500 ° 0,10 = 10 L. 
ar 
a 
Genauigkeitsfragen. 


Es ist nun wichtig zu wissen, ob durch die Verteilung der Gesamtverzerrung auf beide 
Leitern die Ablesegenauigkeit ungünstig beeinflußt werden kann. Es wird dabei die Annahme 
zugrunde gelegt, daß bei der Ablesung von Zwischenwerten die Unterteilung der Intervalle dx 
auf dem Träger x linear vorgenommen wird. Die Funktion $=Pß(a) sei als Doppelleiter 
„=, (a), e=.x,(ß) dargestellt; zu einem Wert «, soll der zugehörige Funktionswert 9, = P(a,) 
mit Hilfe der Doppelleiter aufgesucht werden. Es seien a, und a, durch zwei benachbarte 





Skalenstriche auf dem Träger « dargestellt, so daß la, —a,=%#,. Ferner sei a,=1,= 4a, 
und damit 
v=atm:V, vem<=|), 
ER m’ d u ) 
2.0,)=2,la) + m: d%,:%,(a,) + 5 u A Rn h, m») V= iS 1) : 


Geschätzt wird nun aber beim Eingang in die Doppelleiter nicht der Wert «x, (a,), sondern 
der Wert 

x, (a)= x, (la) tm: [*, (a) — 2, (a,)] 
bzw. 


g2 


(a)=rla)tm N ala) tm.g ala, +t,- 9). 


Beim Eingang in die Doppelleiter begeht man also einen Fehler 1#* = .r; (a,) — x, (a,), der 
nach Durchführung des Grenzüberganges und in Angleichung an die Lage von a, im Intervall 
a, — a, durch 


2 


E23 Ur [23 > 
"= x (a): m-(l m) a RE So 


dargestellt werden kann. Bei der Ablesung des Funktionswertes 5 wird hierdurch allein eine 
Fehlablesung von (näherungsweise) 


Au* 
A =Pa)- Ar =>; Da ee aa 2. Sr 
x, (P) 
verursacht. Doch auch die Verzerrung der -Teilung auf dem Träger x gibt bei der im 
& f 8 8 
Intervall dJe=|x, (P,) — x, (P,)| linear vorgenommenen Ablesung Anlaß zu einem fehler- 


haften Ergebnis. Ganz entsprechend findet man für den durch lineare Unterteilung der 5-Skala 
verursachten Ablesefehler den Betrag 
g2 „! 
3 die ( 
ae Ah a I ı Be | ee | 3 
= %, (Bo 
Der durch linear vorgenommene Einschätzung 


samtfehler 14/ ist also nach (8), (9), (10): 


yei Eingang und Ablesung verursachte Ge- 


1 1 : 
I=5. m (Im) vn Pla) n (ln) Bun): - - - (MM). 


Zur Beurtei einer Beeinflussung des esefehlers Aß durch eine irgendwi ählte 
Zur Beurteilung einer Beeinfl g des Ablesefehl 15 durel rgendwie gewählt 
"erteilung der Gesamtverzerrung der gegebenen Funktion A(a) au ie beiden Teilungen 
Verteilung der ( t g d gel Funkt ( f die beiden Teilung 
x, (a) und #,(P) der Doppelleiter sei mit beliebigem Parameter p die allgemeine Beziehung 
angesetzt: 


2 ; Vxrz( By = P° Vx, (@) » 
Nach (4), (3), (11) ist dann 
v3 (a m: (1 — m) n- i—-n 7. 
Ap Fi . 9: la)+ ir 9-50 
1 u u »; pP (a) 
5 Ba. . e -(1— m) n-(1—ı y 
Ersetzt man zur statistischen Untersuchung die Größen ni = und ( 2 4 welche je- 


weils die Lage des betrachteten Argumentes im Intervall kennzeichnen, durch den (gemein- 
samen) Mittelwert 





YıınAa 


ch. 
941 


it 


de 
ne 
1x 
er 
ft) 
te 


Gy 


rn 


ler 


all 
5) 
ne 


9) 


im 
er- 
la 


je 





WVYııaa 
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) 
"m-(1-— m) 1 
\ - -dm= 5° 
0 4 F 
% 
so erhält man für den im Durehschnitt zu erwartenden Ablesefehler 
d 


9, L , 5 
| ber wlr 
"2 fi LY? 


| 
Z_ | 
er 


v3: da: d5- ? 
! (@) Ja ' ’, je 
l r P3,'P (a) 


Praktisch kommen nun lediglich die Fälle 


l \ ’ 
1. ? =(0 (a-Teilung linear), 
2. p—=0 (P-Teilung linear), 
3. p=1 (gleichmäßig verteilte Verzerrung) 
in Frage. Für die einzelnen Fälle findet man: 
ne 
1. 1p = BD "U3(a) * DE . ß (a)' 
rt 
2. IP= 75 va da Pla), 
u 1 
3. 4p — BR "3(e) * Dr U: 
+ a DE 


u? 3” a fi 0 

9, ( TE (a) 

Da man aus Gründen der Übersicht praktisch wohl stets die Beziehung #; = ß' (a) -#, für 
die mittleren Lagen des Gesamtbereiches anstreben dürfte, sind sich von der Seite der Ab- 
lesegenauigkeit her die beiden ersten Darstellungsarten gleichwertig, während im 


Falle gleichmäßig verteilter Verzerrung der mittlere Ablesefehler 4/5 auf das 0,7-fache sinkt. 


Im übrigen ist zu beachten, daß in den Fällen p>0 — die Ableitung ?’(a) kann ja nach 
einer früheren Bemerkung stets als positiv angesetzt werden — der Ablesefehler 4P stets 


einseitig liegt. Es kann daher zweckmäßig sein, bei der Anlage einer Doppelleiter z. B. mit 
p=1 die 5-Werte um den Betrag des auf der nichtlinearen Teilung der Skalen beruhenden 
mittleren Fehlers 45 zu korrigieren. Wird derart bei der Darstellung der Funktion ß(a) als 
Doppelleiter an Stelle von #,(ß) die Abbildungsfunktion 


r 9; rt a) : 
3 E \ :x (P) 
v2 P (a) 


benutzt, so sind bei der Ablesung sowohl positive als auch negative Fehler zu erwarten, deren 


. Dr 
Ko, (P) _ 12. 


. 0 DR &; .. 
Grenzen bei Ip und + -. liegen, deren Schwerpunkt jedoch bei Null liegt. 


Erweiterung. 

Bisher wurde die (unverzerrte) lineare Teilung als anzustrebendes Vorbild genannt. 
Häufiger jedoch wird in der Praxis die Aufgabe auftreten, daß eine Größe a z. B. an einem 
Meßinstrument abgelesen und eine von a abhängende und einer Doppelleiter zu entnehmende 
Größe b=b(a) an einer anderen Skala einzustellen ist, wobei im Allgemeinfall die Teilungen 
a und b an den betreffenden Instrumenten nicht linear aufgebracht sind. Diese an den 
(Geräten angebrachten Teilungen sollen vielmehr durch die monotonen Funktionen a=a (a) 
und b=b() dargestellt sein, wobei « und 5 die linearen Abmessungen der Teilungsträger 
bedeuten. Dann muß für die Doppelleiter, welche die Beziehung b=b(a) wiedergeben 
soll, natürlich angestrebt werden, ihre Teilungen a und b möglichst den Teilungen a und b 
der Instrumente anzugleichen, denn dadurch wird bei der praktischen Anwendung die 
Übertragung der einzelnen Zahlenwerte erleichtert und der Übertragungsfehler herabgesetzt. 
Dies bedeutet aber andererseits, daß man versuchen muß, die Hilfsgrößen a=a(a) und 
ß = P(b) möglichst unverzerrt wiederzugeben. Damit ist nun der Allgemeinfall auf den schon 
ausführlich behandelten Sonderfall zurückgeführt. Sind also an den Instrumenten die 
Teilungen a und 5b nach den Funktionen a=a (a) und b=b(ß) graviert, so hat man einfach 
die Umkehrfunktionen a=a(a) und 5=P(b)= (a) zu bilden und mit diesen Zahlenreihen 
die oben abgeleiteten Formeln anzuwenden. Die Bezifferung der Doppelleiter erfolgt natür- 
lich nach den Größen a und b. 240 
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Über den Zusammenhang einiger Formeln für den hydro- 
dynamischen Auftrieb. 
Von F. Jung in Wien. 


Es werden zwei Identitäten abgeleitet, mittels derer die Ausdrücke für den 
Auftrieb und das Auftriebsmoment in solche Gestalt gebracht werden 
können, daß sich daraus verschiedene bekannte Formeln auf einfache 
Weise ergeben. 


ine vollkommene (d. h. nicht zähe) Flüssigkeit übt auf einen allseitig von ihr umflossenen 

Körper nur Normalkräfte aus. Sie lassen sich, wie jedes räumliche Kraftsystem, zusammen- 
setzen zu einer Einzelkraft R®) und einem Kräftepaar, dessen Momentenvektor NP abhängig 
ist von der Lage, die man dem Angriffspunkt O0 der Einzelkraft gibt. Für die stationäre 
Strömung einer unzusammendrückbaren Flüssigkeit wurden die betreffenden Formeln entwickelt 
von Kutta-Joukowski, Blasius, Prandtl, Munk, Liebmann und Lagally'), meist 
getrennt für die räumliche und ebene Bewegung. Die nachfolgende Darstellung kommt zu 
dem Ergebnis, daß die erwähnten Ausdrücke sämtlich als Abkömmlinge zweier einfacher 
Formeln (14) und (18) aufgefaßt werden können, in denen statt der Flächenelemente df der 
Randflächen deren Spiegelbilder df* an den Geschwindigkeiten 7 auftreten. 

1. Die Ausgangsformeln. Es sei df ein Element der Oberfläche K des umflossenen Körpers, 
als Vektor aus der Flüssigkeit hinausgerichtet, r der Ortsvektor des Schwerpunktes von df, 
gezogen von einem Punkt 0 des starren Körpers, und p der spezifische Druck der Flüssigkeit, 
dann liefert die Zusammensetzung der an der Körperoberfläche X angreifenden Druckkräfte 
die Einzelkraft (hydrodynamischer Auftrieb) 

Ro=\pdf ee ea ee a 
K 
mit 0 als Angriffspunkt und ein Kräftepaar mit den Momentenvektor (hydrodynamisches 
Auftriebsmoment) 
No —=\rpxdf REN RIP NIT ZT GE 
K 
beide Integrale erstreckt über die Oberfläche A des starren Körpers. Die Flüssigkeit wird 
unzusammendrückbar und homogen angenommen und stationäre Strömung vorausgesetzt. Dann 
genügt der Druck längs einer Stromlinie der Bernoullischen Gleichung. Bei Abwesenheit 
von Massenkräften lautet sie 
Fin . 
a a N EEE ET 


wo v die Geschwindigkeit, « die Dichte bedeuten und die Konstante ce im allgemeinen ver- 
schieden ist für verschiedene Stromlinien. Die Gleichung sagt aus, daß längs einer Strom- 
linie die Gesamtenergie unveränderlich ist. Die Stromlinien gleicher c-Werte liegen auf 
Flächen, Energieniveauflächen, auf denen bekanntlich auch die Wirbellinien verlaufen. Durch 
Finsetzen aus (3) in (1) und (2) folgt 


RW —— u \ 5 df + \ edf, 
ei K 

Nr u \ rX 4 df-+ \ rxXxedf. 
K x K 


Angenommen werde, daß die Körperoberfläche AK Teil einer 
Energieniveaufläche ist, daß also die Strom- und Wirbel- 
linien auf K verlaufen (Bild 1). Dann tritt in den zweiten 
Summanden der Faktor ce vor das Integralzeichen. Die 
Integrale 





\df und \rxXdf Aaısa) 
F F 
: ER Bild 1. 
können nach dem verallgemeinerten Gaußschen Satz be- da 
rechnet werden. Ist @ eine beliebige Feldgröße, dV ein K K_ Körperoberfläche. 
Raumelement und bedeutet das Zeichen n eine beliebige Ri ee 


Multiplikationsart, so ist £ 


1) Vgl. die Angaben bei M. Lagally: Berechnung der Kräfte und Momente, die strömende Flüssigkeiten auf 
ihre Begrenzung ausüben, Z. angew. Math. Mech. Bd. 2 (1922), S. 409. 





XUM 


AUM 
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)- \dVr.-FaA@)=\dfn@ EL ee . 4), 
V F 


wo F die Randtläche des Raumes V ist. Für @ =1, bzw. —r folgt 


\df =\dV.WD 0, 
en F \ 
en = - y ,— 
ru \dfxr=\dV.-WVxXrn)=V0. 
F v 
Es wird also 22 in u n 
en RW—=uR, NP=uN, (5). 
.n- wo fü a. 
ig R=--\sdf, =-Irxadf. ee 
re ‚2 ; 2 
K K 
elt e * 
ist 2. Umformung der Integrale. Um zu ermitteln, wie R und \, von der räumlichen Ver- 
zu teilung der Quellen und Wirbel 
nr ee „*, > VA 
er 


abhängen ?), sind die beiden Flächenintegrale umzuformen in Raumintegrale, in denen q und w 
"Ss auftreten. Zunächst ist nach (4) 


f. r u { ad i v? i v 
it, \> df \a Erz, IX: \a vrx|r z1- 
te BE r io F a f 


Da aber 





es 
und der Rotor von r verschwindet, so ist 
2) _ Pin “ " 
\aixr,=lav.[rg)xr. 
rd F j 
er Aus den bekannten Identitäten 
i 
rd) +leilır, 
3), Wr. ne er 
"5 =eiIr+rXUxXe 
r- 
m- folgt durch Elimination von ("|J )r die Beziehung 
uf „ 
-h er - VZzerWxXnN)XxXr+HU | 
Pr (9). 
=wXxr+gq " | 
aus der man 7 — entnehmen und in die vorhergehenden Raumintegrale einsetzen kann’). 
Als bequemer erweist es sich jedoch, einen etwas andern Weg einzuschlagen. 
Mit Benützung des Einheitstensors ? läßt sich die linke Seite in (8) schreiben 
+ — = v° - , pe > we. « 
FEBSFZETT II SF ZEN ne 
2 ui ui / 
sie erscheint somit als Feldableitung eines symmetrischen Affinors oder Tensors. Der Kürze 
=. 
2 Re. . ’ 
halber sei dieser mit 5 bezeichnet, allgemein also 
ee 5. Se Rue u 5° 
Die Bedeutung des Tensors &* ergibt sich, wenn durch Multiplikation mit ihm einem beliebigen 
Vektor @ der Vektor ä* zugeordnet wird: 
2) Es bedeutet #/5b das Innenprodukt, #Xb das Querprodukt (Vektorprodukt), 7,5 das Allgemeinprodukt (dya 
disches Produkt). 
uf 3) Vgl. die obengenannte Abhandlung von Lagally. 








a2 1a 0 
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«=äae*=äll2e,8—2) | 

=2(alye—[alde (le xXaXN!. . :. -. -».:.::2.. 
—= (a eB)E+ÄÜxXE)XE | 


wogegen 

DE — (ÄXE)XE,. 

Der Vektor @* entsteht also aus @ durch Vorzeichenwechsel des zweiten Summanden. Der 
erste stellt die Komponente von @ dar in der Richtung von 2, der zweite die Komponente 
quer dazu. Diese hat bei &* das entgegengesetzte Vorzeichen, @* ist also das Spiegelbild 
von ä bezüglich der Richtung #. Man könnte den Tensor &* Spiegeltensor*) nennen, # seine 
Achse, und sagen: Innenmultiplikation eines Vektors @ mit einem Spiegeltensor bedeutet 
Spiegelung von @ä an der Achse des Spiegeltensors. Stehen also Vektor und Tensorachse 
aufeinander senkrecht, so kehrt der Vektor seinen Sinn um. Legt man # in die Richtung 
eines Vektors ü, so ist wegen (10) offenbar 


ä=äe=(i 


LEERE RETTET: 
und wegen (9) 


und die Identität (S) nimmt die Form an 
=% 


| =gd+mXd. En DE RE Ir 1 Ze 


—_ 


Die Integration über den Raum gibt für die linke Seite nach dem Satz von Gauß mit Rück- 

sicht auf (11) und (12) 
\dvr.r#=\dfl=|df!*v* 
r F F 


\o’dj* 
R- 
und daher 
a ET REN 
\>d/* —\dl 0 te ME Er BE Er | ° 
F v 
Fällt # mit einer Energieniveaufläche zusammen, z. B. mit der Körperoberfläche K, so stehen 


dl 


die d f senkrecht auf den r, den Achsen der zugehörigen Spiegeltensoren 7*, und dann ist hiernach 
dem oben Gesagten 


df*=-—df 
und das Flächenintegral geht über in den Ausdruck (6) für R, so daß 
— Ren 
R N ° 
K 


Um Momentengleichungen zu erhalten, die den Gl. (13), (14), (15) entsprechen, werde die 
Bet : Ya ee er 9: Ru... i j 
G1.(13) quer (vektoriell) mit r multipliziertund dann die linke Seite [ 5x? so umgeformt, daß 
man gerade eine Feldableitung erhält. Nun gilt allgemein für einen symmetrischen Affinor x 
die Beziehung 


EEE SE ER N T F 
Um sie zu begründen, setze man . 
En Br et ar’ j 
=, (m,W% +Wi ,U), IE: 


i 
wo n mindestens gleich drei ist, eine Form, in die sich jeder Tensor bringen läßt. Die 
Jeiger i sollen im folgenden der Kürze halber weggelassen werden. Zunächst wird 





rı(NW,aN)xr=r| YWWxn)—= ViriWtoa)m)a’xr) 


> 1a)(a” x) + ta |) RNKXTr +" X tw |V)r}). 
Da aber 
(If)r=WV,r=We=W, 
so ist 


(dw, a”)xr= ri +a|n)at)xXr+a’xXW. 


4%, Bei Gibbs heißt er biquadrantal versor, da man @* auch als Umklappung von @ um © ansehen kann. 
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Vertauscht man # mit #” und addiert zur vorigen Gleichung, so folgt 


rıEN@,a +0, xXr=| dr, + xXr, 
also die Gl. (16). Mit Rücksicht auf sie lautet jetzt die quer mit r multiplizierte Gl. (13) 


R, a 2 
J (> xrl=lyE+RXEXKT . Er > RT ae 


Wie bei (13) ergibt sich 


dv. vr @Xn=\dferXF)=\(dF RX F 
4 r 


) F 
=\M"df*xr 
; 
und daher wegen (17) 
v’ 7 _ ; r _ — .n . ) 
\>df"xr- \dav.(ge+wxXerxXr Se Bra ee Re 
F v 
Für eine Energieniveaufläche, z. B. K, wird wieder 
j _ n” z% ( 
N=\rx>df*.... ud 2 ... (19) 


K 
3. Anwendung auf R und N”. Die Raumintegrale (14) und (18) seien nun erstreckt über 
den Raum A zwischen der Körperoberfläche K und einer sehr großen Kugelfläche x, die man 
dann ins Unendliche rücken läßt. Dann wird 
z f %a 7% { r -— 1 —\,->- 
\>d/* +\„ıaf -\dV-(ge+m) v), 


K oo A 


\rx —df* HrxX —df* = \r xXdV-(grtmxXe), 


K 2 A 


a » 


R=-\>d/* +\dV-(yr+wxXe). | 
- i | 
} E n ( (20). 
RR ä 


Mit « multipliziert entspricht die erste Formel dem Ausdruck von Munk-Lagallv, die 
zweite dem von Liebmann-Lagally°). 


Es bedeute 


R,: \ „daft, R,=\dV-(grtwxe), R=R+R .. (eu, 
Be A 
N,=-\rxgaf, Iu=\rxav-getaxd,N=N,+N,.. (1, 
"ae 4 
also 5 En N 34 es Me 
RW) = u(R,+R,), rm, N 0.0.0.2 


Die Ausdrücke für R, und N,, können wieder umgestaltet werden mittels der Formeln (14) 
und (18). Nach ihnen ist 


-R,=\dV-.(ge +wxXe), N,=\rxdV-ge+wxXD), 
J+A J+4 
wo J den Raum innerhalb der Körperoberfläche bedeutet, die Raumintegrale also über den 
ganzen von der Fläche © umschlossenen Raum zu erstrecken sind. Um dies zu ermöglichen 
muß die Strömung in das Innere des Körpers fortgesetzt gedacht werden. 
Im Unendlichen sei die Geschwindigkeit *, überall der Größe und Richtung nach gleich. 
An einer beliebigen Stelle im Endlichen sei sie zusammengesetzt aus v, und einer zweiten 
Komponente »’, 
EEE En en 


5) Vgl. a.a.0,. die Formeln (7). 








AUM 
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Daher ist 
—R=\dV.g@, +) +RXW +), 
I+A4 
\dV-(ge,twxXF,)+\dV .(gvV’ +wxXev) 
J+4 343 
und ebenso 
-N,,=I\rxdV(ge,+twXp,)+\rXdV- (gr +mxXr’). 


J+4 J+4 


vyle,=6, FX,,=®, 


Nun ist 


also nach (7) 6 ? 4 
q=V/|e@,„+Pv) =r/|v" =Yd, 


n„=yxXtw,tV)=/ xvV’=m, 
und folglich 
R,=\dV.(ge,tWxXe,) +dV- ri m'xe), 
+44 J+A 


N.=\rxdV- (ge, twxXe,)+\rxXxdV- (de mxXr). 

I+A4 J+A 
Verwandelt man die Integrale mit den gestrichenen Gliedern nach (14) und (18) wieder in 
Flächenintegrale, so ist 


P ’ 


f e- ul * v 2 > 
\d Ve" +m'xXer) =\ u, 
Fi Fr 


3 A £ 


\Fxar tw) \rX af", 


J+4 o 


wo df’* durch Spiegelung von df an ®’ hervorgeht. Da »’ im Unendlichen verschwindet 
(wie r?), werden die Integrale rechts null, die zweiten Summanden in R, und N,, fallen 
weg und es ist 

R,=\dV-(ye,+wxXr,), | 





J+A 
R 2 (24) 
N,=\rxdV-.(gr„+twxXr,) 
5 J+4 
und daher®) nach (21.1), (21.2) 
R=—\dV.(ge,+wXr,) +\dV-(ge+mxn), | 
J+4 A 
er : (25). 
N,=—-IFXdVv.(gö.+wX5,) + FXAV-(g5+mXD) | 
+ A 
Die Integrale 
\dV.q=%, ae BEE 
V v 


lassen sich als Gesamtquellung, bzw. als Wirbelungsresultierende im Raum V auffassen und 
es ist 


R, —=QAyr4 v.+ W,+ 1X, . . . . . . . . . . (27). 
Da durch die Fläche K keine Flüssigkeit strömt, ist 
09,=0. 


Im Ausdruck für W kann man die Raumelemente so abgrenzen, daß sie mit den Elementen 
der Wirbelfäden zusammenfallen. Ist d! das Bogenelement einer Wirbellinie, d F ein Quer- 
schnitt des Wirbelfadens, so ist dann 

dV=dFjdi, 
und az ie A 

W=\(dFidhDw=\dF|wdl, 

weil w und d/ gleichgerichtet sind. Die Wirbelstärke 

dl’=dFw 


6) Vgl. a.a.0. die Formeln (24). 
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ist konstant längs jedes Stromfadens. Die Integration über den Raum V kann also so aus- 
geführt werden, daß man zunächst das Integral nimmt über jeden einzelnen Wirbelfaden, 
wobei das konstante d /’ heraustritt, und dann über alle Wirbelfäden integriert, 


ae 2 Me 34 © 


Bei der räumlichen Bewegung wird angenommen, daß alle Wirbellinien im Endlichen ver- 
laufen, also geschlossen sind, weshalb 
\ dal=0, W,+4 —D: 

und daher 1 

ET 
Bei der ebenen Bewegung stehen die Wirbellinien senkrecht zur Bewegungsebene und können 
deshalb nicht geschlossen sein. Für eine Flüssigkeitsschicht von der Dieke der Längeneinheit 
(# senkrecht zur Bewegungsebene) wird 


\di=e, W,+4=e\dI!'=el',,,, 
folglich 337 
W ET U MU 2 MPU: > 4 SCHE e 1 
Setzt man im ganzen Außenraum A die Bewegung als quell- und wirbelfrei voraus, so ver- 
schwinden die nur auf ihn bezüglichen Größen 
Q,=0, R,=0,;, N.,=0, 
und es wird für die räumliche Bewegung wegen (29) 
ee ee. ME 
Für die ebene Bewegung wird /',—0, wegen (30) folgt also 
Ro u da an a nr 5 a 
die Kutta-Joukowskische Formel. 
2. 4. Darstellung durch komplexe Größen. Wenn man bei der ebenen Bewegung R, und 


N,, in komplexe Form bringen will, so ist zu setzen 
v=vef, df=df-e«, 


also, wie aus Bild 2 zu sehen 


df=df-t?-o, 
und daher 








r . dz __ 
"df*=(df-e-ie)(v ei’). af 
Deutet man eine komplexe Zahl durch Unterstreichen M Ar K 
an, z. B. 2, ihre Konjugierte durch z’, so lautet die Azıs.2] A2is.3 
vorige Beziehung Bild 2. Bild 3. 
v”df*—=v ee a a en I ER Fi 
Nach Bild 3 ist 
df=-idz, 
die konjugierte Größe folgt durch Vertauschung von i mit —i 
BEE 5 88 5 I, 
und es ist wegen (33), (34) 
a w 
Re\gdfsilsedr.....2222..:8 
Der Kürze halber sei 
.2 1 
= BEITRETEN, nen Ei sa 


gesetzt. Der in N,, auftretende Produktvektor rX dü steht auf der Bewegungsebene senk- 
recht, unterscheidet sich also für verschiedene Stellen nur dem Werte nach. Setzt man 


r=2=rer dü=du=du:ei!, 
so 1st 2 
rxXdü=rdu-sin(y—g_), 
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und dies, mit ö multipliziert, kann man unter Beachtung von (36) auffassen als Imaginärteil 


der komplexen Größe 


rdu-[cos (y — Q)+isin (y 


y))=rdu-iW-N—(re-iv)(duei”) 





oder als Realteil der aus ihr durch Multiplikation mit — i folgenden 


rdu- [sin (y — 9) 


Demnach wird 


2 


N,. -\#x 5 


’ 1 
-icsy— p)=— 


af=\rxau=rl vz’de. 


sh, 54 





(38). 


Nun hat funktionentheoretisch nicht v, sondern*die konjugierte Geschwindigkeit v’ eine ein- 
fache Bedeutung. Man erhält sie in den Ausdrücken, indem man statt R, und N,, die kon- 


jugierten Größen nimmt Ri, N,.- 
wirbelfrei ist, so hat man 
1 


5tru 


\vidz, 


Ri — 


x 


a 5: £ 
NPI=-—-R— u \v’za z 


Wird noch vorausgesetzt, daß die Strömung quell- und 


(39), 


a 


weil die Realteile konjugierter Größen gleich sind. Die erhaltenen Ausdrücke sind die Formeln 


von Blasius. 
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Beitrag zur geschlossenen Näherungs- 
darstellung elliptischer Integrale. 
(Mitteilung aus dem Institut für elektrische Schwin- 
gungslehre und Hochfrequenztechnik an der T.H. 
Berlin.) 

Für die elliptischen Integrale bestehen außer den 


tabellarischen und zeichnerischen Darstellungen 
der Funktionentafeln (z. B. von Jahnke-Emde) 
Näherungsdarstellungen in Form von Reihen- 
entwicklungen, die für kleine Bereiche des 
Moduls k Gültigkeit besitzen. 

Bei manchen technischen Anwendungen ist es 


erwünscht, eine Näherungsdarstellung zu kennen, 
die eine geschlossene Form besitzt und eine mög- 
lichst einfache Funktion ist, so daß sich z. B. 
weitere Integrationen oder Differentiationen nach 
dem Modul k durchführen lassen. 

Im folgenden wird gezeigt, wie man Näherungs- 
funktionen mit sehr kleinem Fehler für vollständige 
und allgemeine. elliptische Integralei. und 2. Gat- 
tung gewinnen kann. 


1. Das vollständige elliptische Integral 
1. Gattung (K). 
a) Einfache Näherung im Bereich 0<k<0,9. 
Es ist 
2 
K (k) = \ Yi 


v0 


dgq 


A? sin? g 


Nach dem Mittelwertsatz der Integralrechnung läßt 
sich eine primitive Näherung finden, wenn man die 
Wurzel mit einem Mittelwert von € vor das Inte- 
eral zieht: 


1 
K(k) = dg= 
yıi-k'sinpn. 
0 


ıj2 


Vi-—k’m 


Hierin ist m — sin’, ein zunächst unbekannter 
Wert zwischen O0 und 1, der sich aus den Tabellen- 
werten von K für verschiedene k bestimmen läßt. 


. (la). 


Die Rechnung ergibt folgende Tabelle: 





@ 0 10° | 20° 40° 60°, 80°| 85° 89° 90° 


0 0,1740,3420,643 0.866 0,985 0.996. 0,9998 | 1,0 


m 0,500 0,5020,5170,5520,627.0,7750,840 0.920 1,0 








Für größere Werte von k ist m nicht konstant, 
sondern verändert sich mit k. Nimmt k nur Werte 
von 0...0,3 an, so isst m=0,5 und 

7 1 

2 y1—05R? 
Für Werte von k zwischen 0. . 
Mittelwert mn — 0,6 einsetzen. 


K(k) (2a). 
.0,9 kann man als 
Dann ist 


T 


K(k)= 2 für k = (3). 


2 y1- 06%? 
Der Fehler beträgt etwa1% 
3ei k > 0,9 versagt (3). 


0,9 


im Bereichk = 0...0,9. 


Beim Differenzieren von K entsprechend (3) 
würde die Näherung m = 0,6 für “x bzw eK USW 
” dk dA? 
zu erheblichen Fehlern führen. 


Es ist besser, m in der Form m = 0,5 + 0,15 k? 
darzustellen. Damit wird eine gute Näherung 





1 


K(k) 
< v1—0,5%: 


Br (3a). 











Der Fehler liegt unter + % %, wenn a den Wert 
60° nicht übersteigt (kS 0,87). 

b) Gute Näherung für beliebiges k durch 
schrittweise Integration mit Fehlerkompensation. 

Eine bessere Näherung als (3) erhält man durch 
folgende Überlegung: Die unter der Wurzel stehende 
Funktion v— 1— k?sin?g ist eine zup = n/4 sym- 
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metrische Funktion. Diese stetige Funktion y kann 
man durch zwei Näherungsparabeln yı und y» er- 
setzen (Bild 1). Yı wird im Bereich 0.../4 und 





—.} 
DZ 

l a ER 

1% W-> L74 


v 


Bild 1. Näherungsdarstellung von 1-- k?sin?Yy. 
ya zwischen x/4 und x/2 benutzt. Die.Fehler werden 
sich, da bei =/4 die Abweichung der Kurven yı und 
Ys von der wahren Kurve y gleich groß ist, nahezu 
kompensieren. 
a) Integration von %, fürg =0...n/4- y,=1-— k?’g? 
a? 
u dg 1 i 
R= \ —— = aresin (krıld). 
2 ] 1 k? U Ai k 
1) 
3) Integration von %, für —=n/4...n/2- 96 
; n\° 


a2 


dg 


K"” 7 2 
| | ö a. 7 k? (4 = ) 
7j4 4 


] 


” i In (rs 7 )+ V5-+k(p -a/2r)) 


1 vs 

1; In s 

; Vd+(kals’— kald 
Damit ist der vollständige Ausdruck für KK’ 
TR 
K(k) = | 


| 6) (4). 


k aresin (k a/4) + In 
Vötikafs®— kaja 


mtö=1—-KR—XR”? 





Dieser Ausdruck ist die Näherungslösung, die für 
beliebige k zwischen 0 und 1 Gültigkeit hat und 
um etwa 1% zu kleine Werte gegenüber den 
Tabellenwerten (Jahnke-Emde, 1. Aufl., S. 55, 57, 
59) ergibt. 

Für numerische Rechnungen ist (4) nur im Bereich 
k=0...09 bzw. a=0...70° bzw. 5=02...1 
geeignet, da für k>1 die Ausreehnung des In 
ungenau wird. In diesem Fall entwickeln wir 
den Nenner 








Br ER FR 
Vd+(kapae kalız pl rn 
In se nt2_ 1-(V? \ 
Yd+(kn/s)’—kajd y® \kn/2 
ka/2 1 
ch + ka 


yö r N 


da | h) < kisr/2. 


Damit wird im Bereich 
k— 0,9...1 bzw. a= 70°,...90° bzw. d6—=0...0,2 
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K(k) = | 

Bi 2, (ka14) ka 1 ( (4a) 
7; |are sin (ka 4) —+In nz "er | 
y® | 


oder mit k=sina: 6 = cos?a; yö cosa=—k’ 


r 1 
K(e)= — 
| sin a ’ (4b), 
are sin (n/4 sina)+1In(tza-/2) 4 | 
| ba EIERN "tea-a/2)?) 


Wie man aus (4b) erkennt, wird 

K mit k>1 (a>nx/2) logarithmisch x. 
Andererseits ist der Grenzwert von 

K mit k>0 nach (4) 772, 
wie man aus den Reihenentwicklungen von (4) er- 
kennt: 


aresin(kn/4)=kafa+nk?+... 


In(Vö + (ka/s®) — kal4 
In(!—-kn/a+zk°’+...) 
kalA—zk? 


Der Fehler von 1%, den die Lösung (4) aufweist, 
läßt sich um mindestens eine Zehnerpotenz ver- 
ringern, wenn die ziemlich grobe Näherung nach 
Bild 1, die an der Sprungstelle immerhin eine Ab- 
weichung von +25 % bedeutet, durch eine bessere 
Näherung nach Bild 2 ersetzt wird. 


c) Näherung für beliebiges k mit Fehler < 1. 

In Bild 2 wird die Funktion 1—k?’sin’g in drei 
Teilabschnitten angenähert. Außer den beiden 
Näherungsparabeln verwenden wir in dem Bereich 
von 27,8° ...62,2° eine gerade Linie als Nähe- 
rungsfunktion %,, die Y bei 45° tangiert. Die 
Sprungstellen sind so gewählt, daß die Ab- 
weichung der Näherungswerte vom genauen 
Funktionswert gleich groß ist. Es liegt z. B. 
bei 278° die Gerade v; um 0,0177 über u, 
während vyı um den Betrag 0,0178 unter y bleibt. 
Die Abweichung von 0.0177 bedeutet bei 27,8° 
einen Betrag von 2,3%, bei 62,2° den Betrag von 
8,1%, während die erste Näherung im Abschnitt b) 
eine Abweichung von 25% ergab. 














T 

| 

k? 

| 

dk? 

nen ZZ 
99 2 450 20 30° 
Az 4 y—  Nf 


Bild 2. Gute Näherung der Funktion 1 — k?sin?y. 


a) Integration von Yı im Bereich $=0...0,485 
(27,89). 
0,485 
TRR | d« 1 ü 
K' \ . = are sin (0,485 A). 
; ) 1— Ko? k 


0 
f) Integration von Y» im Bereich $ = 1,086... 1,5708 
62,2°...909% 
1,5708 
o ii 
ver b 


\yö+R@— a9? 
1,086 
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ö 
= - In } 


k  Y8-+0,485 k)?— 0.485 k 
=1—k—=Kh". 


y) Integration von Ys im Bereich $ = 0,485... 1.086 





21.5°...085° 
1,086 1,086 
" dg i dgq 
K'” = _ 
\ 1 — k?sin?g \; 1—05%A?’(1—cos2Yp) 
0,435 0,485 
und mit =n/4—+ y 
v—03 
\ dy 
IN1-05Rr— Rey 
v— — 0,3 
” 2 
K”—= 3(V1-02%— Y1—08%)- 
Damit ist die Lösung K=K’+K”-+-K'” 
i 1 ® a 
K(k)== n are sin (0.485 k) 
+2) 1— 0.2 — y1— 0,8%? 
k mr: 
ö 
—+In 
Vd-+ (0,485 A)? — 0,485 k 
=] k?— k’? 
Diese Lösung besitzt einen Fehler unter 1/oo. 
Für numerische Rechnungen im Bereich k 


0.9.. 1 kann man analog (4a) und (4b) entwickeln 








k—09...1 bzw. a — 70°...90° bzw. 6 —0...02 
E 1 h w 
Kik) k arc sin (0,485 A) 
1 - 09 x 2 
+>V1 0,2% y1- 08% ai 
k 
VITA, 1 
+97, + osTA® 
(y3) 
K (a) ; arc sin (0,485 sin a) 
sin a 
+2/1-92sina — Y1-0,8sin?a | p) 
sin a 
1 


+-1n(0,97tga)—+ 





(0,97 ga)? 


2. Das vollständige elliptische Integral 
2. Gattung (E). 
Es ist 
z]? 


E(k)=\dpyYi-#ksin”e. . . . (6). 
0 


a) Gute Näherung auf Grund des Mittelwertsatzes. 

Der Mittelwertsatz führt im Gegensatz zum Inte- 
gral 1. Gattung (K) hier zu einer sehr guten An- 
näherung, weil E nur zwischen =/2 und 1 schwankt, 
wenn k von 0...1 geht, während K im gleichen 
Bereich die Werte von x/2...oo durchläuft. Es wird 


Vi-k:m\dy- 


E(k)= aj2y1—k?m (7). 


Für kleine Werte von k errechnet man m zu 0,50. 
während m für k=1 den Wert 0,595 annimmt. 
Demnach ist m nur geringfügig von k? abhängige. 
Eine sehr gute Näherung, die für k -0...1 einen 
Fehler unter % % liefert, stellt die Funktion dar: 














1 
E(k)=Yy1-— %?(0,5 + 0,095 AP) 5 
= (8). 
a 
E(k)—= Y1- 0,5 k?— 0,095 A175 
b) Näherung mit Fehler < 1 "/o. 
Die Methode der schrittweisen Integration mit 


Fehlerkompensation läßt sich leicht auf die Be- 
rechnung von E anwenden und führt zu folgendem 
Ergebnis: 


E (k) — 0.2425 (5 + (0,485 A)? 
+ Y1-— (0,485 k)?) 
‚Vi 022" — YV1--08%:" 
1.5 A? 


(9). 


+; n arc sin (0,485 A) 


1) 
] Ö u: (0,485 k)? _ 0.485 k 


—+öln 





Die Entwicklung für Zahlenrechnungen im Bereich 


k—09...1 bzw. « = 70°...90° bzw. ö6=0...02 
ist 
E (k) = 0,2425 (y5+ (0,485 k)? 
+ Y1-— (0,485 k)?) 
.V1-02#’"—y1 0,8%" 
1,5 R? 
(9a) 


1 
+ >% arc sin (0,485 k) 


0,97 k ö 
+ ln + 097% 


er 7) 


E (a) = 0,2425 (Y 5 -+ (0.485 sin a)? 
+yı1 





oder 


(0,485 sin a)?) 


8 3 
= V1-02sin’a — YVI1—0,Rsin?a 
1,5 sin? a : (9b), 
- are sin (0,485 sin «) 
Y 2sina N 


cos? a 


32 97 
—+ cos? a In (0,97 tga) + (0,9 tga)? 





Meist wird man die hohe Genauigkeit der Formeln 
(9), (9a), (9b) nicht benötigen. sondern auf den 
einfachen Ausdruck (8) zurückgreifen. 
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3. Das allgemeine Integral 2. Gattung. 
Die angegebenen Näherungsmethoden bleiben 
keineswegs auf die vollständigen Integrale be 
schränkt, sondern führen zur Darstellung der all- 
eemeinen Integrale 1. bzw. 2. Gattung, wenn man 
sinngemäß statt =/2 als obere Grenze 9 einsetzt. 
Denn es ist definiert 


4 
E dgq 
F(k,Q)= EEE." 
‚ Prag er 
u 
4 
Ek,y)=\deyi-Rsing - - . (MM. 
0 


\uch «lie Lösung mit Hilfe des Mittelwertsatzes 
1äßt eine geschlossene Darstellung des allgemeinen 
Integrals zu, wie für E(k,g) gezeigt werden soll. 


Berechnung von E(k,_). 
Der Mittelwertsatz führte für g a/2 zur Lösung (8) 


Vi — k? (0,5 + 0.095 k®) = . 


E (k, ııl2) 
Es ist nach (11) 
4 
E(k,y)=)\dgy1-—#sin?g. 
u 


Zunächst seien die Grenzwerte für k— 0 und 1 be- 


rechnet: 


a)k—0: 
E (0.9) q ee Er Re 9 EEE 
b) k=1: 
4 
E(l,g)=\deYV1-sin?g sing. . (BB). 
0 
Es ist also 
un 06 bei ELad)=1.. (1. 
p=nl2 


dg 


Wir setzen daher den Mittelwertsatz im Anschluß 
an (8) in der Form an: 


Y 
E (kg) = Vi -- A? (0,5 + 0,095 AM) /(g)\da 
u 


E(k,Q)=Py1— RR? (0,5+0,095 A) (p) - (DB). 
Dieser Ansatz erfüllt für k — 0 selbsttätig (die Be- 
dingung (12). Die Funktion fiy) muß man so be- 
stimmen. daß bei k—1 auch E (k,p) den Vor- 
schriften (13) und (14) genügt. 

Es ist 


E (1.9)= p y1—0,595 /(9) (15a). 


Man kann vermuten, daß fig) nur Potenzen von 9? 
enthält, weil in der Definitionsgleichung von E nur 
sin’ auftritt. Wir versuchen, die zwei Bedingungen 
(14) mit dem Ansatz zu erfüllen: 


(p) (5) +3) SE 

Für 9 = r/2 muß f(p) den Wert 1 annehmen: 
le 5 RE (>; 
Eine zweite Beziehung für «a und f erhalten wir aus 


der Bedingung (14), daB = 0 sein soll für = a2: 


lg 


dE . 7-0,595 f (g) 
= 11—05%/(p) — 


de 2 V1—059 (4) 
r e) i u} /20595 P’ (p) 
do/e=an Y1— 0.595 ayı 0.595 


0,637 — 0,733 f’ (p). 


Also hat /’(gp)xj2 den Wert 0,869 anzunehmen, damit 
1E 


(© —0 wird. Aus (16) folgt 
d Y/y=al2 
, 2 a w. 
P(p)ar = Zu E= a 0,869. 
Nach (17) sta+3=1 
also 3 0,317 
und a= 1.317. 


Damit foigt 
2 
3917| 9 917[ 9 
/(y) 1,317 (75) — 0,317, 


”) 
Eine bessere Annäherung für mittlere Werte von 
g liefert die Funktion 


( 2 
Ko) 134 (275 


0,543 2? — 0,056 q* 


u | 
084775) ü .. (18). 


Damit wird die Lösung 





E (k,y) 
7 l 1 - k? (0.5 + 0,095 k®) [0.543 p? 


E(k,g): = 


0,056 4° 


2 4 
oh/ 2 - znsilie v © 4 ] 
l 1 — A? (0,5 + 0,095 A |134(, 5) 0345) 











f (19). 
Diese Lösung liefert für beliebige € zwischen 
0...2/2 und k=0...1 Fehler unter 0,5 %. 

Für k=1 stellt E (1.9) gleichzeitig sing dar. 

Tatsächlich ergibt die Funktion 

E(1,r)=gY1-— 0,3233 1? + 0.0333 p* 
gegenüber sing zwischen 0 und =/2 Fehler von 
höchstens 2 %/oo. 

Im Anschluß an die Näherungsdarstellungen 
elliptischer Integrale seien als einfache Anwen- 
dungsbeispiele die Schwingungsdauer eines Pen- 
dels bei großen Ausschlägen, ferner der Umfang 
von Ellipsen beliebiger Exzentrizität in geschlosse- 
ner Form berechnet. 

Schwingungsdauer eines Pendels. 

Läßt man ein Pendel zwischen den Grenzlagen 
— a und +.a nach Bild 3 ausschwingen, so ist die 
Dauer der vollen Schwingung 
l \ dg 
2g. V eos [1 

u 


T-4| (20). 


cosa 


Bild 3. Bezeichnungen 
zur Schwingungsdauer 
des ebenen Pendels. 
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Nun führt U=2nay1-0,58. 0,095 e° 

a a2 oder 

4 dg f dy 

\ —#> \ RR a? + b? „_„[a? — b?\5 (23). 

| Veos g — cosa 2 ec =2a 0.095 a? | 
od vo | 1 — sin’ — sin?’y 2 \ ea 

—y2.K(k) Fehler < 0,5%, 


auf X, das vollständige elliptische Integral erster 
a F FR... u . 
Gattung mit dem Modul k=sinZ - Für Schwin- 
gungsweiten a des Pendels bis etwa 120° kann 
man näherungsweise mit einem Fehler von weni- 
gen °/o nach (3a) K durch 





2 y1— (0,5 0,15 A?) A? 
PR. 1 
2 11-—-05®— 0,15% 
darstellen. Damit wird die Dauer der Halbschwin- 
gung Tı,=2 | ! K(k) 
9 
1 
Fun | / 1 O5 ke 0.15 kr’ 
g yv1—05 ‚15 M 
k=sin s ‚as120° 








bzw. die Frequenz der vollen Pendelschwingungen 


Tel 
2% | I 
9 


a< 120° 





ORT er 
1-0,5sin?, —0,15sin! 2 


22). 











Für sehr große Schwingungsweiten (a > 120°), die 
allerdings selten vorliegen werden, kann man 
ebenfalls die geschlossenen Näherungsformeln (4) 
bzw. (4b) mit großer Genauigkeit heranziehen. 


Umfang der Ellipse. 


Bezeichnet man 


| a? — b? 
E= 
a? 


als die Exzentrizität einer Ellipse mit den Halb- 
achsen a und b, so ist der Umfang U der Ellipse 
eegeben durch !) 


U=4aE (n]2,e) 


2 
E (aj2,e)=\Y1--esin’gpdg 
0 


ist das „vollständige elliptische Integral zweiter 


Gattung“ mit dem Modul e. Nach (8) läßt sich E 
durch die gute Näherung 


A n 
E(a2,)= 5 V1 0,58 — 0,095 &1° 
mit einem Fehler <0,5% darstellen. Damit wird 
der Umfang von Ellipsen mit der Exzentrizität & 
bzw. den Halbachsen a und b 


ı), Hütte I. 26. Aufl. S. 119. 


Dieses Ergebnis kann leicht gedeutet werden: 

1.Große Exzentrizität: für e>1(bla>0) 
ist U=27ay0405=4a. Die Ellipse ent- 
artet zum Doppelstrich mit der Länge 2a. 

2.Kleine Exzentrizität: für e<08 





(b>06a) ist U=2n | "Te (23a). 


> 





nit einem Fehler unter 7°/oo gegenüber (23). 














a f) > R 

3 y% 
War 
3 + u > 

11. | 





Bild 4. 

Die Diagonale des der Ellipse umschriebenen 
Rechtecks ist nach (23a) gerade y2mal so groß 
wie der Durchmesser des Kreises, .der gleichen 
Umfang wie die Ellipse hat (Bild 4). Die Ellipse 
schneidet nun in Bild 4 die Diagonalen des un- 
schriebenen Rechtecks an Punkten P, deren Ko- 
ordinaten die Länge «=a/Y2 und y=b/Y2 
haben?). Daher gewinnt man entsprechend Bild 4 


a?+b? 
> 


Ellipse und umfangsgleicher Kreis. 


in der Strecke OP mit der Länge 2 


unmittelbar den Radius des Kreises von gleichem 
Umfang, wie ihn die Ellipse besitzt. 
Berlin. 0. Zinke. 220 
2) Über weitere bemerkenswerte Eigenschaften von 
Pı...Pa s. O. Zinke, Umfang von Ellipsen und Eigen- 
schaften der Ellipse an den Diagonalenschnittpunkten, 
Z. math. u. naturw. Unterr. 1941. 


Über Kurven konstanten Bahndruckes. 

1. Einleitung. Wir betrachten einen Massen- 
punkt, der in einem geneigten Rohr, siehe Bild 1, 
reibungslos hinabgleitet. bei B mit der Anfangs- 
geschwindigkeit v, in ein gebogenes Rohr eintritt 
und an diesem entlang der Rohrwand reibungslos 
hochgleitet. 





Bild 1. Bahndruck, am Massenpunkt und au der Rohrwand 
angreifend. 
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Der Massenpunkt verliert dabei an Geschwindig- 
keit, und zwar folgt diese in einem beliebigen 
Punkt P zu 


ye—=Uur —29y=2g(H—y . . (N), 


falls man darin die „Geschwindigkeitshöhe*“ H — ,' 
-9 
einführt. 

Für die Bewegung gilt mit den Bezeichnungen 
von Bild 1 


- dv 
DH-N+l—m =, 
un d ,) 


oder wenn man nach dem Tangentenvektor t und 
Normalenvektor n zerlegt 


dv, i » ..MVUy? 
m + gsinp)t—+|N — Jcosy|n=V, 
de 79 r) t | R ) N 
Daraus folgt 
dv ir yo M vg? Q (2 
— — gsin g; N= +Qcosp . (2). 
di eig u. 


Dabei ist der am Massenpunkt angreifende Bahn- 
druck N positiv, wenn der Richtungspfeil zum 
Krümmungsmittelpunkt weist. Es erweist sich in 
der weiteren Rechnung als zweekmäßig, den Bahn- 
druck auf die Gewichtseinheit zu beziehen. In 
Verbindung von (1) mit (2) folgt dann der dimen- 


h N 
sionslose Bahndruck Q zu 
N 2(H — y) Pr 
—— + COS cd O2) 
Q R r.c6 fi . . . ( . 


Aus der Abhängigkeit des Bahndruckes von der 
(reschwindigkeitshöhe erkennt man, daß N positiv 
und auch negativ werden kann. Im ersten Falle 
wird der Bahndruck von der äußeren Rohr- 
wand, im zweiten Falle von der inneren Wand 
auf den Massenpunkt ausgeübt. Für den Fall, daß 
N>0 ist, wird die innere Wand überflüssig '). 

In der vorliegenden Arbeit sollen nun diejenigen 
Kurven untersucht werden. auf denen bei der be- 
trachteten Bewegung in jedem Bahnpunkt der- 
selbe Bahndruck ausgeübt wird. also N = const ist. 


Aus (3) folgt dann mit const = 
DJ ö 
:E NL —=C I DA 


C folgt aus den Anfangsbedingungen R -- R, 
und 9 = 9% für y=0 zu 


C= 


R, OB: - - ...» (98). 
2. Differentialgleichung für die 
Kurven konstanten Bahndruckes. 
Um die Differentialgleichung für Cartesische 
Koordinaten zu erhalten, wenden wir die bekann- 
ten Formeln 


aUtyW, dı ) 
ee. RER, 


COS 4 
/ ds ] 1+-y’? 
an und erhalten 


2 (H — y)-y” . . . 
t -yi+y? 5) 
[a rise fit... 8 


Das ist die Differentialgleichung der Bahn- 
kurven. Wählt man y als unabhängige Veränder- 
liche und setzt die Differentialquotienten 

1) Man vgl. z. B.: „Vorträge über Mechanik“, I. Teil, 


8. Aufl. von Prof. Dr.-Ing. Kaufmann, — Hannover 19.27 — 
S. dw, 
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dl? r 
dy 1 d’y dy? 
dr dr d.r? hi dzr\?’’ 
dy la 3 
so wird 
2(H — y)-”" — r(1+r’)+C-1-+ 229) — 0 (6). 


3. Lösung der Differentialglei- 
chung. 

Mit Hilfe der Substitution 
It 
Iy 
erhält man nach einigen leichten Zwischenrechnun- 
gen ein erstes Integral von (6) in der Form 


E Sin t & r” zum Evi «> 
( 


A 
PR, 
VH—-y 

dr= ER 4: 250 

m +C 

IYH —y | 
Die Integrationskonstante A bestimmt sich aus 

der Randbedingung „=0 für F Ä t2g. zu 
AL’o 


A=YH-(osp—C). 
Mit diesem Wert für die Integrationskonstante 
erhält man nach einer Umformung 


Ü— cos To 


C 
Y 
Ki 
dz= e :dy .(Ta). 
U >08 e 
1 er cos fo 
Y 
| . 
Hieraus erhält man nach der Formel 
a _ At 
E 4 


folgenden einfachen Ausdruck für den Krümmungs- 
halbmesser 


1) 


5 H R Yy 3 P} 
u CC -cos To ee R,(1 H «+ 8). 
2H 


1. Charakteristische Werte. 
Aus (8) in Verbindung mit (4) erhält man: 
U — c08 9% 


cspg=Ü— DE ° 
Y 
1: 
| H 
und weiter 
I /C- cos?) | 
y—H-|l e Fa . (10). 
Endlich kommt 
Ü "OS ( 51? 
PER 9 „Athen Ge 1 
IK cos | 


In der Zahlentafel 1 sind Ordinaten und Krüm- 
mungshalbmesser gemäß (10) und (11) für einige 
charakteristische Winkel zusammengestellt. 


Zahlentafel 1. 








2 Y R 
| C — cos 9o\°] C — cos gu\® 
oe Jalı-(777 )\ ln) 
1?) 
n C — 608 9,\?] C — cos g,\' \HAI 
s a|ı a u | in. We De 


= 
at 
— 
_ 


Ü — cor go\] ‚„jC — cos 
| C+1 °)| Ru| 3 re 
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>». Lösung der Differentialglei- und erhalten für das umgeformte Integral folgen 


ehung (7a). 
Die Lösung liefert die ge- 


uns 3 Sonderfälle. 


trennt zu behandeln sind. 
1.C >1 (Bahndruck größer als das Gewicht). 
2.C = 1 (Bahndruck gleich dem Gewicht). 
3.0 <1 (Bahndruck kleiner als das Gewicht). 
1.’ Fall: CO >L 
Wir wenden die Substitution an: 
C C—csp_Cz-+1 
i = Er. 
| H 
HC 08 gu)” Ri 2; | 
r zollesu... Jsc\_‘ Entf. 8 sel. 
(? — Dr 1 |C— cos H 
I © 1 2 
+24. 1-17 —— V 1-4 
\Ü — 608 9, H 
an Bier. Y ‚| \ 
+93U-are cos | . an ER: , | B 
es |C — cos y, | H + f 


Da es sich hier um eine Schwerkraftbewegung 
handelt. bei der immer ein Scheitel vorhanden ist. 
genügt es, die Kurven für Ü >1 mit 9, =O (unte- 


den Ausdruck 


2H (C —eos 9,” (Z+C)(Cz +12 1: 
. 3 7 


dz: (C? — Y)Pla ] 1 2? 


Die Ausführung der Integration mit nachfolgen- 
der Rücksubstitution vermittels 


BER; 


C?— 1 \ Yy 
z = C ver | 42 
—C0SYp% H 


das gesuchte Integral von (6) zu 


liefert 


C? 1 Y 1? 
. « = 19 
Bd 2 





Mit 
ren 
chung 


aus (13) die im unte 
vereinfachte Bahnglei- 
Integrationskonstante B 


erhalten wir 
berinnende 
wobei die 


Fo 0 
Scheitel 
(13a). 


rer Scheitel) und die Kurven mit C<1mit | = 7 aus der Randbedingung x 0 für y=-0 zu Null 
(oberer Scheitel) beginnen zu lassen. folgt. 
H | y 1 2 12 
d <CIC+N1 1— ( 1- |(( ) 
cC+.-ya—ı (7 ' H | se 
. (13a). 
F2(+Nn.)/ 1ı-Ic+n V + is ea ern | ı- I cl 
| H | | H f 
Da wir in (13a) eine periodische Lösung vor ben H 14a) 
uns haben, ermitteln wir die Periode der Kurven. u. (C + D?-Vc? — u ER 
Zunächst bestimmen wir die Koordinaten 7. %, .. 
der oberen Scheitelpunkte; hierfür haben wir in Für C>1 folgt aus (14a) a>x. d. h. die 
RE Kurven werden von € — 1 ab aperiodisch, denn für 
(13a) für | 1 = —=- zemäß (9) einzuführen CU’ < 1 wird a imaginär. 
H C+1 
und erhalten 2. Fall: C=1. 
Hier wenden wir die Substitution an: 
= ” ii \ 1 S 1 
%»= ci) ; | —C08 Yo _ 
EU AR, y au 
SC(2n 1) ar 
L, Bar ns n| | a 
(C+1)?-yC?— 1 ) 


Die Periode der Bahnkurven folgt dann aus (14) zu 


z=2H( 


-cosg 0) 


| 2 R PER 
14 I1— cos | H 


B folgt aus den Anfangsbedingungen nach eini- 
gen leichten Umrechnungen zu 


1 „Fo 1 0 = 
B= 2, 00 ar cotg : . (15a). 
— J Y |"'2__ 
r 24H] 2 j 
che } 2H 
Mit C —=1 folgt aus (4a) com —1- ‚man 


0 
. r + - »r ® 
erkennt daraus. daß die Kurven aus (15) für ein 
beliebiges R, nicht mehr mit 9, 0 beginnen 


und erhalten für das umgzeformte Integral die Lö- 


sung 


1/a 1 >) Yy |5 2 \ ” 
1 —..| — . 1-7 —1| +B} (5). 
20 |1— cosp, | H Bali 
Aus dem bereits erläuterten Grunde beschränken 
wir uns auf den Fall 9% = = und erhalten die 
Gleichung für die im oberen Scheitel (einzigen 
Scheitel) beginnende Kurve in der Form 
| Y ne | f 
1—; 1 (16). 
] H 1-8 
können. da dann obıge Gleichung nicht erfüllt 
wäre. Für RR >» wird cosg, —1 und nach (11) 
Ro = &: das bedeutet aber. daß die Bahn- 


kurve eine horizontale Gerade ist, für welche der 
Bahndruck gleich dem Gewicht, oder C —=1 ist. 
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3.F:;: .6 ö a r N 
4 1 = Su man aus den Wurzelzeichen Y— 1=i heraushebt 
Hier ist eine erneute Ausrechnung überflüssig. und für arceos ||] iArC&oj[] setzt. die 1.ö- 


denn man erhält unmittelbar aus (13), indem sung 














H -(Ü — 608 gu)? 1— (0? | “ 12 
ke sr fe-|. = E ) +e) | Wi ER ep 
( .— )>'e [LU —C0S@ H | (& cos ya = °® | 
vo 1? = 
+] | 1 Yıc| =. ct ARE, 
Ü 608g, H | 
er Y | 
+ 30 -Ar&of |; | 1- Z+0C|+B 
( COS fo H | | 
wobei B aus den Randbedingungen folgt zu 
1 —( cos sing, : Tl UcCoSsgQ, 
= rn el Fe, Ban ac .,,. ee 
| ( cos yo | U — COS Yo ( COS Yo 
Die Bedingung für den Scheitelpunkt lautet nach ’ 4-14 F ar) . | 
(Th 
k Ys Ü-C0SYyo } H(Ü— cos po)" 
(12) mit« . R I ; damit folgeı 2; B- 
ee. uch ee 5 u (1 — ©) | 
Wir beschränken uns auch hier auf die Kurven. 
die Koordinaten des Scheitelpunktes zu welche mit 90 = = beginnen und erhalten hierfür 
mit B=0 
H Ä ıy/ ar: RE 
1-0%-yı-ca\ | H\) | H\ 
IS). 
| ] 
2(C +9 ae a; ER ..‘& 7 WORRE . 1 +30-Ar Ef C+ (1 co 1 y\ 
H | Zul) 
Für cosp = wird nach (10) U — 0; das be- Ar rn 
deutet aber, daß die Kurven für C<1 eine Asym- ZU4N___17 __Äohentmie 
ptote (Tangente an den unendlich fernen Punkt) I Sl iv”)  Zeitgeradt 
haben. deren Neigungswinkel gegeben ist durch I] IL, N 
Yo —arccosC. Die Kurven für C—=1:g, 0 ! | vi” Nr 
.® : : . “ | f - 
und © =0:9,.=, bilden die beiden Grenzfälle. x . x“ 
2 Ss ı % 
Sonderfall: C=0. S | R 
Aus (17) erhält man mit C=0 und B=2tg m Ri 5 
; rg > 
. / Y ’ 
zr=2Hcos’g,— eo? - rigo)|- I } en 
To / tg” oo H cos? gu Yo My / = %, 
we [aa] Js 
Drücken wir y durch x aus. so erhalten wir, in- Bild 2. 
: v, . 
dem wir noch H=,_ setzen: u 4 
29 y=Hsn®’p=35 sin? go: 
e g2° va > 
yz{ı\W89— 2 02 608? g, ee De |. z,—=2H sin 9, C08 Yo = = SINYGCOS go. 
Die Kurven für C — 0 sind somit Wurfparabeln. AR ’ 
2 H f Hohen ME 
Wie man aber sofort erkennt, kann € = ( R, dy ag 
+ eos, nur für 90 > , bestehen, da sonst € nicht 
Null werden könnte. Die Geschwindigkeitshöhe 
folet daraus zu 
R, cos « R, €08 po 5 2 
H=—- ut, siehe Bild 2. 
Die Geschwindigkeitshöhe ist mithin die Hälfte Bild 3. | 
der Projektion des Krümmungshalbmessers auf die 
Vertikale durch den Krümmungsmittelpunkt M. 
Weiter kommt für die Koordinaten des Scheitel- \ 
punktes der Wurfparabeln aus (17 b) mit © — O und re 4 
B=2tg% A1983 
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Der Krümmungshalbmesser im Parabelscheitel 
ist bekanntlich der Parameter, also folgt nach (12) 
mit q 7 


8, =P: - R, cos? g, = Ro cos? yo. 


’ 
R, COS To 
) 


Mit H erhält man für die Abszisse des 


Scheitelpunktes weiter 


2, = — Ro: Sin Yo cos? go: 
Man findet demnach die vertikale Achse der Pa- 
rabel. indem man die Vertikale durch P mit der 
Horizontalen durch M im Punkt B zum Schnitt 
bringt und von B eine Senkrechte auf den Krüm- 
mungsradius fällt. Der Punkt 4 ist dann ein 
Punkt der Parabelachse. denn es ist 


2, = BE - 


(R,CoS Yo )-SIN Yo "COS go 
R, sin v o cos? 7 B% 


eine wenig bekannte Parabelkonstruktion aus den 
Bestimmungsstücken: Krümmungsmittelpunkt. Rieh- 
tung der Parabelachse und Parabelpunkt mit Taı- 
ventenriehtungz in diesem Punkt. 

Brinst man die Tangente in P mit der Parabel- 
achse im Punkt E zum Sehnitt und halbiert die 
so erhaltene Strecke DE (DS=ES),so ist S der 
Scheitel der Parabel. Als letztes erhält man den 
l’arameter » der Parabel in der Strecke AD. da 
bekanntlich die Subnormale der Parabel konstant 
und gleich » ist. Es folgt weiter das interessante 
Ergebnis. daß die Höhenlinie die Leitzerade der 
Parabel ist. denn es ist 





. I 1 
H.cos?’y =; Recon =». 


LS=H— 5 


Im  Scheitee der Parabel ist 
Urs == 0, COS pP —=const; d%y, muß aber genau so 
sroß sein, als wenn der Massenpunkt die Höhe 
(H — y,) frei durchfallen hätte. Es muß daher sein 


Uys 4 7 


Uxs 27 ) x ee 
H -y, 29 BE huge 9a,=#-008 g,; 


oder 
= H-sin? yo. 

6. Graphische Konstruktion der 
Kurven. 

Ist in einem Kurvenpunkt ? die Richtung der 
Tangente ? bekannt. so kann man die Kurven 
konstanten Bahndruckes nach folgender Konstruk- 
tion graphisch darstellen: man trägt auf der 
Kurvennormalen zum Krümmungsmittelpunkt zu 
die willkürliche Strecke PD -— 2a und vom Krüm- 
mungsmittelpunkt weg die Strecke PR a auf. 
Dann fällt man vom Punkt B ein Lot auf die 
Ordinate durch P, so daß PA — acosp wird. Von 
A nach aufwärts trägt man die Strecke AC = aU 
ab; damit wird PC—=aC-—acosg 
Nun zieht man durch E eine Parallele zu CD und 
erhält im Schnittpunkt F den Krümmungsmittel- 
punkt für den Kurvenpunkt P; EF\CD,PF=R. 











a(C — COS P). 





Beweis: 











PR PC 
PF PD’ 


oder 
H—y _ a(C—cosp) _U—cosg 


RB: 3% 2a 2 


oler 
2(H Yy 
R 


COS 1 Ü 


siehe Gl. (4). 


Man kann nun die Kurve nach dem Krümmungs- 
kreisverfahren durch Kreisbogenelemente ersetzen. 
indem man in jedem neuen Kurvenpunkt obige 
Konstruktion zur Bestimmung des neuen Krüm- 
mungshalbmessers anwendet. In der Nähe eines 
Scheitels fallen Ordinaten und Kurvennormalen zu 
sammen; man wendet dann diese Konstruktion 
zetrennt. am besten in einem rechtwinkligen Drei 
eck an und überträgt den gefundenen Krümmungs 
halbmesser in die Zeichnung. 

Die graphische Konstruktion gilt für alle Werte 
von €: ist C, R, und 9, gegeben. so findet man 


U— cos g, 
aus (da) H= es; 


S R,. In Bild 4 ist nach dem 


xraphischen Verfahren eine Kurve für U —3 ze 
zeichnet. wobei zur Kontrolle der Genauirkeit 
Koordinaten und Krümmungshalbmesser im oberen 
Scheitel gemäß 


Ü 3 
R, Rul,. 


nachgeprüft und eine sehr gute Übereinstimmung 
festgestellt wurde. Bild 5 zeigt die nach (16) be 


rechnete Kurve für ( 1, während Bild 6 auch 
.) 
zwei nach (18) errechnete Kurven für Ü — , und 
.) 
. a EVER A 
C=, zeigt. Bill 7 zeigt vier Kurven für ( 0 


und (Yo)ı 120°, (2) = 135°, (90); = 150° und 
Fo) 180°. wobei die Konstruktion nach Bild 2 
angewendet wurde. 


7. Zusammenfassung. 


befaßt sich mit der 
konstanten Bahn- 


Arbeit 
Kurven 


Die vorliegende 
Untersuchung der 
druckes. Die aus der Bedingung C erhaltene 
Differentialgleichung wird zur Erzwingung einer 
zeschlossenen Lösung so umgeformt. daß y die 
unabhängige Veränderliche wird. Die Lösung ent- 
hält drei charakteristische Fälle: 


1.für C >1 folgen periodische Kurven, deren 
Periode für kleiner werdendes C immer größer 
wird, 


2.für C=1 folgt eine Sonderlösung, wobei diese 
Kurven den periodischen vom aperiodischen 
Bereich abgrenzen, 


3.für C<1 folgen aperiodische Kurven mit 
einer charakteristischen Eigenschaft einer 


Pi arccos (Ü ist. Ein 
0. die als 


Asymptote, wobei 
Sonderfall sind die Kurven für C 


Wurfparabeln erkannt wurden. für welche 
wi, 
Io —= ist 


Zum Schluß wird noch ein graphisches Verfah- 
ren entwickelt, wonach die Kurven für beliebiges 
C sehr einfach gezeichnet werden können. 

Karlsruhe. Johann Fadle VDI. 198 
z. Zt. im Felde. 
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Zum Verfahren der schrittweisen 
Näherung. 


Es wird im folgenden eine Lösung der Differen- 
tialgleichung Y” + F(2)y=0 unter gewissen Vor- 
aussetzungen über F(®) nach der Methode der 
schrittweisen Näherung berechnet; für die so er- 
haltene Reihe werden bestimmte Eigenschaften 
ihrer Glieder festgestellt, und es wird eine Rest- 
abschätzung gegeben. 


Vorgelegt ist die Differentialgleichung 


te ee A; 9 


F (x) sei in einem abgeschlossenen Intervall (a, b) 
stetie — und damit beschränkt — und von einer- 
lei Vorzeichen (wobei Nullstellen zugelassen sind; 
es soll jedoch Fix) in keinem Teilintervall iden- 
tisch Null sein). Es werde die Lösune (x) mit 
den Anfangsbedingungen 


A EE 


e ist eine Stelle aus (a, b mit dem Ver- 
fahren der schrittweisen Näherung berechnet. Als 
\usgangsnäherung wird entsprechend den Anfangs- 
bedingungen , = 1 gewählt. Man erhält für Y (7) 
lie Reihe 


Y = Tr Yı + 2 A 


mit Yu 1 und 
(2) 
x 
Yu \\F&-„.-, dd&dn,n=1,2,3---- 
ec 
Nach einer wohlbekannten Schlußweise ergibt 
sich 


Yul)>V für nox 


gleichmäßig im Intervallazr sb. 
Die Glieder der Reihe (2) haben folgende Eigen- 


schaft: 
Yn\X) Yn-ılf) 
4 (X) [ olT) ” 
r In-ı\ In-e\ et (3). 
fürrasz<c 
und c<r<b, n=2,3,4---- 
Der Beweis werde nur für e<r=zb durch- 


geführt, da er für as x<c analog verläuft, und 
läßt sich durch vollständige Induktion führen. Für 
2 folgt die Richtigkeit auf Grund der Mono- 


tonie von | Yı (X) 
x 
y))=\\I FO yld)ididn 
+ a 
an 
n ’i{E 4 rn 2 » 
Yıladı-ı \ Fe&e)dSdn=yi(e). 
cc 


/weimalige Anwendung des  verallgemeinerten 
Mittelwertsatzes der Differentialreehnung auf den 
Quotienten 

Yn \€) 

ı (ec) 


Yn (*) 
(7) 3 Yn 


Yn\r) a 
Yu-ı (7) Yn 
liefert: 
Yn (.£) 2 Yn (Ak) M Yn (k) u; Yn (ec) h Yn (d) 
Yn-ı\?T) Yu ı (Ak) Yıh ı (Ak) — Yı ‚(c) Yu ( 
Dabei ist 
EER<2. c<iIich; 


Aus (2) folgt 


Bd.21 Nr.2 Apr. 1941 
MU) _Mm-iı (WM 
In-ı) Yn-s(d' 
also: 
Yn\X) Yn-ı (Ü) 
YIn-ı(?) Yn-2(D 
W BR 3 “ Yn-ı (8) e 
enn gezeigt werden kann, daß „, eine 
n-a\ 


monoton steigende Funktion ist. gilt 


Yn-ı  ‚YIn-ı\E) 
Ye ld)  \..52) 
und damit 
Yn(z) Bi: N. 
Yn-ı(E) Yn-2\#) 


Der Schluß von (rn -1) auf n ist damit durchgeführt. 
Es wird daher im folgenden das Vorzeichen der 
Ableitung von 


Yun 1 (Fr ) 


Yn-2 (FT) 


mit Hilfe einer Identität bestimmt. Dazu wird der 
folgende Ausdrvek unter Benutzung von (2) um- 
geformt: 
x 
\F-yu 3’In ‚ds 
5 
X 
m, Yı 2"’In f ds 
$- 
2 


’ 1 ’ ’ - 
= Yn-ıIu-2 7 \Yn-ı'Yn-2d8 


( 
x 


= 4 ' i ” LE 
ai Yn-ıIn-2Tt Yn-ı' Yu \Yn-2'"Y4n-ı Us 


( 
N e 

277 \Fy, „ds 
ce 


’ N ’ 
In-ı Yn-2T YIn-ı Yu 


Nach Zusammenfassung der beiden Integrale er- 
gibt sich: 

x 

\ F-(yn 3’ Yn-ı - Ya „)ds (4). 


c 


mn ’ ’ 
= Yn-ı' Ya  Yn-ı Yn-a 


Ist F (x) im Intervall (a, b) nichtnegativ, so haben 
die 4;(r) wechselndes Vorzeichen. und es ist 


Yn-ı| _ __Yn-ı 


Yn 2 In ® 
Die Ableitung wird gleich 
m * ı"Yn-2— In-ı"Yn-: 
Yn-2 
Unter Verwendung von (4) folgt: 
x 
Yn-ı (2) 1 [ ‚ ’ e 
=— F-(4 1 — Y-„)dE. 
Yn »(£) Y, “ Yn-s’ Yn-ı YUn-» ds 
2 


Aus der bis (#-1) als richtig angenommenen Un- 
gleichung (3) ergibt sich, daß 
Yn-s’I/n-ı Yn »<Q— 
ist. Das Integral ist infolgedessen negativ und so- 

Yn 1 (E) ER 
positiv. 
Yn-a\X) 
Für im Intervall (a, b) nichtpositive F (x) erhält 
man das gleiche Resultat. 


mit die Ableitung von 





YıUM 


ne 


rt. 
ler 


m- 


Jen 
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Le) 
or 





Yn-ı (2) 
Yn-2(£) 
tion, was zu zeigen war. 

Nach (3) ist die nach unten beschränkte Folge 


ist also eine monoton steigende Funk- 


Yı Ya Ya 

„| mn! | 
bei festem = monoton abnehmend. Sie hat also 
einen Grenzwert: 


WET 


Yi (x) ni 


lim 
i>« Wi-ı(2) 


Hierin ist: 


BEER. nr es 
Denn es ist 
YIn|__ Yı : 2 ER Yn n Yu E gqn 
Y 0 Y o Y 1 Y n—i 4 n 1 


und für n > geht die linke Seite gegen Null. 
Mit Hilfe der beiden Ungleichungen (3) und (5) 
läßt sich eine Restabschätzung der n-ten Partial- 
summe der Reihe (2) geben. 
Es ist 
Y =Y + Yı ++ + Yu Ru +1» 
wobei 


Rna+ı =Yn+1ıT+Yn+24+°°°° 


re ” Y 
Dabei sei n so gewählt, daß " 


< 1list;nach 


IN ' 
(5) ist das bei festem @ +. möglich. Dann läßt 
sich R, +, mit (3) folgendermaßen abschätzen: 


Ru +1 
la ı Inte! ,  Yn+s ] 
Yu+ıl)l+ $ ER, 
DEN | Yn+i Yartı 
| U P} Ya D ih | 
Yatı 1 m Hn+2 4 Yn+s li Yn+? es 


HIn+ı 3 Yn+a Yn+ı 


< A| im FRE 


Yn-ı Yn-ı Hu 1 
| 1 
Hu 1 1 er 
In ı 


Es hat sich also ergeben: 
Ya) 
Un 1 (2) 
Yu (X) (6). 


Ya-ı(?T) 


Ru+ıld)! < | yu(®) 


für jedes r+c aus (a,b). 


Diese Formel gestattet, mit dem letzten und vor- 
letzten Glied der n-ten Näherung 


Yt YıtYt''+t Yu 
den Rest abzuschätzen. 
Ya) 
Yn-ı(2) 
ergibt sich aus (6) leicht eine Abschätzung, die für 
alle aus (a, b) gleichmäßig gilt: 


Infolge der Monotonie von und 4 (x) 


Ru+ı(2)| < Max 


Yn (a) Yı (b) 
(a) / N. 
Yna) In -.\e . Yun (b) In-ı (b) ( ’) 
Yn (a) Yn (b) 
YAn-ıla) Yn-ı (6) 


Es soll noch die gegebene Restabschätzung mit 


der üblichen verglichen werden, die sich aus der 
Beschränktheit von F (x) herleiten läßt: 


F(x)I<M in asesb, 


Ru+ı(%) | 
(rt — e)?n+? 1 
< Mn+1. X (7); 
‚ (2n-+ 2)! M (2 — c)? 
2n +3) (2?n—+ 4 


Ein Vorzug von (6) besteht darin, daß die Ab- 
schätzung auch möglich ist für Anfangsbedingun- 
gen an den Stellen 2= +. Die Güte der ze- 
gebenen Restabschätzung möge an zwei Beispielen 
nachgewiesen werden. 

Vorgelegt ist die Gleichung 4” +2-y=0, Ge 
sucht ist die Lösung mit den Anfangsbedingungen 


y(0) 1, K(0)=d0. 


Das Näherungsverfahren liefert: 


r> +6 
yval setag 56 
zn 


ie Ru+,(2) 


2.3.5.6- + -(3n--1)-3n ' 

Die Restabschätzung (6) ergibt die obere 
Schranke: 

Ryu+ ı(?) 
r 3n--3 1 

a e 

2.3.5.6 ++ -(3,_,)Bn-(In—1)-In 2 
für | 2 |< (83n-—-1)-3n. 


Abschätzung (7) liefert: 
r 3NT3 1 
Rn+, (z)|< —- . 
ht (2n-+2)! x |? 
(2n—+3)(2n—+4) 
für | 2) <(2n+3)(2n +4). 
Man erkennt, daß die erste Schranke kleiner ist 
als die zweite. 
Ein weiteres Beispiel zeige den numerischen 
Unterschied der beiden Abschätzungen. 
Es ist die Lösung der Gleichung 
„ ke? 


4 1 4 0 
ud (1- Pu 


mit den Anfangsbedingungen 
y)=1, y(O)=0 


anzugeben. Die strenge Lösung lauset: 


y=Vi+a:cos(y1 


Die zweite Näherung yy = 4 + Yı + Y, ist nach 
der Substitution 2==tgg gleich 


k?-aretg.r). 


k? Ai 
snm1l-zrtgEpr 


8 


tgP—p?). 





Restabschätzung 
T Y Yı Y— Yı - = 
nach (6) nach (7) 





1.0000] 0,6279] 0,6284 1—-0.0005] 0.0012 0.0014 
2,1445 |— 0,0794 | - 0,0732 | - 0.0062] 0.019 0,147 
3,0777 |- 0,6630 |— 0,6477 |- 0,0153] 0,05 1.43 














In der Zahlentafel sind an drei Stellen die Funk- 
tionswerte (für k = 1) und Abschätzungen nach (6) 
und (7) zusammengestellt. 


Prag. J. Bitterlieh-Willmann. 237 
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über Schwerpunkts-Parallelweiser- 
Tafeln. 
Wir gehen aus von drei Gruppen von m bzw. 
n bzw. r Skalen 


x! - H mg . (Su ), yı — I m w.\s ı) | { 
3); 
w=1,2,.--,m | 
A HM; 
7 Kn q " (9.). y"! Ln y, N») | " 
02) 
„——L2, +++, 
I Kr 111 £ a. Lı y Im (£&o) \ = 
(8 
0 1.2 ! 


wählen darauf m bzw. n bzw. r Punkte (gegeben 
durch u bzw. 7, bzw. £e) und bestimmen zu jeder 
der drei Punktgruppen den Schwerpunkt 


I I I I 


T En ze YıHt En ze 7) 
Pl: Pe l nm } y! u 

m " m 
pi . rl yH 

n F n 

zUN N z yu' ya 

pi. pi „ya 

r ; r 
Hierauf lege man dureh Pl zur Geraden PIl pi 


die Parallele und bestimme deren Sehnittpunkt 
P ic”, y®) mit einer weiteren Resultat-Skala 


r"’—=H g (d ,; PP J yp ( u) (V). 
Alsdann besteht wegen der geometrischen Schlüssel- 
gleichung 

au 


ii 


y' yl ie u 
x! ‚rl ri 


unter der Voraussetzung 
HiIE—=K/L 


zwischen der gesuchten Variabeln 9 und den ge- 


gebenen Veränderlichen £&u, 7,. £e der Zusammen- 


hang 
y(®) 2 Y (Zu) b yi (m) — 2 ya (£e) 
q (9) -2 ol (Eu) IE) - Eye) 


der durch eine bemerkenswert große Anzahl von 
unabhängigen Veränderlichen ausgezeichnet ist. 

Zwecks praktischer Verwendung empfiehlt es sich, 
die Skalen zu spezialisieren: Man lege die Skala (0) 
auf die y-Achse, 7° = 0, die Skalen (2) geradlinig 
und vertikal mit zU=0, die Skalen (3) geradlinig 
und horizontal mit „N=0. Dadurch erhält man 
eine Darstellung für den (immer noch sehr allge- 
meinen) Zusammenhang 


„ I 

\ ’ J 
(9 2 - < u m) 
ya)=<yulu) T-9uH) Si, 
-dY, (Se) 


Weitere nahe. Z.PB. 


Spezialisierungen liegen B 
kann man, m = m’—+ m” gesetzt, die Skalen (1) 
für u=1,2,---, m’ auf vertikale Geraden der 


Mittelabszisse Null und für u=m’+1,. m’+2, 
+, m’+-m’ auf horizontale Geraden der Mittel- 
ordinate Null legen, wodurch man einen Zusammen- 
hang erhält, in dem vier voneinander unabhängige 
und vollkommen getrennte Gruppen von m’, m”, n,r 
Funktionen je einer Variablen auftreten. 

Bei der großen Allgemeinheit der hier darge- 
stellten Zusammenhänge darf man sich nicht wun- 
dern, wenn die Darstellung nicht so bequem ist, 


wie die Nomogramme für nur 2 oder 3 unab- 
hängige Veränderliche. Jedoch dürite es sich 


zwecks leichter Bestimmung der Schwerpunkte emp- 
fehlen, die Zahlen m, n,r sämtlich =4 oder allen- 
falls noch =6 zu wählen!). 

Zum Schlusse sei darauf hingewiesen, daß man 


an die Punkte ?, Pl, PU, Pl statt des oben be- 
nutzten Parallelweisers auch einen Kreuzweiser 
anlegen kann: Man erhält dann Schwerpunkts- 


Kreuzweiser-Tafeln. 


Münster i. W. ludwig Neder. 268 
1) Auch ist es zweckmäßig, die 4 Skalengruppen (0) bis 
(3) in 4 verschiedenen Farben zu zeichnen. 


BUCHBESPRECHUNGEN 


Handbuch der Werkstoffprüfung. Herausgegeben 
unter besonderer Mitwirkung der Staatl. Material- 
prüfungsanstalten Deutschlands, der zuständigen 
Forschungsanstalten der Hochschulen, der Kaiser- 
Wilhelm-Gesellschaft und der Industrie sowie der 
Eidgenössischen Materialprüfungsanstalt Zürich. 
I. Bd.: Prüf- und Meßeinrichtungen, bearb. von 
R. Berthold, Berlin, A. Eichinger, Zürich- 
Düsseldorf, W. Ermlich, Berlin, G. Fiek, 
Berlin, L. Föppl,. München, R. Glocker, 
Stuttgart, E. Lehr, Augsburg, E. Siebel, 
Stuttgart, OÖ. Vaupel, Berlin, herausgeg. von 
Prof. Dr.-Ing. E. Siebel, Techn. Hochschule 
Stuttgart. XIV + 658 S. m. 763 Abb. Berlin 1940, 
Verlag Julius Springer. Preis geb. 69 M. 

Das vorbildlich ausgestattete Werk ist ein zu- 
verlässiger Ratgeber für den Werkstoffprüfer und 
darüber hinaus eine Fundgrube für den Versuchs- 
ingenieur, den Meßgerätebauer, den technischen 
Physiker und den Maschinenkonstrukteur. Die 
ausgezeichnete Darstellung wird durch Schema- 
und Schnittzeichnungen, durch perspektivische 
und Lichtbilder vorzüglich unterstützt. 

Der umfangreiche Stoff ist nach den verschie- 
denartigen Prüfverfahren unterteilt und jedes ge- 
schlossen von einem maßgebenden Fachmann dar- 
gestellt: Grundlagen der Werkstoffprüfung (Sie- 
bel), Prüfmaschinen für ruhende Belastung (Fiek), 
für stoßartige Belastung (Lehr), für schwingende 
Belastung (Lehr), Sondergeräte zur Prüfung von 


Federn, zur Härteprüfung, zur technologischen, 
Verschleiß-, Abnutzung-, Lager- und Drahtseil- 


Prüfung (Eichinger). Weitere Abschnitte befassen 
sich mit der statischen und dynamischen Deh- 
nungsmessung an Konstruktionsteilen (Lehr), mit 
der spannungs-optischen (L. Föppl) und röntgeno- 


graphischen Spannungsmessung (Glocker), mit 
der zerstörungsfreien Werkstoffprüfung mitteis 


strahlender und magnetischer Durchdringung 
(Vaupel). Beachtlich ist auch die Abhandlung über 
die Geräte und Methoden zur Nachprüfung der 
Kraftanzeige und zur Erhöhung der Genauigkeit 
(Ermlich). Alles in allem ein grundlegendes Werk. 


3raunschweig. G. Niemann. 235 
Dr. LUDWIG HÄNERT, Geschütz und 
Sehuß, eine Einführung in die Ge- 
schützmechanik und Ballistik. Dritte, 


verbesserte u. erweiterte Aufl. VIII +418 S. m. 
177 Textabb. Berlin 1940, Verlag Julius Springer. 
Preis geb. 27 M. 

Die erste Auflage ist in dieser Zs. im Bd. 9 (1929) 
S. 253, die zweite im Bd. 16 (1936) S. 253 bespro- 
chen worden. Die Gliederung des Inhalts ist im 
eroßen ganzen die gleiche geblieben (Il. innere 
Ballistik, II. Geschützmechanik, III. Außenballi- 
stik), jedoch sind die einzelnen Abschnitte bei un- 
veränderter Grundhaltung des Gesamtwerkes um- 
gearbeitet, ergänzt und besonders nach der phy- 
sikalischen Seite hin vertieft worden. 
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Im Abschnitt I wird u. a. jetzt auch die Gestalt 
der Drallkurve unveränderlichen Leistendruckes 
bestimmt, wenn der Pulvergasdruck in Abhängig- 
keit vom Geschoßweg im Rohr bekannt ist und 
die Reibungsarbeit berücksichtigt wird. 

Die Geschützmechanik bringt neu bzw. erweitert 
die Stabilitätsbedingungen beim Abschuß, Betrach- 
tungen über die Mündungsbremse, einen Abschnitt 
über die Bewegung einer federnd gelagerten Masse 
unter dem Einfluß einer Kraft mit bekanntem 
zeitlichem Verlauf und einen Abschnitt über die 
Modellgesetze der Physik. 

Am wenigsten verändert ist der außenballistische 
Teil. Neu aufgenommen wurde hier die von Popoff 
angegebene Lösung der Bewegungsgleichungen 
des Geschosses unter Benutzung eines schiefwink- 
liven Achsenkreuzes. 

Das Buch ist auch in der neuen Bearbeitung 
seinem ursprünglichen Ziel treu geblieben, leicht 
verständlich und zuverlässig in die mathemati- 
schen und physikalischen Grundlagen der neu- 
zeitlichen Waffenlehre einzuführen. Bei der augen- 
blicklich zu beobachtenden (und hoffentlich auch 
später anhaltenden) Anteilnahme an wehrtechni- 
schen Fragen in Kreisen der Industrie und der 
Hoch- und Fachschulen ist dem Buch weite Ver- 
breitung zu wünschen, damit es so mithilft, dem 
Nachwuchs das Verständnis für die wissenschaft- 
lichen Grundlagen der Wehrtechnik zu erschließen. 

Daß sich das Hänertsche Werk in der vom Ver- 
fasser getroffenen Auswahl des Stoffes und bei 
der ihm eigenen Klarheit der Darstellung bewährt 
hat, beweist die rasche Aufeinanderfolge der drei 


Auflagen. 233 
Essen. K. Stange VD!I. 


Dr. techn. VACLAV HRUSKA, VenkovniElek- 
triekä Vedeni, Po&@itanä Jako Pruznä Rete- 
zovka. (Elektrische 
rechnerisch als elastische Kettenlinien 
aufgefaßt). 60 S. Prag 1940, Verlag Elektro- 
technicky Svaz Üeskomoravsky. Preis brosch. 
24 Kronen. 

Gewöhnlich werden Hochspannungsleitungen 
unter Zugrundelegung des Bildes unelastischer 
Kettenlinien berechnet, wodurch sich gegenüber 
dem wirklichen Verhalten Abweichungen sowohl 
hinsichtlich der Spannungen, als auch besonders 
hinsichtlich der Gestalt der Seilkurve ergeben. 


Die Mitbeachtung der Elastizität bedeutet aber 
eine beträchtliche Komplizierung des Problems. 


Diese nun meistert der Verf. mit sichtlichem Er- 
folg durch Einführung eines von der Spannung im 
Tiefstpunkt abhängigen Parameters, wodurch es 
ihm gelingt, sämtlichen Gleichungen für die Span- 
nungen und Koordinaten eines laufenden Seil- 
punktes eine sehr übersichtliche und für die nume- 
rische Auswertung geeignete und brauchbare 
Form zu geben. 

Durchgängig werden zwei für die Praxis beson- 
ders wichtige Problemstellungen behandelt und als 
Probleme I und II unterschieden, nämlich die Be- 
stimmung aller Spannungen und der Gestalt aus 
der Forderung, daß unter den ungünstigsten Wit- 
terungsverhältnissen (Reif- und Windbelastung) 
die höchste Spannung am oberen Aufhängepunkt 
den vorschriftsmäßigen Höchstwert der zulässigen 
Beanspruchung gerade erreicht (Problem T) bzw. 
Ermittlung der „Montagespannungen“ unter Be- 
rücksichtigung der zur Zeit der Montage ,herr- 
schenden Witterungs- und Temperaturverhält- 
nisse, welche aber so zu erfolgen hat, daß bei 
späteren evtl. viel ungünstigeren Verhältnissen 
doch niemals die zulässige Höchstbeanspruchung 
überschritten wird (Problem I). 

Nachdem Verf. im ersten Abschnitte aus der 
Betrachtung des Gleichgewichtes des Seilelementes 


Überlandleitungen. 


die allgemeinen Grundgleichungen der elastischen 
Kette in der oben erwähnten Parameterdarstel- 
lung gewonnen hat, betrachtet er im zweiten Ab- 
schnitte die unelastische Seilkurve und zeigt, wie 
die numerische Auswertung der bezüglichen Glei- 
chungen durch (von ihm bereits an anderer Stelle 
veröffentlichte !)) Nomogramme weitgehend unter- 
stützt werden kann?). Auch die Relationen zwi- 
schen den Hyperbelfunktionen werden, soweit 
später benötigt, angegeben. Wo die Genauigkeit 
der so gewonnenen Ergebnisse nicht ausreicht, 
bedient sich Verf. des am gegebenen Orte klar 
erläuterten Iterationsverfahrens, dessen einzelne 
Schritte in Berechnungstabellen übersichtlich an- 
gegeben werden und dessen Konvergenz durch 
eine geschickte Modifikation sehr beschleunigt 
werden kann. Im dritten Abschnitte wird ebenso 
die elastische Kettenlinie behandelt und die durch- 
gerechneten Beispiele dürften dem praktisch täti- 
gen Ingenieur eine anschauliche Anleitung dazu 
geben, wie er in ähnlichen Fällen vorzugehen hat. 
Im vierten Abschnitte wird der Einfluß an den 
Enden der Seilkurve etwa gelenkig angeschlosse- 
ner als starre schwere Stäbe gedachter Isolatoren 
in analoger Weise berücksichtigt, und so darf 
wohl behauptet werden, daß die interessante 
Schrift jedem, der sich mit ähnlichen Fragen zu 
beschäftigen hat, großen Gewinn bringen kann. 
Bet einer Neuauflage wäre vielleicht eine weiter 
sehende Entwurfsbegründung der Nomogramme 
empfehlenswert (z. B. Erklärung der Gl. 2, 46 auf 
S. 14), damit auch jener, dem weitere Literatur 
nicht zur Verfügung steht. das Nomogramm leicht 
entwerfen kann. In Gl. 2, 43. S. 14. fehlt bei 4 
im Nenner der Faktor q und in Abb. 3, S. 18, 
wäre die Ordinate b des Punktes Ms anzugeben. 
Prag. K. Karas. 3281 
', Laska-Hruska: Pocet graficky a graliekomeeha 
nieky. Prag 1922. 
2) Hruska-Kelbieh: Pocitaci pravitko k urceni pru 
hybu venkovnich elektriekyeh vedeni, Elektroteehn Obzor 
Bd. 25 (1986) S. 661. 


Grimsehls Lehrbuch der Physik, 
neubearb. von Prof. Dr. R. TOMASCHEK, Direktor 
d. Physik. Instituts d. Techn. Hochschule München. 
II. Bd.: Elektromagnetisches Feld/Optik. 9. Aufl. 
X+867 S. m. 1209 Abb. und 1 farbigen Tafel. 
Leipzig und Berlin 1940. Verlae B. G. Teubner. 
Preis geb. 26 M. 

Das Grimsehlsche Lehrbuch der Physik in der 


Neubearbeitunge von Tomaschek hat sich längst 
unter den physikalischen Lehrbüchern einen 


festen und allgemein anerkannten Platz errungen. 
Infolge der anschaulichen Behandlungsweise auch 
schwieriger Gebiete findet der Studierende hier 
eine wertvolle und gwlückliche Ergänzung vieler 
rein formal mathematischer Darstellungen, wie sie 
in den Lehrbüchern der theoretischen Physik 
üblich und oft unvermeidbar sind. 

Die vorliegende 9. Auflage ist gegenüber der 
vorhergehenden im großen und ganzen unver 
ändert geblieben. Im einzelnen enthält sie jedoch 
eine Reihe von kleinen Verbesserungen. So ist 
der Abschnitt über die Eigenschaften starker 
Elektrolyte vollkommen umgearbeitet und der Ab- 
schnitt über Ferromaenetismus auf den neuesten 
Stand gebracht worden. Ferner sind die wich- 
tigeren technisch-physikalischen Hilfsmittel der 
Kernphysik und anderer Gebiete neu aufgenommen 
worden. 

So wird jeder, der sich auf dem Gebiet der all- 


gemeinen Physik und ihrer wichtigsten Anwen- 
dungen umfassendere Kenntnisse aneienen will, 


dieses Buch mit größtem Nutzen durcharbeiten. 
Dresden. H. Stuart. 31 
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Nachrichten 


Z. angew. Math. Mech. 
Bd. 21 Nr. 2 Apr. 1941 





THEODOR VAHLEN, Abstrakte Geome- 
trie. (Deutsche Mathematik, Beihefte, im Auftrage 
der Deutschen Forschungsgemeinschaft herausge- 
geben von Theodor Vahlen. Zweites Beiheft, 1940.) 
2. neubearb. Aufl. X + 224 S. m. 92 Abb. Leipzig 
1940, Verlag S. Hirzel. Preis kart. 8 M. 

Die Prüfung des logischen Aufbaues der Geo- 
metrie bedarf einer Anordnung von Grundsätzen 
(Axiomen) und aus diesen abgeleiteten Lehrsätzen. 
in der kein Grundsatz von den übrigen abhängt 
und jeder Lehrsatz gerade an der Stelle eingereiht 
ist, an der er bewiesen werden kann. Die hierbei 
nötigen Nachweise müssen es streng vermeiden, an- 
schauungsmäßige Begriffe hinzuzuziehen, ehe sie 
logisch eingeordnet sind; und dazu dient mit be- 


sonderem Erfolg das Verfahren, das von dem 
geometrischen Inhalt der Sätze absehend sie in 


Aussagen über Systeme nicht-geometrischer Dinge. 
insbesondere über Zahlensysteme, umdeutet. So 
entsteht die „Abstrakte Geometrie“. An ihren An- 
ang gehört Teil I eine Darstellung der 
„Grundlagen der Arithmetik“, in der von Mengen, 
von Zahlensystemen und von Größensystemen ge- 
handelt und der für die Geometrie so wichtige, all- 
gemeinere Begriff der Stetigkeit geklärt wird. Hier- 
auf erfolet in Teil II und III der Aufbau der „Pro- 
jektiven Geometrie“, des Inbegriffes der Lageeigen- 
schaften geometrischer Gebilde. Unter Il werden 
die Sätze der Verknüpfung entwickelt und die 
grundlegenden Beziehungen der Sätze von Desar- 
Pascal aufgedeckt: aber erst die 


eues und von 
unter III eingeführten Sätze der Anordnung im 


Verein mit dem Begriff der Stetirkeit und seinen 
Abarten gestatten die vollständige Begründung der 
projektiven Geometrie. Als mächtiges Werkzeug 
dieser Untersuchungen wird zugleich die Koordi- 
natengeometrie begründet und durch die „Wurf- 
rechnung“ vollendet: sie gewährt auch die Mög- 
lichkeit, die imaginären Elemente 
folgerichtig einzuführen. In Teil IV fließt aus der 
Einsetzung Begriffes der .uneigentlichen 
Punkte“, der unterstellten Schnittpunkte sich nicht 
schneidender Geraden einer Ebene. die .„Affine Geo- 
metrie“, der Inbegriff der mit dem Parallelismus 
zusammenhängenden Eigenschaften, und in Teil V 
aus dem Begriff der Strecke und «den an ihn an- 
schließenden Grundsätzen der Kongruenz die .Me- 
trische Geometrie“. Hier nın treten nach der Zahl 
der uneigentlichen Punkte auf einer Geraden oder 

was unter gewissen Bedingungen auf dasselbe 


des 


der Geometrie ' 


hinauskommt — nach dem Verhalten der Winkel- 
summe im Dreieck geren zwei Rechte drei logisch 
sleichberechtigte Annahmen auf. aus denen das 
uns gewohnte Lehrgebäude der Euklidischen Geo- 
metrie und die beiden Niecht-Euklidischen Geo- 
metrien — wie vorher unter umfassender Prüfung 
der Abhängigkeiten zwischen den einzelnen 
Sätzen — entwickelt werden. 

Dresden. W. Ludwig. 235 


Dr.-Ing. R. HUGERSHOFF, Prof. a. d. Techn. 
Hochschule Dresden, Ausgleichungsrech- 
nung, Kollektivmaßlehre und Korre- 
lationsrechnungimDienste von Tech- 
nik. Wissenschaft und Wirtschaft. 
(Sammlung Wichmann, Bd. 10.) VI+86 S. Berlin 
1940, Verlag Wichmann. Preis geb. 6 M. 

Das vorliegende kleine Buch, das insbesondere 
für Forst- und Kulturingenieure wie auch für 
Naturwissenschaftler geschrieben ist, gibt zunächst 
als rechentechnische Grundlage der Kollektivmaß- 
lehre eine kurze Einführung in die Ausgleichungs- 
rechnung nach der Methode «der kleinsten Quadrate. 
Hier wird u. a. auch die Abgleiehung einer Beob- 
achtungsreihe durch eine ganze rationale Funktion 


gereben. Der zweite Teil bringt dann Mittelwerte 
und Streuungsmaße eines Kollektivs, normale, 


asymmetrische und mehrgipfeligee Häufigkeitsver- 
teilungen und den Vergleich mehrerer Kollektive. 
Der dritte Teil schließlich, die Korrelations- 
rechnung, benutzt im wesentlichen den Korrelations- 
koeffizienten, behandelt aber auch kurz die Korre 
lation von Zeitreihen. Die zahlreichen Beispiele 
sind vor allem für den Studierenden der Forst 
wissenschaften von Interesse. 

Das Buch scheint mir als erste Einführung in das 
behandelte Gebiet recht geeignet zu sein. 


Dresden. Willers. 249 
Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 


bleibt vorbehalten): 


Dr. HANS GEORG MÖLLER, o. Prof. a. d. Univ. 
Hamburg, Grundlagen und mathematische 
Hilfsmittel der Hochfrequenztechnik. 
(Lehrbuch der drahtlosen Nachrichtentechnik Bd. 1. 


herausgegeb. von Nicolai v. Korshenewsky und 
Wilhelm T. Runge. XVI+293 S. m. 353 Abb. 
Berlin 1940. Verlag ‚Julius Springer. Preis «eh. 


28.80 M. 


NACHRICHTEN 


Neu berechnete Tafeln. 

J. B. Russel: Journ. of math. and phys. (Massa- 
ehusetts). Bd. XII (1932/3) S. 291 bis 297: A table 
of Hermite functions. Werte der Funktionen des 
parabolischen Zylinders 


Yun (2) =e”*lz H„(z) 
mit 5 geltenden Stellen für n„=012...11 und 


2 =—=0,04... 1,00 (Schritt 0,04), = 1,1... 4.0 (Schritt 
01,2 =42...T0 (Schritt 02), z= 75, = 

Dabei sind die Hermiteschen Polynome defi- 
niert durch 


an | en 
. +2xt 
H„(s) = e ), . 


dam\ 


(22)" — 2(s) (2rjn—2 


+ 1-3. 22 4) rn 4 


Mit den Bezeichnungen bei Jahnke-Emde.3. Aufl. 
1938, S. 32, ist 
Pu )=Ynt2n Y2a'wn(ey2). 

(Für die g„ mit großem n ist der Schritt zum 
Teil sehr groß. Die sechsten Differenzen sind noch 
dreiziffrig und weisen schon Vorzeichen-Schwan- 
kungen auf.| 

Stuttgart. Fritz Emde. 263 


Persönliches. 

Dr. phil. Ludwig Föpp| o. Prof. für Festigkeits- 
lehre an der Techn. Hochschule München wurde zum 
ordentlichen Mitglied der Bayrischen Akademie der 
Wissenschaften ernannt. 


Das von Prof. W. Müller geleitete Institut für 
Theoretische Physik an der Universität München 
wurde in „Institut für Theoretische Physik und 
Angewandte Mechanik“ umbenannt. 
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